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ATHABASCA FORMATION OF WESTERN CANADA 
By WILttAm CARRUTHERS Gussow 


ABSTRACT 


Geological reports over the past 50 years have so confused the age designation of the 
Athabasca formation that logic has been abandoned even in the face of new evidence 
that the current Precambrian designation is untenable. 

The absolute age of the Athabasca formation cannot be fixed without fossil evidence. 
Herein, it is tentatively assigned to the Middle Devonian? or Late Cambrian? on the 
basis of structural relationships and lithologic characteristics, with full awareness of the 
risk involved in such a practice. The name Martin Lake series is here introduced for 
the ancient folded and faulted sedimentary-volcanic rocks of the Martin Lake area, 
which have been erroneously correlated with the Athabasca formation. Recent dating 
indicates that they are late Early Precambrian and are much older than the ‘‘Late 
Proterozoic” to which they have been assigned. 

This is an example of the need for a revision of the Precambrian time classification and 
the abandonment of terms such as Proterozoic, Algonkian, and Archean, which are 
carry-overs from the old concept of a twofold subdivision of the Precambrian—now 
known to be at least five to six times as long as the fossil record. 
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INTRODUCTION it as the “Athabasca Sandstone” and desig- 


nated it as Cambrian. Dowling and Tyrrell (77 
Tyrrell, 1897) concurred with McConnell’s age 
classification but correlated the Athabasca 
RBS sandstone with the Keweenawan, which at that 
’ Athabasca i is both the original ‘and the recently time was also considered to be Cambrian. In 
revised spelling; Athabaska is now obsolete. 1914, McInnes placed the Athabasca sandstone 


The extensive flat-lying sandstone along the 
south shore of Lake Athabasca! (Fig. 1) was first 
described by neem (1893). He referred to 
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in the Upper Huronian. However, in 1915, 
when his map was published, the Athabasca 
sandstone indicated 
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Reasons ior this Precambrian classification are 
the apparent absence of fossils from the Atha- 
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FiGureE 1. 


compiled from published and unpublished sources. 


and considered to be Precambrian(?). Camsell 
(1915) and Alcock (1915; 1917) followed Mc- 
Innes and placed the Athabasca sandstone in 
the Keweenawan(?). Since then it has been 
referred to the Precambrian, not, however, with- 
out question. J. B. Mawdsley (1949) and Blake 
(1956) attempted to determine its age but were 
unable to find fossils, so it retained its Pre- 
cambrian designation. 


INDEX Map SHowinG LAKE ATHABASCA AREA 


Regional distribution of the Athabasca formation and the Cambrian, Ordovician, and Middle Devonian 
basal sandstones is indicated. Map base traced from Geological Survey of Canada Map 1045A. Geology 





flat-lying sandstones with a thick series of 
tightly folded and deformed Precambrian sand- 
stones and volcanic rocks in the area of Beaver-f 
lodge Lake, and the correlation of these withf 
“Proterozoic” rocks at Great Slave Lake, 
Dubawnt Lake, and Baker Lake (Fig. 1). To 
avoid further ambiguity, the Precambrian 
conglomerate-sandstone-volcanic rock series in 
the Beaverlodge area is here renamed the 
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Martin Lake series after the type locality at 
Martin Lake. The name Athabasca formation is 
retained for the flat-lying strata south of Lake 
Athabasca, which were originally named the 
“Athabasca Sandstone” by McConnell. 


preciation is also expressed to R. T. D. Wicken- 
den, C. K. Bell, B. C. MacDonald, J. S. Ker- 
meen, and K. D. Hannigan for their interest. 
Grateful acknowledgment is made to J. B. 
Mawdsley for arranging for the loan of aerial 
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Figure 2.—Map oF THE East ENp oF LAKE ATHABASCA AND Fonp pu Lac 
Map shows the contrasting topography of the folded Precambrian of the north shore and the flat-lying 
\ Athabasca formation of the south shore. Traced from Geological Survey of Canada Maps 364A and 365A 
\ (Alcock, 1936). 








ow |] 
95° 
The present author questions the Pre- photographs of the Trout Lake area. J. B. 
evonian? Cambrian age of the Athabasca formation and Mawdsley and P. S. Warren kindly read the 
Geology, attempts to re-evaluate the published data. The manuscript and offered valuable sugges- 
results of this study were presented in the dis- tions. Editorial comments by F. J. Alcock 
cussion of a paper by I. C. Brown and G. M. and J. C. Sproule are gratefully acknowledged. 
ries off Wright in a symposium on “The Proterozoic in 
n sand-— Canada” (Gill, 1957). A short preliminary out- REGIONAL SETTING 
Beaver-f line by Gussow was published in 1957. 
se with The original Athabasca sandstone area is 
Lake, ACKNOWLEDGMENTS bounded on the north by Lake Athabasca, Fond 
1). To du Lac, Fond-du-Lac River, and Black Lake. 
mbrian The writer expresses his sincere thanks to D. ~The first base maps of the region made from 


A. W. Blake and Mrs. Blake for making avail- 
able unpublished photographs, original field 
notes, and thin sections of rock samples. Ap- 


acrial photographs show that the sandstones of 
the south shore are essentially  flat-lying, 
whereas the topography of the north shore indi- 
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REGIONAL 


cates the highly folded and contorted nature of 
the Precambrian rocks there (Fig. 2). This 
sharply contrasting topography indicates a 
major angular discordance between the flat- 
lying rocks of the Athabasca formation and the 
highly folded, contorted, thrust-faulted, and 
metamorphosed rocks of the Churchill 
georogen.? Alcock’s larger-scale maps (1936) 
further emphasize this striking difference. The 
absence of major thrust faults from the highly 
contorted Precambrian area—an area which 
seems to have undergone mountain building— 
induced the author in 1936 to re-study these 
maps, with the result that a well-defined north- 
east fault, the Black Bay-Tazin River master 
fault (northwest side up) along the southeast 
shore of Liack Bay and extending northeast 
bevond the Tazin Lake Sheet (Map 363A) and 
across the northwest corner of the Fond-du-Lac 
Sheet (Map 364A) became apparent. In addi- 
tion, the east-west Tazin River fault, the 
Mickey Lake (Tom) fault, and the St. Mary’s 
Channel (Gunnar) fault across the Cracking- 
stone Peninsula were recognized. In the eastern 
part of the north-shore region, the Grease River- 
Chambeuil Lake and the Black Lake master 
faults were postulated (Fig. 3). 

These observations indicate that all the for- 
mations of the north shore, including the young- 
est so-called Athabasca series of the Martin 
Lake and Tazin Lake areas are of pre-fault age. 

No structural trends are apparent in the flat- 
lying formations of the south shore. Along the 
narrow east arm of Lake Athabasca, Pine 
Channel, Fond du Lac, and Fond-du-Lac 
River, the deformed rocks of the north shore 
terminate abruptly against the flav-lying strata 
of the south shore (Fig. 2). This angular dis- 
cordance represents a profound unconformity 
and indicates that the Athabasca formation of 
the south shore is post-faulting and post-de- 
formation and very much younger than the 
contorted Precambrian formations of the north 
shore. The Tectonic Map of Canada (Derry, 
1950)? and the new Geologic Map of Canada 


* Rice (1949, p. 283) defines an orogen as “...a 
regional unit involved in mountain-making.” A 
georogen is considered herein to be a tectonic unit 
of continental proportions. 

3 The Tectonic Map of Canada shows a series of 
northeast-trending faults in the Athabasca forma- 
tion south of Black Lake. There is no known evi- 
dence for these, and nothing in the topography sug- 
gests their presence. Furthermore, it is stated on 
the map that the sandstones in this area are “... 
intruded by basic sills.” There is no justification for 
this statement, and this appears to be an erroneous 
belief that has become firmly established. 


uw 


SETTING 


(1955) show the same unconformable relation- 
ships along the southern contact of the 
Athabasca sandstone area. 


INVESTIGATIONS SINCE 1917 


Alcock (1936) carried out a detailed recon- 
naissance of the north shore and initiated a 
geological mapping program which has con- 
tinued ever since. He named and subdivided 
the Precambrian formations and established a 
pattern for the work that followed. The 
Athabasca formation was correlated with the 
Martin Lake series, although evidence to the 
contrary had been published at the time. Both 
were placed in the Proterozoic and assigned to 
the Keweenawan. 

Sproule, mapping the southern margin of the 
Athabasca formation in the Cree Lake area 
(Fig. 3), described it as ‘‘probably of Palaeozoic 
age” (1938b), and in 1941 he described it as 
“early Palaeozoic age” on preliminary maps, 
but on his final maps it is designated as Late 
Precambrian, Proterozoic (?)*, and the marginal 
notes state it may be Paleozoic. Tyrrell’s (1897) 
diabase dike is again shown on the final maps 
as cutting the Athabasca sandstone on Diabase 
Peninsula, on the west side of Cree Lake (but 
with a north-south strike), although this was 
not shown on Sproule’s preliminary map. The 
later Geological Survey maps of Saskatchewan 
(1947) and Alberta (1951) retained the Pre- 
cambrian designation for the Athabasca forma- 
tion and assigned it to the Proterozoic. The 
Geological Map of Canada (1955) assigned it 
to the Late Proterozoic. 

Furnival (1940; 1941a; 1941b), mapping in 
the extreme northeast corner of the Athabasca 
sandstone area followed Alcock (1936) in cor- 
relation and dating. However, he did not 
mention any basalt flows or dikes in the 
Athabasca formation. 

During the past 8 years, a new series of 
1-mile maps has been published by the Geologi- 
cal Survey of Canada. Christie (1949; 1953) 
mapped both the Athabasca formation and the 
Martin Lake series as Precambrian, although 
abundant structural evidence casts serious 
doubt on the validity of this assumption. 

J. B. Mawdsley (1949; 1956) searched for 
fossils in the sandstone in the vicinity of Pine 
Channel (Fig. 2) without success. He observed 
that “...it rests with profound unconformity 


4Sproule (personal communication, 1958) ad- 
vised the writer that this was an editorial change 
made without his concurrence. 
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on a surface of ancient, probably Archaean 
rocks” (1949, p. 22). 

Age determinations (Collins and Freeman, 
1951) on uranium minerals from veins which 
cut the Martin Lake series indicated that these 
rocks are older than 1680 million years®, and 
as a consequence the rocks of the type area of 
the Athabasca formation were dated as late 
Early Precambrian. 

Nevertheless Fraser (1954) placed the 
Athabasca in the Proterozoic. He differentiated 
as to lithology but not as to age, although rocks 
cf two definite ages can be easily separated 
on his maps. 

Hale, in his Black Bay Map (1954a), also 
placed the “Athabasca Series” in the Pro- 
terozoic but dated it as pre-diabase. The rocks 
mapped as ‘“‘Athabasca Series” are undoubtedly 
equivalents of the Martin Lake series and 
accordingly are Early Precambrian. No true 
Athabasca formation is present in the Black 
Bay map area. Hale (195-+4b; 1955) differentiated 
between an Upper and Lower part in the 
“Athabasca Series’; these are separated by 
an angular unconformity. Although he noted 
the fresh appearance of the Upper Part, he 
placed both in the Late Precambrian (Pro- 
terozoic). Now, however, he indicates that the 
age of the diabase relative to the Upper Part 
is uncertain (1954b). From his maps and de- 
scriptions, it is obvious that both the Martin 
Lake series and the Athabasca formation are 
represented. 

J. C. Mawdsley (1954, B.Sc. thesis, Univ. 
Alberta) made a mineralogic study of the 
Athabasca sandstone and observed that the 
Martin Lake series has been correlated with the 
Athabasca formation 


*.,.in spite of the lack of any lithologic basis for 
such a correlation. The two types of rock were 
quite obviously deposited under very different en- 
vironments.” 


Mawdsley found a remarkably small per- 
centage of heavy minerals in the Athabasca 
sandstone and identified minor amounts of 
tourmaline, zircon, and rutile, and traces of 
magnetite and garnet. On the basis of optically 
oriented authigenic tourmaline, he concluded 
that the Athabasca sandstone must be marine. 

In 1952, Blake undertook a reconnaissance 


5 More recent work by Eckelmann and Kulp 
(1956) concludes that the time of initial uranium 
deposition was about 1900 + 20 m.y. ago. The 
validity of absolute age based on pitchblende, how- 
ever, is apparently still open to question. 
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of the whole areca. He made the startling dis. 
covery of, and described, the Trout Lake 
dolomite member, 40 miles south of Lake 
Athabasca (Fig. 3). In addition, he recorded 
valuable data and collected many specimens, 
Ill health prevented publication of his pre. 
lirniaary report until late in 1956; it represents 
an ‘nterpretation of his field notes. Blake came 
» clear understanding of the age of 
this interesting formation, and given an oppor- 
tunity to thoroughly review the literature in 
the light of his field discoveries, he probably 
would have solved the many difficult field 
problems. With the help of Blake’s field notes 
and photographs of outcrops, the loan of thin 
sections, and the present author’s own field 
experience, the author is taking the liberty of 
reinterpreting the evidence which Blake 
collected. 

Blake described “folding” at four widely 
separated localities. The dips east of Turnor 
Point, and south of Riou Lake, and along 
Fond-du-Lac River east of Black Lake are 
probably cross-bedding and foreset beds. Alcock 
(1936, Pl. VII B) shows a similar occurrence 
south of Fond-du-Lac River. The locality de- 
scribed by Blake on Fir Island, however, sug- 
gests that these quartzites belong to the base- 
ment rocks and are Tazin series? Blake states 
that ‘‘...in no other part of the sandstone 
area have dips of more than a few degrees been 
observed ...”’ (1956, p. 8). He also states 
that ‘‘. .. no faults have been noted in the main 
area of the Athabasca formation” (1956, p. 8). 

Hale (1954b) mapped a series of dips in the 
Athabasca (?) formation on Charlot Peninsula 
and along the shore south of the mouth of 
Charlot River. In all probability, these repre- 
sent depositional dips or are foreset beds. At 
all places the Athabasca formation appears to 
be flat-lying or with low depositional dips. 

Blake observed that 


clases 


**.. the rocks north and south of the lake [Atha- 
basca] differ lithologically in many respects and 
have apparently formed under different conditions. 
Rocks south of the lake are composed almost en- 
tirely of well-rounded quartz grains. With few ex- 
ceptions, conglomerate beds and feldspar grains are 
virtually absent and the red color so prevalent 
north of the lake is lacking. These facts indicate 
that the sediments have been transported a great 
distance and deposited on a broad, nearly level 
surface, and suggest that the correlation is unwar- 
ranted” (1956, p. 9). 


In spite of this, he attributed the difference 
to a transition from the arkosic type north of 
the lake to the feldspar-free sandstones south of 
the lake. 
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Alcock recently made the same observation 
(personal communication, 1956) and con- 
cluded that there is no reason for correlating the 
two rock groups. 

During 1952, 1953, 1954, and 1955, Tremblay 
(1955; 1957a; 1957b) mapped the Martin Lake 
type area and vicinity. On four of his five 
preliminary map sheets (scale, 1 inch to 800 
feet), the Martin Lake series is referred to as 
the “Athabasca Series’ and is called Pro- 
terozoic. 

A swarm of gabbro dikes and sills which cut 
the pre-Martin Lake complex northwest of the 
Black Bay-Tazin River master fault is shown 
as post-Martin Lake (post-Athabasca by 
Tremblay’s terminology) in age because at two 
places in the southwest corner of the map area, 
gabbro cuts rocks of the Martin Lake series or 
is intruded between. them and the underlying 
Tazin rocks. Trembiay admits that the “late” 
gabbros also resemble a series of sheets, sills, 
and dikes intrusive into the lower part of the 
Martin Lake series, below the series of inter- 
bedded extrusive lavas and arkosic sediments. 
There is thus no evidence that any of the 
igneous rocks are post-Martin Lake in age. 
They undoubtedly welled up through fractures 
in the pre-Martin Lake basement rocks during 
early Martin Lake time and were intruded into 
the lower Martin Lake sediments and, later, 
poured out as extrusive lavas interbedded with 
arkosic sediments. They are all here considered 
early Martin Lake in age. Tremblay’s basic 
information was used in compiling Figure 4 
but has been reinterpreted with the help of 
aerial photographs. 

S. J. T. Kirkland (1953, unpub. manuscript, 
Queen’s Univ.) concluded that 


“Apart from the clastic nature of the rocks in the 
two areas (north of Lake Athabasca and south of 
Lake Athabasca), there is little ground for correla- 
tion The so-called Athabasca (Martin Lake) rocks 
seen by the writer (Kirkland) at the Eldorado Ace 
mine, near Beaverlodge Lake, bear little resemblance 
to those in the Middle Lake area ... They are of a 
decidedly different nature ... and the correlation of 
such rocks with the type Athabasca series is ques- 
tionable.”’ 


In the summer of 1956, Fahrig (1958)* 
mapped the eastern boundary of the Athabasca 
area and designated it as Proterozoic. The same 
year, C. K. Bell carried on detailed geological 
mapping in the Crackingstone Peninsula; the 
results of Bell’s investigations have not been 
published and are not available to the author. 





§ Published after completion of this manuscript. 


The basal member of the Athabasca forma- 
tion, because it resembles the arkosic sandstones 
and conglomerates of the Martin Lake series, 
has led geologists (Alcock, 1920; 1936) working 
on the north shore to correlate these two on the 
basis of lithology. This correlation has become 
so conventional that even though Hale (1954b) 
discovered a profound angular unconformity 
within his “Athabasca Series”, he placed both 
parts in the Precambrian. Christie (1949; 1953), 
mapping in the Beaverlodge Lake area, con- 
fidently correlated the Martin Lake arkose and 
conglomerate with the basal sandstone and 
conglomerate of the Athabasca formation on 
the islands south of Crackingstone Peninsula. 


MartTIN LAKE SERIES AND ATHABASCA 
FORMATION REDEFINED 


Martin Lake Series 


The type area of the Martin Lake series is 
defined as the vicinity of Martin Lake in the 
Beaverlodge Lake area, Saskatchewan, and 
extending northeast and southwest along the 
Black Bay-Tazin River master fault for about 
18 miles (Fig. 4). Structurally, the main part 
of the area is synclinal, pitching northeast, with 
dips on both limbs averaging 65°. The total 
preserved thickness of this series is estimated 
as 8000 feet and consists of siliceous arkosic 
sandstones and conglomerates, interbedded 
with basaltic lava flows which make up ap- 
proximately 1000 feet of the section. The 
original thickness is unknown because of ero- 
sion. The whole series is offset by numerous 
faults’ and is cut by veins containing pitch- 
blende. 

Three other areas in which the Martin Lake 
series crops out are (Fig. 5): (1) a synclinal area 
at the east end of Tazin Lake, northwest of the 
type area; (2) the Thluicho-Gulo lakes basin, 
north of Camsell Portage; (3) the Charlot River 
basin, between Ellis Bay and Wellington Lake. 
A remnant occurs along the southeast shore of 
Black Bay against the Black Bay-Tazin River 
master fault at Langley Bay, and another 
remnant occurs on Slate Island and Peninsula. 
These formations are undoubtedly late Early 





7 A major fault known as the St. Louis fault cuts 
the main syncline. It is easily traced on the de- 
tailed maps of the area (Christie, 1949; 1953; 
Tremblay, 1957a; 1957b) and on aerial photo- 
graphs. The fault trends east-west and passes 
through the narrow cleft used by the portage and 
road at the northern end of Beaverlodge Lake. 
Movement on this fault is north side up, relative 
to the south side, offsetting the synclinal axis and 
the formational contacts. 
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Precambrian as shown by age dating in the 
type area (Collins and Freeman, 1951; Eckel- 
mann and Kulp, 1956), and can probably be 
correlated with the sandstone-volcanic series 
at Dubawnt Lake, Baker Lake, and Great 
Slave Lake (Fig. 1). 

Most of the rocks are of continental fresh- 
water origin and consist predominantly of 
arkoses to feldspathic sandstones, buff to 
red, with angular to rounded grains, coarse 
angular conglomerates, and associated volcanic 
flows and dikes. The finer-grained beds are 
cross-bedded, ripple-marked and exhibit sun 
cracks. 

The Martin Lake series lies unconformably 
on the deeply eroded Tazin group and the 
associated post-Tazin granitic rocks, which 
underlie most of the north-shore area. The 
Martin Lake series has been folded and faulted 
and deeply incised prior to deposition of the 
Athabasca formation. The folded Martin Lake 
series and older rocks are cut by many small, 
black, dark-green, or brown basic dykes. These 
are probably early Martin Lake and _pre- 
volcanic in age. 


Athabasca Formation 


The flat-lying Athabasca formation (Fig. 3) 
underlies an area of about 28,000 square miles’; 
it extends 260 miles from the Athabasca River® 
cast along the south shore of Lake Athabasca, 
Fond du Lac, Fond-du-Lac River, and Black 
Lake, to the vicinity of Wollaston Lake, and 
south for 135 miles to Cree Lake. The south- 
eastern boundary has not been defined and may 
lie somewhat farther west. The total thickness 
of the Athabasca formation is unknown; a 
vertical section 400 feet thick is exposed in the 
escarpment south of Lake Athabasca. The 
southern contact at Cree Lake lies nearly 900 
feet above the northern contact on Lake 
Athabasca. The flatness of the regional dips" 
may indicate a total onlapping stratigraphic 
section of nearly 1000 feet. The Athabasca 
formation crops out on several islands (Long, 


5 The total area of New Brunswick. 

®An unmapped exposure was observed by the 
author near the mouth of the Athabasca River on 
the main navigation channel. 

‘0 Dowling (in Tyrrell, 1897) believed that the 
sandstone beds dip slightly to the west; more re- 
cent workers note very gentle south dips. Possibly 
there is a regional southwest tilt similar to that 
recognized in the Plains; if so, then older beds 
should be exposed in the east, with progressively 
younger beds exposed westward along the present 
erosion surface. 


juhnston, Stewart, Grouse, Feather, Course. 
Bigfowl, and numerous unnamed - smaller 
islands) near the north shore of Lake Athabasca 
(Fig. 5), and some isolated patches or outliers 
occur on the north shore—notably, a large 
area north and south of Spring Point, on Slate 
Island, Charlot Island and Peninsula, as well 
as in several places described by Alcock as 
“sandstone dikes’”"—showing that this forma. 
tion underlies most of Lake Athabasca and that 
it once extended even farther north. 

The Athabasca formation is almost entirely 
sandstone, commonly siliceous. A dolomite 
member, recently discovered at Trout Lake, 
about 40 miles south of William Point, is de. 
scribed separately. Typically, grains are 
rounded to subangular. The rounded quartz 
grains have optically oriented authigenic over- 
growths of quartz, many with sutured bound- 
aries. The sandstone is fine- to medium-grained, 
although grit and pebble bands occur; these 
mark bedding planes that are otherwise difficult 
to observe. The weathered surface of the rock 
is white, but locally it is buff. The rock is 
massive-bedded to thin-bedded, and _ cross- 
bedding is common. In the lower part of the 
section, the sandstone is iron-stained and 
ripple marks are common. Locally, red shale 
layers up to 10 inches thick occur near the base 
of the formation. In describing the Athabasca 
formation along the southern margin, Sproule 
(1941) mentioned numerous, irregular light- 
colored carbonate pellets up to 1 inch in di- 
ameter. A hand specimen from a boulder of 
Athabasca sandstone in the Mudjatik area 
contains numerous white, soft, kaolin chips. 

It has been suggested that the Athabasca 
sandstones may be aeolian, which would explain 
the absence of fossils. However, the authigenic 
tourmaline overgrowths, the very regular bed- 
ding over considerable distances, both hori 
zontally and vertically, the even cross-bedding, 
the clay chips and irregular carbonate pellets, 
and the Trout Lake dolomite member seem to 
show that the sandstone and the dolomite mem- 
ber are of marine origin. 

TROUT LAKE DOLOMITE: The author believes 
that the Trout Lake dolomite member is part 
of the Athabasca formation. 

Blake (1956, p. 5) noted high ridges of dolo- 
mite which crop out in the vicinity of Trout 
Lake, 40 miles south of William Point (Fig. 3) 


1 Sandstone preserved in cracks and fissures in 
older rocks. One of these, in the vicinity of Elliot 
Bay, has been traced more than 7000 feet (C. K 
Bell, personal communication, 1957). 
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and are clearly visible in aerial photographs. 
They extend at least 20 miles in an arcuate 
direction, east and west from Trout Lake 
(Fig. 6). They are now known to crop out in a 
complete circle and are 300-400 feet high (K. D. 
Hannigan, personal communication). 

Blake described the carbonates as exceedingly 
fine-grained and highly contorted. They may 
be, he noted, cream, buff, pink, or reddish; 
the laminae are extremely fine to very coarse. 
The dolomite is fine-grained or aphanitic. In 
places it is folded into major anticlines and 
synclines with dips up to 80°. In addition, the 
rocks have been deformed locally on a small 
scale. There are many zones of intense breccia- 
tion, intricate contortion, and rupture. 

Blake noted thin sections of the carbonate 
showing delicately laminated beds which are 
moderately brecciated. Rounded “‘floating”’ 
sand grains of less than 1 millimeter in diameter 
were found in two sections. The dolomite is 
unmetamorphosed. 

Blake’s description, together with an ex- 
amination of the thin sections, indicated to 
this author an evaporitic carbonate facies 
which has undergone soft-rock or penecon- 
temporaneous deformation, or ‘‘autoclastic 
brecciation”’, due to subaqueous slumping after 
partial lithification. There is no evidence of 
metamorphism. The rock is unknown in the 
Precambrian of the north shore. Similar rock 
types which are common to the Middle 
Devonian Elk Point formation of the Interior 
Plains and to the Ghost River formation of 
the Front Ranges of the Canadian Rockies 
will be discussed in following sections. 

Extremely fine microlaminae, micro-oolites, 
microstylolites, and primary orientation of 
minute, acicular crvstals normal to the laminae 
were observed. 

Blake (1956, p. 6) stated ‘‘... it is possible 
that the Trout Lake (carbonate) is a part of 
the Athabasca series. ...’’ The present writer 
concurs. It appears to be the youngest member 
of the Athabasca formation. 

Shrock (1948, Fig. 236, p. 276) illustrated 
that subaqueous flowage and deformation can 
cause sliding on slopes of less than 6°; this may 
explain the contortion of the sediments. The 
descriptions of Blake’s traverses (unpublished 
field notes, August 10, 1952) seem to fit such an 
explanation, and the contortion is interpreted 
as penecontemporancous deformation of partly 
consolidated calcareous oozes, causing folding 
and brecciation shortly after deposition of the 
sediments. 


Finally, the east-west strike of the Trout 
Lake outcrop area is sharply discordant with 
the prevailing northeast grain of the basement. 
This suggests it is not Precambrian. 


Sus-ATHABASCA UNCONFORMITY 


The basal contact oi the Athabasca formation 
is exposed at several places along the northerm 
edge. The basal conglomerate lies on the base. 
ment complex with a profound angular dis. 
cordance. This unconformity has been described 
at the following places: (1) on the islands south 
of Crackingstone Peninsula; (2) at several 
points between Stony Rapids and Middle Lake; 
(3) near ‘e southeast shore of Black Lake 
(J. S. Ker....en, personal communication); (4) 
along Pine Channel; (5) east of Poplar Point; 
(6) north of the mouth of Charlot River; and 
(7) on Slate Island (Figs. 2, 3, 5). 

Tyrrell (1897) assigned the underlying steeply 
dipping sedimentary rocks (Martin Lake 
series?) on Slate Island (Fig. 5) to the Huronian 
and placed the flat-lying strata (Athabasca 
formation) in the Cambrian. 

Alcock (1936) described the unconformity on 
several of the islands south of Crackingstone 


Peninsula and mentioned the ‘‘sandstone 
dikes.” 
Blake (1956, p. 8) described the contact 


“ 


as “...an angular unconformity of great 
hiatus.” East of Poplar Point, a flat-lying sand- 
stone rests on steeply dipping, ferruginous 
quartzite of the Tazin series. A few inches of 
quartzite and vein quartz conglomerate lies 
at the base of the sandstone. Between Stony 
and Middle lakes, the unconformity is exposed 
intermittently for 5 miles at the base of a 25- 
to 50-foot sandstone escarpment. There, too, a 
thin layer of quartz conglomerate appears at 
the base of the sandstone. The sandstone in this 
area is highly hematitized and cross-bedded 
and contains many grit and pebble bands. The 
Tazin-type rocks 8 feet below the unconformity 
were so weathered in pre-Athabasca time that 
their original character cannot be ascertained. 
The weathering products are almost entirely 
quartz and kaolin, with some iron oxide— 
“typical lateritic residual soils from long- 
continued weathering or decomposition of the 
underlying rocks” (S. J. T. Kirkland, 1953, 
unpub. manuscript, Queen’s Univ.). 

Hale (1955) described the unconformity at 
the west end of Slate Island and observed the 
fresh appearance of the upper part of the series, 
adding, 
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“Only negative evidence, the apparent absence of 
jossils, supports the assignment of this part of the 
series to the Precambrian.” 


Hale probably observed the _ flat-lying 
Athabasca formation (his Upper Part) lying 
unconformably on the Martin Lake series 
(his Lower Part) on Slate Island and north of 
the mouth of Charlot River. 

Along the southern margin of the Athabasca 
formation, numerous inliers of crystalline base- 
ment rocks and outliers of Athabasca formation 
have been mapped. No actual contacts have 
been observed, but Sproule (1941a; 1941b; 
1941c) described the basement rocks as rotted 
and disintegrated near the contact by long 
weathering. They contrast sharply with the 
freshly glaciated outcrops. 


DIABASE DIKES 


No diabase dikes are known to cut the 
Athabasca formation. A note on the Tectonic 
Map of Canada (Derry, 1950) states that the 
sandstone area is intruded by basic sills, but 
this statement is correct only for the lower 
Martin Lake series of the north shore and has 
been applied erroneously to the whole area. 

At one locality on the west shore of Cree 
Lake, near the southern margin of the sand- 
stone area, Tyrrell (1897) described an outcrop 
of light-green uralitic diabase striking N.65°E., 
which he believed cut the Athabasca 
sandstones. More probably this represents a 
small embayment in the southern edge or a 
window of basement rocks. The strike conforms 
to the regional strike in the basement complex, 
and the diabase was not actually found in 
contact with the Athabasca formation. The 
near-by quartzites could be Tazin group. 

Blake (1956, p. 10) mapped a small outcrop 
of diabase on the south shore of Lake 
Athabasca, approximately 4 miles east of 
Turnor Point. Its relationship to the Athabasca 
formation is not known. If it does cut the 
Athabasca sandstone, it would be the only 
known example. It could be a Precambrian 
inlier. 


CORRELATION OF THE ATHABASCA FORMATION 
General Statement 


The Martin Lake series is a thick series of 
late Early Precambrian interbedded continental 
arkosic sandstones, conglomerates, and basalt 


flows, which have been folded and meta- 


morphosed and are cut by numerous faults 
and quartz veins. The Athabasca formation 
is a relatively thin series of flat-lying, unmeta- 
morphosed, white to pink marine sandstones 
and minor quartz-pebble conglomerates, with 
an evaporitic carbonate member. It lies on the 
Precambrian basement complex with profound 
unconformity. To date, it has yielded no fossils. 

The Martin Lake sedimentary rocks are 
predominantly arkoses and feldspathic sand- 
stones which exhibit all the characteristics of 
rapid continental accumulation; the Athabasca 
sandstones, in marked contrast, are almost 
entirely pure quartz, the product of long and 
extensive chemical decomposition of an ancient 
land surface, followed by erosion and sorting 
of these residual soils. The end products consist 
of quartz, kaolin, and minor amounts of iron 
oxide, with only the most resistant of the ac- 
cessory minerals. 

Regionally, the Athabasca formation and 
the flat-lying Paleozoic formations of the 
Interior Plains are structurally conformable, 
and both lie with profound angular discordance 
on the crystalline rocks of the Churchill 
georogen (Fig, 7). Obviously, the Athabasca 
formation is much younger than even the 
youngest strata of the basement complex, the 
Martin Lake series, and is closer in age to the 
Paleozoic section of the Plains. 

Until diagnostic fossils are found, no precise 
dating is possible. Three possible correlations 
are: (1) correlation with the basal Cambrian 
sandstones in the subsurface of the Interior 
Plains, which are Middle and Upper Cambrian; 
(2) correlation with the basal Ordovician 
Winnipeg sandstone; or (3) correlation with 
the basal Middle Devonian sandstones, which 
overlie the Precambrian basement immediately 
west of the Athabasca area, in the northern 
half of Alberta (Fig. 1). 


Basal Cambrian 


In the subsurface of the Plains, the lower 
part of the Middle Cambrian consists of a 
poorly sorted, basal, white, pink, and red- 
mottled quartz sandstone with interbedded 
shales. The sandtones become pebbly and 
arkosic near the base and are glauconitic. They 
contain an abundant marine fauna. They grade 
into Upper Cambrian beds to the east by onlap. 

The Athabasca formation lacks the abundant 
marine fauna and the glauconite and mica 
which characterize the Middle Cambrian de- 
posits of the region. Furthermore, the Cambrian 
formations of the Plains thin eastward by 
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onlap and become more shaly. They have a 
zero depositional edge in the subsurface (with 
only slight pre-Ordovician erosion) and do not 
crop out along the margins of the Canadian 
Shield. 

Recent drilling southwest of Lac la Ronge 
near the north end of Montreal Lake pene- 
trated 650 feet of basal Cambrian sandstone 
fore passing into the Precambrian basement 


[his is most significant in de age ol 


\thabasca sandstone 


In Manitoba, basa ( crop out 
along the eastern edge ol \ 
zoic section where it borders the Canadian 
Shield. These, known as the Winnipeg forma- 
tion, conformably underlie fossiliferous Middle 
Ordovician rocks (Baillie, 1952). They lie on 
the Precambrian basement complex with pro- 
found unconformity. In the subsurface, these 
sandstones lie on rocks of Early Ordovician 
age. North of Lake Winnipeg, the Winnipeg 
sandstone strikes west into Saskatchewan and 
disappears beneath Middle Devonian rocks, 
which here lie unconformably on Ordovician 
rocks. 

In the vicinity of Great Bear Lake, and 
extending south to Great Slave Lake, sand- 
stones which conformably underlie — the 
Ordovician section and overlie the Precambrian 
complex with pronounced unconformity are 
probably age equivalents of the Winnipeg 
formation, although maps of the Geological 
Survey date them as Precambrian (Protero- 
zoic). The nearest outcrops of Winnipeg sand- 
stone in both areas are about 200 miles from 
the Athabasca sandstone area. 

Between Lac la Ronge and Great Slave 
Lake, Ordovician and Cambrian formations 
are absent because of pre-Devonian erosion, and 
Middle Devonian rocks lie directly on the 
granite. 

Wallace and McCartney (1928) studied the 
Winnipeg sandstone microscopically. The sands 
are well sorted, rounded, and pitted, but unlike 
the Athabasca sandstone, are cemented by 
calcite. Minor heavy accessories include, in 
decreasing order of abundance, zircon, magne- 
tite, tourmaline, staurolite, and rutile. This 
assemblage and the fossils in the conformably 
overlying strata indicate a marine environment. 

Lithologically, the Winnipeg sandstone is 
similar to the Athabasca sandstone in hand 
specimens, and they could be separated only 


runcated Paleo- 


by fossil evidence. The presence of evaporites 
in the Ordovician, however, is unknown. 


Basal Devonian 


Between Lac la Ronge and Great Slave 
Lake, along the truncated east edge of the 
Paleozoic and the margin of the Precambrian 
Shield, Middle Devonian strata lie directly on 
the Precambrian basement. Outcrops have 
been described at only a few places—at Lac la 
Ronge, on the Clearwater River, on the Firebag 
River, along the Slave River south of Fort 
Fitzgerald, at the junction of the Slave and 
Peace rivers, and on Great Slave Lake. 

McInnes (1910; 1914) and Dowling (in 
Tyrrell, 1897) described the Middle Devonian 
on the south shore of Lac la Ronge, Saskatche- 
wan. This is approximately 175 miles south and 
slightly east of Cree Lake, at the northern edge 
of the truncated Paleozoic. The rocks are fos- 
siliferous magnesian limestone or dolomite, 
sandstone, and conglomerate, lying directly on 
the Precambrian. The conglomerate contains 
limestone, pebbles which must be the remnants 
of eroded Ordovician about 25 miles to the 
south. This would date the basal conglomerate 
as post-Ordovidian. 

At Contact Rapids on the Clearwater River, 
Saskatchewan, at the truncated eastern edge 
of the Interior Plains, R. deWit (1952, un- 
pub. manuscript, Geological Survey of 
Canada), described the Paleozoic section im- 
mediately above the Precambrian. A_ thin 
regolith overlies the granite and is uncon- 
formably overlain by flat-lying sandstones of an 
entirely different character than the under- 
lying “‘wash.”’ No fossils have been found in the 
basal sandstone. The maximum thickness of 
the sandstone is 10.5 feet; beds are 4-12 inches 
thick. The sandstones grade upward into 
greenish-gray silts and sandy marls which 
contain fossil remains similar to the fossils 
in the overlying dolomite. The silts grade up- 
ward into sandy light-gray to buff, crypto- 
crystalline, dense dolomite and _ limestone 
known as the Methy dolomite which contains 
an abundant Middle Devonian"? fauna, gen- 
erally poorly preserved but rich in forms and 
containing some small brown spores. 

Throughout western Canada, the Devonian 
lies on older rocks with a sharp erosional un- 
conformity. Thus the whole paleozoic section 


" The first fossil collections from this area were 
tentatively assigned to the Silurian, but more com- 
plete material has established them as Middle 
Devonian. 
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until fossils are found in the Athabasca forma- 
tion, an correlation is impossible. 
Even in the absence of fossils, the presence of 


on the Clearwater River must belong to the 
Middle Devonian, as no unconformity is evi- 
dent within the sedimentary section; the only 


absolute 


TABLE 1.—TENTATIVE REVISED TABLE OF FORMATIONS FOR NORTHERN SASKATCHEWAN 


Geologic ‘Time Formation Lithology Remarks 
PALEOZOIC 
Middle 
and/or Upper 
Cambrian? 


than 500 mil 
lion years 


Devonian? Athabasca Forma- | Pure white quartz sandstones and — Less 


tion occasional quartz-pebble beds, 
with a dense evaporitic dolomite 


member 


Profound Unconformily 


PRECAMBRIAN 


Early Precambrian Greater than 2000 


million years 


Martin Lake series Interbedded arkose, feldspathic 
sandstone and conglomerate, and 
basaltic lava flows 

Some intrusive sheets, sills and dikes 


of early Martin Lake age 
Major Unconformity 


Post-Tazin 
intrusives 


Lamprophyre dikes 
Granite and gneiss 


Intrusive Contact 


(Archean?) Tazin group Metamorphosed sedimentary and 

volcanic rocks 
Conglomerate 
Dolomite and dolomitic quartzite 
Ferruginous quartzite 
Quartzite 

Note that the writer does not subscribe to the term ‘Proterozoic’? which is commonly intended to imply 
“immediate pre-Cambrian age’. Obviously, none of the formations of Northern Saskatchewan fits this 
definition. The term Archean is believed to serve a useful purpose to designate the ancient rocks of the pri- 
mordial crust. Howevei, it is not certain if the Tazin group is Archean or somewhat younger. 





unconformity is that at the base of the sand- 
stone. The whole conformable section may be 
assigned to the Middle Devonian and is cor- 
related with the upper part of the Elk Point 
of the Interior Plains. 

The Clearwater section is about 50 miles 
south of the Athabasca sandstone area. Sproule 
(1951) correlates this sandstone with the 
Athabasca sandstone. R. deWit (personal com- 
munication) has compared hand specimens of 
the Athabasca sandstone and_ the _ basal 
Devonian sandstone of the Clearwater River 
and agrees that they seem identical. However, 


an evaporitic member in the Athabasca formaf 
tion supports such a correlation. The type of 
environment postulated for the formation o. 
the authigenic tourmaline described by J. C. 
Mawdsley (1954, unpub. B.Sc. thesis, Univ 
Alberta) in the Athabasca formation corre- 
sponds with that of the Middle Devonian of 
the Interior Plains. 

The Slave River exposures are small, but the 
rocks are similar to those on the Clearwater, 
and the contact with the Precambrian is identi- 
cal. On Great Slave Lake, the Middle Devonian 
s well exposed and abundantly fossiliferous. 
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It is subdivided into the Presqu’ile and Pine 
Point,"= and is correlated with the Middle 
Devonian Elk Point section of the Plains. The 
base is exposed only at Gypsum Point on the 
North shore, where Silurian? and Ordovician 
rocks separate the Devonian from the Pre- 
cambrian. North of Lesser Slave Lake, through- 
out the northern half of Alberta, the Middle 
Devonian immediately overlies the Precam- 
brian. Locally developed basal sandstones have 
a maximum thickness of 60 to 100 feet. 


CONCLUSIONS 


Evidence favors a_ tentative Middle 
Devonian? or Upper Cambrian? age for the 
Athabasca formation. It is also possible that 
the sandstone is Cambrian and the Trout 
Lake dolomite is Middle Devonian. The Martin 
Lake series is late Early Precambrian and 
much older than the flat-lying Athabasca 
formation. 

The foregoing deductions should be investi- 
gated further, and an attempt should be made 
to find evidence of organic remains such as 
fossil fish fragments (which have recently been 
found in the Ghost River formation of the 
Front Ranges and in the basal Devonian sand- 
stones of the Interior Plains), Chara, and plant 
remains. 

The Athabasca formation is either flat-lying 
or with slight depositional dips; the Martin 
Lake series is steeply dipping to overturned. 
Pebbles in the Athabasca are fresh-looking and 
undeformed, whereas in the Martin Lake series 
and older formations, they may be gneissic. In 
the vicinity of the Martin Lake series, the 
basal conglomerate of the Athabasca formation 
contains pebbles of the Martin Lake series. 





8 Although E. M. Kindle (1942) places the over- 
lying Slave Point formation in the Upper Devonian, 
until recently (Law, 1955), it has generally been 
regarded as Middle Devonian. It was originally 
assigned to the Middle Devonian by Cameron 
(1922). No diagnostic fauna has been collected 
from the type area on Great Slave Lake. However, 
new fossil collections (W. G. Crook, personal com- 
munication, February 1957) from Slave Point 
equivalent (?) in the subsurface, and identified by 
P. S. Warren, have established its age as Late 
Devonian and correlate it with the lower part of the 
Waterways and the Beaverhill Lake formations. 

Fish remains from thin shales interbedded with 
the basal Middle Devonian granite wash at the Red 
Earth oilfield, Alberta, have now heen identified as 
Coccosteus by Robert H. Dennison (G. FE. Meade, 
personal communication, October, 1958). The same 
forms have been found in the Ghost River formation 
where it crops out on End Mountain, in the Front 
Range of the Rocky Mountains, Alberta. 


The Martin Lake sandstones and conglomerates 
are commonly associated with volcanic flows 
and have been faulted and folded and cut by 
numerous basic dikes; no known dikes or vol- 
canic rocks are associated with the Athabasca 
formation. The Martin Lake series is of con- 
tinental origin, the Athabasca is marine. 

A tentative revised table of formations for 
the geology of northern Saskatchewan is given 
in Table 1. 


Economic GEOLOGY 


The possibility of uranium mineralization 
in the Athabasca formation should be investi- 
gated. The Beaverlodge area testifies to the 
presence of rich primary uranium mineraliza- 
tion in the basement rocks, and possibly the 
basal sandstones and conglomerates of the 
Athbasca formation are the host for Blind 
River-type mineralization. 

Low-grade uranium mineralization has al- 
ready been discovered in the Athabasca forma- 
tion on Middle Lake (Blake, 1956; Hriskevich, 
1949) and on Stewart Island, and it should be 
remembered that the absence of ore-grade 
values at the outcrop does not preclude the 
possibility of ore-grade mineralization at depth, 
beyond the zone of surface leaching. 
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SIMPLE GEOLOGIC STRUCTURES 


By ALLAN R. SANFORD 


ABSTRACT 


An analytical method is derived which gives the elastic response of a homogeneous 
rock layer to two-dimensional distributions of vertical displacement applied along its 


lower boundary. 


Displacement field, stress distribution, and distortional strain-energy 


density diagrams are presented for three types of displacement applied at the lower 
boundary of 5-km thick layers possessing average sedimentary rock properties. These 
three types of displacement are: (1) sinusoidal vertical displacement and no horizontal 


displacement; 
displacement; 


(2) an approximate step in vertical displacement and no horizontal 
and (3) sinusoidal vertical and horizontal displacement (horizontal 


displacement 90° out of phase with the vertical). Displacement fields and stress dis- 
tributions for each type of applied displacement are nearly independent of the elastic 
characteristics of the layers. The magnitudes of displacement necessary to initiate frac- 
ture at some point in the layer are small (3.7-8.2 m) for the three types of applied dis- 
placement. For applied displacement (1) and (3), the initial fracture is a vertical tensile 
crack at the crest of the fold. For applied displacement (2), the initial fracture is either 
a vertical crack at the upper surface or a shear fracture at the lower surface. 
Displacement fields and fracture patterns for scale-model experiments of two problems 
similar to the analytical examples are presented. For a sinusoidal vertica! displacement, 
the fracture pattern is a complex zone of normal faults which taper inward toward the 
axis of the fold and die out at depth. For a step in vertical displacement, the fracture 
pattern is (1) a series of reverse faults which start vertically at the base of the layer, curve, 
and intersect the upper surface at low angles, and (2) a series of normal faults dipping 
toward the convex side of the reverse faults. Of particular interest are the reverse faults 
which show that vertical movement at depth can generate low-angle faulting at the sur- 


face. 


Displacement fields found in the elastic analyses are good first-order approximations 
of displacement fields in the scale-model experiments. Points of initial fracture observed in 
the model experiments agree closely with those computed in the elastic analyses. The 
line of fracture for the curved reverse faults in the model experiments can be predicted on 
the basis of the Mohr fracture criterion and the stress distribution from an elastic analysis. 
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INTRODUCTION 
Much of the structural deformation in 
sedimentary basins results from differential 


vertical movements of underlying basement 
rock. This paper describes an analytical and 
experimental investigation of deformation 
resulting from simple distributions of vertical 
displacement along the bases of homogeneous 
rock layers. A special form of the general theory 
of elasticity was used in the analytical work, 
and scale models with homogeneous layers of 
dry sands and clay were used in the experi- 
mental work. 

Purposes of the investigation were to de- 

termine: 

(1) The characteristics of the folding and 
faulting produced by simple two-di- 
mensional distributions of applied verti- 
cal displacement along the bases of 
homogeneous layers. 

(2) The degree of similarity between a two- 
dimensional elastic analysis and a two- 
dimensional scale-model experiment of 
the same physical situation. 

The physical situations investigated were 
simple idealized representations of the geologic 
conditions in sedimentary basins. However, 
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investigation of these simple situations is a 
necessary first step toward an understanding of 
the more complex problems occurring in 
nature. 

The writer is indebted to Prof. C. H. Dix, 
under whose guidance this investigation was 
carried out at the California Institute of 
Technology, and to Prof. C. R. Allen, also of 
the California Institute of Technology, who 
read the manuscript and made many useful 
suggestions. Finally the author expresses his 
appreciation for two years of financial aid re- 
ceived through a Standard Oil Company of 
California fellowship. 


ANALYTICAL STUDY 


Method 


In the analytical method employed in this 
study, simple distributions of vertical displace- 
ment are specified along the base of a homo- 
geneous elastic laver. The method is similar to 
one used by Hafner (1951) to describe the 
response of a homogeneous elastic layer to 
applied stresses. The two methods give similar 
results when applied to simple geologic situa- 
tions such as folding of a layer. However, for 
certain structures, one method may be superior. 
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Block faulting of sedimentary layers above a 
rigid basement can be represented best by 
specifying displacements. In simple examples of 
this type of structure, adjacent blocks of base- 
ment rock are moved uniformly up or down with 
respect to each other. The overlying layer is 
undeformed by the movement except in the 
region where the two blocks meet. The mag- 
nitude and distribution of stress along the base 
of the layer cannot be determined in the zone of 
deformation. The best analytical representa- 
tion can be made by using the only available 
information—the displacements. 

An analytical method based on displace- 
ments along the base of an elastic layer has an 
additional advantage. Displacement can be 
easily measured in scale-model experiments. 
Therefore a comparison between model ex- 
periments and elastic analyses is possible when 
displacements similar to those used in the 
experiments are specified in the elastic analyses. 


Derivation of Analytical Method 


The analytical method is based on the theory 
of elasticity (Timoshenko and Goodier, 1951; 
Love, 1944; Muskhelishvili, 1953). Several 
assumptions made in this theory restrict its 
application to idealized representations of 
geologic structures. The displacements and 
strains must be small. The material undergoing 
deformation must be perfectly elastic (able to 
recover its initial form after removal of stresses), 
linearly elastic (possessing a linear relationship 
between stress and strain), homogeneous 
(possessing the same specific properties through- 
outs its volume), and isotropic (possessing 
elastic properties which are the same in all 
directions). 

Even with these assumptions, application of 
the theory of elasticity to three-dimensional 
geologic problems remains difficult. However, a 
simplification is possible within the framework 
of the theory. If the lower boundary of a layer 
undergoes displacements in the x- and y- 
directions only, and these displacements are 
identical for all xy-cross sections along the 
z-axis from minus to plus infinity, then all 
xy-cross sections are in the same condition 
(see Fig. 1). Therefore, a two-dimensional 
analysis of a single cross section will describe 
the stresses and displacements along the length 
of the z-axis. With this procedure, two-dimen- 
sional situations similar to cross sections near 
the center of elongated geologic structures can 
be analyzed. 


The nomenclature used is as follows (notation 
is the same as that of Timoshenko and Goodier 
(1951, p. xvii)): 


Normal components of stress parallel 
to the x-, y-, and s-axes 
Tay Shearing-stress component 


Or ,Ty, Gz 


o1,02,0; Maximum, intermediate, and minimum 
principal stresses 
u,v Components of displacement in the 
x- and y-directions 
Ou ov 
€é¢ = —, € = — Unit elongations in the x- and 


- 7] 
Ox : oy 
y-directions (strains) 


ou ov 
xy = — — Shearing strain 
oy Ox 
G Modulus of rigidity 
v Poisson’s ratio 
r Lamé’s constant 
: \+ 2G 2(1 — v) ; 
Ki, = * = ———— Constant used in this 
G 1 — 2 
paper 
: : 1 age 
Ks = 1 — K, = — ——— Constant used in this 
1 — 2p 
paper 
p Density 
g Gravitational acceleration 
ae oe 
Vv? = — + — Laplacian operator 
Ox oy" 
a = mr/L Coefficient 
Ea Distortional strain-energy density 
i Thickness of the elastic layer 
gE Distance between points of zero vertical] 
displacement 
B. Critical displacement—the maximum 
amount of displacement without frac- 
ture 
o Normal stress across a potential frac- 
ture plane 
T Shearing stress along a potential frac- 
ture plane 
TO Shear stress necessary to overcome the 
initial shear strength due to cohesion 
o Angle of internal] friction 
0 Angle between line of fracture and 


principal stress 


The sign convention used for normal and 
shearing stresses is shown in Figure 1. Positive 
normal stresses are tensile, negative normal 
stresses are compressive. Positive shearing 
stresses have the directions shown in Figure 1, 
negative shearing stresses have directions 
opposite from those shown. The displacements 
wu and v are positive if in the direction of 
increasing x and y, respectively. 

The solution of any two-dimensional elastic 








22 A. R. SANFORD—STUDY OF SIMPLE GEOLOGIC STRUCTURES 





—— 
-_— 


— Re ee) aga Oe 





x 


Ley (+) 














L ; 
—X _ aes ¢) 
Fy Nees 




















/ 
4 
yi 
4 
Ps 
y, 2 
- ° fe) 
) Leo) Fo 
—— Ft 
- 
ee ae —— a H 
Le ee — a 
Jy 








y 


FIGURE 1.—C0-ORDINATES AND STRESS SIGN CONVENTION FOR THE ELAstTIC PROBLEMS 


For two-dimensional problems, displacements are identical for all xy-cross sections along the z-axis from 
minus to plus infinity. Positive normal stresses are tensile; negative normal stresses are compressive. Posi- 
tive shearing stresses have the directions shown; negative shearing stresses are opposite to those shown. 


problem must satisfy the following equations 


(Timoshenko and Goodier, 1951, p. 11-27): 
do OTe, P 
— “=40, (1) 
Ox oy 
do, OT x, 
ees c 4 of (2) 
oy Ox 
Va, + ay) = 0. (3) 


(These equations and those which follow omit 
the terms pertaining to the weight of the 
layer. The effect of the weight of the layer is 
incorporated into the solutions in the next 
section of the text.) Equations (1) and (2) are 
the differential equations of equilibrium which 
insure a state of static equilibrium. Equation 
(3) is the compatibility equation (expressed in 
terms of stress) which insures a proper relation- 
ship between the three components of strain, 
so that two-dimensional deformation can occur 
without discontinuities in displacement (Musk- 
helishvili, 1953, p. 95-97). In addition to 
satisfying equations (1) through (3), the 
solution must satisfy specified stress and/or 
displacement boundary conditions. If the 
solution of the problem gives stresses and dis- 
placements which fulfill all these conditions, 
then it can be shown that the solution is 
unique (Timoshenko and Goodier, 1951, p. 
236; Muskhelishvili, 1953, p. 66-71). 

The stresses in equations (1) through (3) 
are related to the strains in the following 
manner: 


Ou | Ov Ou 
At— + - + 26 —, (4) 
Ox oy Ox 
Ou Ov , Ov re 
= AES + SF + =, (5) 
Ox OV ov 
_ fou ov 
ty = G{—+—). (6) 
oy Ox 


Rewriting equations (1) through (3) in terms of 
strains gives: 


Ls Ou Ov a 
(r a G) — (% + \ aa GV*(1) = 0, (7) 


Ox \Ox oy 





— ou ov 
(A + G) — (* + - ) + GV%r) = 0, (8) 


dy \ax ay 


: Ou Ov 
v3 | 20 + G) ( - + *) = 0. (9) 
OX oy 


The displacements can be expressed as algebraic 
sums of partial derivatives of a scalar potential 
@(., y) and a vector potential V(x, y), in the 
following way (Phillips, 1933, p. 186): 


Ob av ; 
“t= - ag (10) 
Ox oy 
ab OV 
po + —, (11) 
ov Ox 


Substituting derivatives of equations (10) and 
(11) into equations (7) and (8) gives the 
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equations of equilibrium in terms of the two 
potentials: 


0, (12) 


o 0 
(. + 2G) — (Wb) + G — (vw) 
Ox oy 


a) 0 
(\ + 2G) — (Vb) — G— (Vw) = 0. (13) 
oy Ox 
Differentiating equation (12) with respect to x 
and equation (13) with respect to y and adding 
the two equations gives: 


VV) = 0. (14) 


Similarly, differentiating equation (12) with 
respect to vy and equation (13) with respect to 
rand subtracting one equation from the other 
gives: 


v7(Tw) = 0. (15) 


Equations (14) and (15) require that @ and W be 
either harmonic (Vb = Q) or biharmonic 
functions (Vb = 0). Substitution of partial 
derivatives of equations (10) and (11) into the 
equations for o, and o, (equations 4+ and 5) 
and summing the results gives: 


or to, = 2(. + G)V*. (16) 


The sum of o, and o, is zero when ®@ is a 
harmonic function. This condition cannot be 
satisfied by deformed layers (for example, at a 
free surface where o, = 0 everywhere and o, 
has finite values). Therefore ® must be a 
biharmonic function. 

The solutions for ® and WV must satisfy the 
biharmonic differential equations and_pre- 
scribed boundary conditions. Of interest here is 
the deformation of an elastic layer of thickness 
H which has the following stresses and dis- 
placements specified on its boundaries: 


Atv =0, o, = 0,and7,, = 0 (stresses zero); 


Aty = H, » = —B cos ax, and u = BP sin ax, 


where oO< P< 1. 


The lower surface of the layer undergoes vertical 
and horizontal displacements proportional to 
the cos ax and sin ax, respectively. Assume 
that the vertical and horizontal displacements 
within the layer are also proportional to the 
cos ax and sin ax. Equations (10) and (11), 
which express displacements in terms of the 
potential functions, indicate that the above 
distribution of displacements may be possible 
with potential functions of the form: 


(17) 
(18) 


© = cosax F(y), 


Vv 


sin ax f(y). 


At this point, there is no guarantee that 
equations (17) amd (18) will lead to a solution of 
the problem. Equations (17) and (18) are the 
correct expressions for the potential functions 
only if the stresses and displacements derived 
from them satisfy equations (1) through (3) 
and the prescribed boundary conditions. 

F(y) and /(y) can be determined by sub- 
stituting equations (17) and (18) into equations 
(14) and (15) and solving the resulting fourth- 
order differential equations. This leads to gen- 
eral expressions for ® and W. 


® = cos ax[Ayet” + Aoyer 
+ Ase + Agyee], (19) 
W = sin ax[Asee” + Agvet” 


+ Aze-ev + Asye%], (20) 


where A, through As are arbitrary constants. 

Equations (19) and (20) contain eight 
arbitrary constants, whereas only four bound- 
ary conditions are specified in the problem. 
Four extra arbitrary constants appear in the 
equations for the potential functions as a result 
of (1) expressing the displacements as partial 
derivatives of potential functions (equations 
10 and 11), and (2) differentiating the equations 
of equilibrium to obtain the two biharmonic 
equations. Some of the constants in the general 
expressions for ® and W are related. Substitut- 
ing equations (19) and (20) into equations (12) 
and (13) gives: 


1+ 2G (1 — ») 
Z = Ao = oe Ae = ¢y. 9 
{6 ( Ps ) seal edie 


r G 
Ag = ( ‘se = ') A, 
G 


_ _21 — ») 
we 2s 


(22) 
A, = —K,Ay,. 


These relationships eliminate two of the extra 
constants. 

The remaining two extra constants were 
eliminated by trial and error until equations for 
® and W were found which satisfied the four 
boundary conditions and the equations of 
equilibrium and compatibility. The equations 
for ® and WV satisfying the requirements for a 
solution of stress and strain in the elastic layer 
are: 


© = cos ax[Ai(e™ + e%™) + Aoye™ + Ayye*], 
(23) 


W = sin ax[As(ev + e-@) 


+ Ky Asver¥ om K,Agye~2"}. (24) 
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The constants in the above equations are equal 
to: 

A, = 0, 


ee Bal» ae Ba 
“NN 
{Kl — Phe” — (1 + PY + alle *| 
~ 2K, — Pye@*}, (25) 

Ba Ba 

A; = 7 = ; 

‘ A 

KSC + Pleo” + (1 — P)QaHl — 1)e*"| 


— AK, + Pre*"}, (26) 


1, BCs _ BKs 
——_—° - oa 
-[((A, + all PR;)(e0" + & @") 
+ (P — allK;)(et — e& @*)|, (27) 
where 
N = aN; 
= a[K3(2 — e&4 — ¢ 2H) — Ki K(2 + eat 
+ eH) + 4e7//72?K32 + 4 K,), (28) 
and 
: ; 1 : 
K; = (1 — Ki) = — ——__, (29) 
1 — 


The general equations for v7, #, 7, , o,, and 
Try are: 


B cos ax 


ay 


9 


[C2e@™ayKs + 1) 


— Cye™ayK; — 1) — Clee” + e@)], (30) 


B sin ax : : : 
4 = 3 [—Cre(ayK; = K)) 

™M 
— Cye a(ayK; + K)) a Clee = ¢ ay), (31) 

2GBa cos ax 
1: ell > F 
N; 
[Ce Kzaye! + Cy Kz; ave av — C,(e — e ay) | (32) 
2GBa cos ax 
oO; = —d, = 
™M 


-[—2K3(Cse — Cye-@)], (33) 


2GBae sin ax , : 
Try = y [C;(e% + e av) 
i 
_ Cs K;e*(ay -+ 1) + C;K3e Kay — 1)]. (34) 


Also of interest in elastic analysis are the 
equations for the principal stresses. The prin- 
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cipal stresses at a point within a stressed layer 
are the stress values normal to three per. 
pendicular planes across which the shearing 
stresses vanish. In a two-dimensional analysis, 
the intermediate principal stress (a2) parallels 
the z-axis. The equations for the principal 
stresses in terms of o,, o,, and 7, are: 


a; o By V ‘ww be 2, (35) 


o,+o, or — o,\? / 
ie ns = V ( = ‘) oe 
? / ? a 


The response of an elastic layer to general 
distributions of displacement along its lower 
boundary can be determined by superposing 
equations of the type which have been shown 
herein (Timoshenko and Goodier, 1951, p. 235), 
Combination of equations is accomplished by 
means of a Fourier series, which can be used to 
approximate a variety of geologically reason- 
able displacements at the base of a sedimentary 
layer. 

Superposition can also be used to incorporate 
into the equations an initial state of stress due 
to the weight of the layer. In this work, a state 
of hydrostatic stress: 

or = —pgy, o, = —pgy, and r,, = 0, 
which satisfies the equilibrium and com- 
patibility requirements, has been superposed 
with stresses arising from the deformation of 
the layer. This superposition of equations gives 
the correct stress values for layers which are 
assumed to be in hydrostatic equilibrium prior 
to deformation. 


Application of Analytical Method to Examples 


The internal stress distribution of a deformed 
elastic layer can be portrayed by a set of orthog- 
onal curves called stress trajectories. These 
stress trajectories are tangent to the directions 
of the principal stresses at all points within 
the layer. The equation for the direction of one 
of the principal stresses (0; or oy) in the two- 
dimensional case is: 

2rzy 


tan 28 = ? (37) 
Or — Gy 


where @ is the angle between the positive .r-axis 
and the principal stress measured in the direc- 
tion of the positive y-axis. In the numerical 
examples discussed in following sections, a 
graphical construction of the stress trajectories 
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was made on the basis of 8 values computed for 
a grid of points within the layer. 

The only factors which influence the orienta- 
tions of the stress trajectories are the shape of 
the layer (ratio of H to L) and Poisson’s ratio. 
Rigidity (G), maximum applied vertical dis- 
placement (B), and hydrostatic stress cancel 
out of equation (37. The shape of the layer 
(H/L) and Poissor:’s ratio appear in the stress 
equations through the constants Ai, K3, 
C. Ce, and Cs. 

The distribution of displacement in a de- 
formed layer can be portrayed by a field of 
displacement vectors. The orientation and 
length of each vector indicate the direction and 
magnitude of the net displacement at the point 
from which the vector originates. 

The only factors which influence the orienta- 
tions and relative magnitudes of the displace- 
ment vectors are the shape of the layer and 
Poisson’s ratio. Rigidity (G) does not appear 
in the displacement equations. The shape of the 
layer (7/L) and Poisson’s ratio enter into the 
displacement expressions through constants 
Ki, Ks ’ Ce ’ C4 ’ and C5 ° 

Stress distributions and displacement fields 
of different examples can be compared if each 
numerical example sustains an amount of 
applied displacement which is just necessary to 
initiate fracture at some point within the layer. 
In the numerical examples which follow, dis- 
placements necessary to initiate fracture were 
computed on the basis of a special case of the 
Mohr fracture criterion. The equation for this 
fracture criterion is: 


T=) —o tang, (38) 


where 7 and o are the shearing and normal 
stresses across a potential fracture plane, 79 is 
the shear stress necessary to overcome the 
initial shear strength due to cohesion, and 
¢ is the angle of internal friction. The shear 
strength is dependent on the type of stress 
(tensile or compressive) and the amount of con- 
fining (hydrostatic) stress. The criterion is 
diagrammed in Figure 2. 

This special form of the Mohr criterion was 
selected because it satisfactorily describes the 
shear strength of many sedimentary rocks under 
the confining pressures encountered at depth in 
sedimentary basins (1000-2000 kg/cm?) 
(Jaeger, 1956, p. 82). For zero or low confining 
pressures, the criterion may not explain tensile 
fracture of rock subjected to tensile stress. 
Inasmuch as present experimental data are 
inadequate for establishing a fracture criterion 


bo 
mn 


for rocks under tension, the assumption was 
made that the Mohr criterion does describe 
fracture under tensile stresses. 
Critical displacements in the 
examples were determined as follows: 
(1) Principal stresses were found for a grid of 
points by substituting computed values of 
or, 0,, and Tz, into equations (35) and 
(36). However, in computing these stresses, 
B (which appears in the equations for 
O,,0,, and 7,,) was not given a value. By 
carefully examining these computed prin- 
cipal stresses, all but a few points in the 
layer could be eliminated as probable 
points of fracture. 
A displacement (B) necessary to initiate 
fracture at each of the probable points of 
fracture was found by substituting the 
values of the principal stresses (expressed 
as fractions of B plus a hydrostatic stress 
term) into a rearranged form of equation 
(38): 


fe sll ‘ 
oa = 7 COS ? =e (4 ) sin ?, (39) 


numerical 


) 


<> 
4 


and solving for B. 

(3) The point giving the smallest value of B 
was selected as the point of fracture, and 
the computed B was taken as the critical 
displacement B, . 

The Mohr criterion explains satisfactorily the 
shear strength of brittle materials which are 
generally elastic up to the point of fracture. 
However, some materials, after an initial period 
of elastic deformation, yield by plastic flow 
rather than by fracture. In the elastic region 
of the deformation, the analysis of stress and 
displacement is the same regardless of the type 
of failure which occurs. The initiation of plastic 
flow, however, is probably controlled by a dif- 
ferent condition than fracture in brittle sub- 
stances. Flow is believed to occur when the 
strain energy of distortion per unit volume 
reaches a value characteristic of the material 
(Jaeger, 1956, p. 93). The formula for this 
criterion in the two-dimensional plane strain 
case is (Timoshenko and Goodier, 1951, p. 
149): 


1 
Eq = iG {for — o,]? + lo, — v(o, + oy)? 
+ [v(oz + oy) — o,/? + 6r,,7}. (40) 
This criterion, credited to R. von Mises, differs 


from the Mohr criterion in that initiation of 
plastic flow is independent of confining stress 
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and occurs regardless of whether the stresses are 
tensile or compressive. 

In the numerical examples, values of dis- 
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placement specified along the lower boundary 
of the elastic layer (see Fig. 3). In Category I, 
the lower boundary of an elastic layer undergoes 
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Figure 2.—DIAGRAM OF THI 


= Mower FRACTURE CRITERION 


7 and a are the shearing and normal stresses across a potential fracture plane, 7» is the shear strength due 


to cohesion, ¢ is the angle of internal friction, and @ 
of fracture. 


using the values of stress dictated by the Mohr 
fracture criterion. The values obtained in this 
manner are not related to any characteristic 
values for rock materials. However, the dis- 
tributions of distortional strain-energy density 
do show the regions susceptible to flow in 
materials which favor this mode of failure. 


Numerical Examples 


the 
are 


Introduction —Table 1° summarizes 
numerical calculations. The examples 
categorized according to the nature of the dis- 


is the angle between the stress trajectory and the line 


horizontal displacement. Equations (25) 
through (34) (with P = 0) were used for the 
numerical examples in the first category. In 
Category II, the lower boundary of an elastic 
layer undergoes an approximate step in vertical 
displacement and no horizontal displacement. 
Eight solutions of the type given by equations 
(25) through (34) (with P = 0) were super- 
posed by means of a Fourier series to obtain the 
numerical example in this category. In Cate- 
gory III, the lower boundary of an elastic layer 
undergoes sinusoidal vertical and_ horizontal 
(90° out of phase with the vertical) displace- 
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ments. Equations for the numerical example in 
this category are given in a later section. 

The boundary conditions at the upper sur- 
face of the elastic layer are 0, = 0 and Tz, = 0. 


0 L/2 


The cohesive strength, 250 kg/cm’, was cal- 
culated from crushing-strength data on sedi- 
mentary rocks (Birch ef al., 1942, p. 116; 
Marin, 1952, p. 381; Mills et al., 1955, p. 310) 
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CATEGORY IL 
FicurE 3.—CLASSIFICATION OF NUMERICAL EXAMPLES 
Applied vertical displacement in Category I is sinusoidal; applied vertical displacement in Category 


II is an approximate step. 


These conditions specify a free upper surface 
which corresponds to the surface of the earth. 
In Categories I and II, the boundary conditions 
on the lower surface are (1) no horizontal dis- 
placement and (2) a specified distribution of 
vertical displacement. The latter conditions 
approximate the situation along the lower sur- 
face of an elastic sedimentary layer which is 
deformed by vertical movement of rigid bed- 
rock beneath it. The contact between rigid 
bedrock and elastic sedimentary rock is 
“welded” so that horizontal movement along 
the contact is prevented. 

Average values of the properties of sedi- 
mentary rocks are used in each example. The 
density, 2.5 gms/cm*, and Poisson’s ratio, 
0.25, are based on the values listed in Birch 
et al. (1942, p. 8-37, 76). The rigidity value 
adopted, 2.0 X 10" dynes/cm?, is Gutenberg’s 
estimate of rigidity for Tertiary sandstone at a 
depth of 2 km (Gutenberg, 1951, p. 367). The 


. angle of internal friction, 35°, is based on 


experimental evidence cited by Hubbert (1951). 


and the adopted angle of internal friction of 
3S. 

Results of the numerical calculations are 
diagrammed in Figures 4-6 and 10. These 
figures are drawn for a short portion of a layer 
between the limits , = 0 and x = L. However, 
the symmetry of the problem permits visualiza- 
tion of results for longer layers if the drawings 
are rotated 180° about vertical axes in the plane 
of the figure. 

The displacements shown in these figures are 
only those resulting from the application of 
vertical displacements on the lower boundary 
of the layer. The stress values, on the other 
hand, include the initial hydrostatic stresses as 
well as the stresses arising from the deformation 
of the layer. 

Category I examples.—Displacement fields for 
Category I examples are shown in Figures 44 
and 5A. Theoretically, two factors influence the 
relative magnitudes and orientations of the 
displacement vectors in these figures—the 
dimensions of the layer and Poisson’s ratio. 
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TABLE 1.—SUMMARY OF THE NUMERICAL EXAMPLES Fiat 
= — . - ‘ ——— y 
Category I Category II Category III rg 
saninaaite laniianeiadied ie ae 
At y = 0 oy = 0; 72, = 0 o, = 0; 72, =0 oy=0;7ra9 | 
itis rTy phe. wit 
aed | v= = > oe Tx “i 
Conditions | ie it oe <eliton = Sa w  (2n—1) |oy = —A aa; oi 
' -cos(2n — 1) =: =0 Try =0 | rati 
_ L 0 
Examples Ia Ib Ila Illa 1B 
a ah , stre: 
| fn 
1. Dimensions of layer mat 
H—km 5 5 5 5 und 
L—km 25 15.7 25 15.7 3B. 
2. Elastic properties poin 
v—Poisson’s ratio 0.25 0.25 0.25 0.25 type 
G—10™ dynes/cm? 2.0 2.0 2.0 2.0 is a 
3. Fracture properties max 
o—degrees 35 35 35 35 disp 
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y 0 0 0 0 ing | 
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maximum stresses 4Ca 
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« 0 0 0.451 L/2 the 
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o,—kg/cm? — 1282 — 1324 — 2229 — 1292 fold 
x 0 0 0.451 L/2 fract 
y iH H H H Frac 
Tzy—kg/cm2 74 103 489 63 high 
x L/2 L/2 L/2 0 comy 
y i I H H/2 suffic 
6. Magnitude and location of Moh 
maximum displacements at Ca 
the upper boundary displ 
v—meters 738 3.9 Se 5.7 amp! 
x 0 0 0 L/2 maxi 
u—meters 3.8 2.4 3.0 y given 
x L/2 L/2 L/2 0 Fis 
7. Magnitude and location of Cate; 
maximum displacements at First 
the lower boundary hotte 
v—meters 8.2 4.5 3.7 5.8 ends 
x 0 0 L/2 block 
u—meters 0 0 0 2.4 regio 
x - = 0 are n 
disres 
Figure 7 shows the relative changes in the The stress distributions for Category I ex- detor 
maximum displacements at the boundaries with amples are shown in Figures 4B and 5B. The ra 
changes in these factors. The shape of the orientations of the stress trajectories in these _" 
layer hasa strong influence on the displacement figures are determined by the shape of the layer | “Y*' 
field, whereas Poisson’s ratio has little effect. (H/L) and Poisson’s ratio. Comparison of — 
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Figures +B and 5B shows the influence of shape 
on the direction of the stress trajectories. 
Stress distributions calculated for two other 
examples (layer dimensions same as Fig. 4B) 
with Poisson’s ratios of 0.167 and 0.333 were 
nearly identical in appearance. Between these 
two values, orientation of stress trajectories 
isnearly independent of variations in Poisson’s 
ratio. 

Of interest in the stress diagrams (Figs. 
4B and 5B) is the area of horizontal tensile 
stress centered over the crest of the deformed 
laver and extending to a depth of approxi- 
mately 1 km. Inasmuch as rocks are very weak 
under tensile stress, point F on Figures +B and 
5B, where the tensile stresses are greatest, is the 
point of initial fracture. For brittle rock, the 
type of fracture under the stress conditions at F 
is a vertical crack (Jaeger, 1956, p. 74). The 
maximum applied displacement or critical 
displacement before fracture occurs at F in 
Figures +B and 5B is 8.2 and 4.5 m, respec- 
tively. 

Two regions are equally susceptible to yield- 
ing by plastic flow according to the distortional 
strain-energy density distributions in Figures 
4C and 5C. One region is centered over the crest 
of the deformed layer, the other is centered in 
the trough of the deformed layer. The high 
distortional strain-energy at the crest of the 
fold reaches a maximum value at the point of 
fracture determined by the Mohr criterion. 
Fracture does not occur in the other region of 
high distortional strain-energy, inasmuch as the 
compressive stresses in this region are in- 
sufficient to initiate fracture according to the 
Mohr criterion. 

Category II example—Figure 6A shows the 
displacement field for the Category II ex- 
ample. The influence of layer shape on the 
maximum displacements at the boundaries is 
given in Figure 8. 

Figure 6B shows the stress distribution for the 
Category II example. Two features are notable. 
First, most of the vertical stress along the 
bottom of the layer is concentrated toward the 
ends of the adjoining uplifted and downdropped 
blocks. Second, the stresses in the central 
regions of the uplifted and downdropped blocks 
are nearly zero if initial hydrostatic stresses are 
disregarded. Both these features indicate that 
deformation of the layer is restricted to the 
region where the two blocks meet. In addition, 
the lack of stress in the central regions of the 
laver indicates that a uniform displacement 
superposed on any of the examples will not 
change the stress distribution. Therefore the 


results shown in Figures 4B, 5B, and 6B are 
also solutions for examples in which all vertical 
displacements applied along the lower bound- 
ary are in one direction. 

In Figure 6B, a region of nearly horizontal 
tensile stress lies just to the left of x = L/2, 
y = 0 and extends to a maximum depth of 1 
km. The point of initial fracture (F\) is at 
x = 0.375L, vy = 0, where these tensile stresses 
are greatest. The type of fracture under the 
stress conditions at F; is a vertical tensile crack. 
The amount of applied vertical displacement 
required to initiate fracture at F; is 3.7 m. 

The same applied displacement will cause 
fracture at Fy (x = 0.53L, y = H) as well as 
F, if the fracture properties (7) and @) of the 
layer are changed. Figure 9 shows the conditions 
under which fracture will occur simultaneously 
at F, and F»,. The stresses at 2; are compressive 
in both the x- and y-directions. The type of 
failure for brittle rock under these stress con- 
ditions is a shear fracture. 

Two separate regions of high distortional 
strain-energy density occur along the upper 
boundary of the laver (Fig. 6C). The maximum 
value in each of these regions is equal to the 
maximum values found in the Category I ex- 
amples. The highest distortional strain-energy 
concentration occurs at the bottom of the layer 
directly above the point at which the rate of 
change of applied displacement is greatest. 

Scope of numerical examples in Categories I 
and II.—Theoretically, two factors influence 
the orientations of the displacement vectors and 
the stress trajectories—the dimensions of the 
layer and Poisson’s ratio. The calculations 
which have been described indicate that layer 
shape has a strong influence on these orienta- 
tions, whereas Poisson’s ratio has little effect. 
As a result, the numerical solutions given here 
are applicable to any homogeneous elastic layer 
providing (1) its shape is the same as the shape 
of the layer in the example, and (2) an ap- 
propriate scaling factor is used to determine the 
magnitudes of stress in the stress diagram. The 
scaling factor is based on the rigidity of the 
layer and the maximum applied displacement. 
The maximum applied displacement must be 
less than the amount of displacement necessary 
to initiate fracture. 

Elastic properties of the layer do not have to 
be known for a comparison of elastic analyses 
and scale-model experiments. Displacements 
and fractures (see ‘‘Prediction of fractures’’), 
the only items easily observed in model ex- 
periments, are independent of the clastic 
constants in the elastic analyses. Therefore, 
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-EXAMPLE Ta DiaGrams: A—DISPLACEMENT FieLp; B—Stress DistripuTIon; C—DistorTIONAL STRAIN-F 


FIGURE 4.- 
In example Ia, the lower boundary of an elastic layer 5 km thick and 25 km long undergoes half a wave length of sinusoidal vertical displacement. 


geometrical similarity of layer shapes is the 
only requirement for a comparison between a 
two-dimensional elastic analysis and a two- 
dimensional scale model. 

Category III example.—A “welded” contact, 
which is specified along the lower boundary in 
Categories I and II examples, is one extreme of 
many possible contact conditions. The opposite 
extreme is a “frictionless” contact. The fric- 
tionless lower-boundary case was analytically 
investigated to determine the effect of contact 
conditions on deformation resulting from a 
sinusoidal distribution of vertical displace- 
ment. 

A frictionless condition at the lower boundary 
means that the shearing stress is zero along that 
boundary. The magnitude of applied horizontal 
displacement which gives zero shearing stress 
along the lower boundary is difficult to com- 
pute. The easiest approach is to specify stresses 
on the lower boundary as Hafner (1951) did in 
his work. 

The solution of problems in which all bound- 
ary conditions are expressed in terms of stress 
reduces to the solution of a single biharmonic 
equation which will lead to stresses satisfying 
equilibrium (equations 1 and 2), compatibility 
(equation 3), and the prescribed boundary 
conditions. Of interest here is the deformation 
of an elastic layer of thickness H under the 
following boundary conditions: 
Aty=0, o, =0, 71, =0 (stresses zero) ; 
Aty = H, o, = —Asinax, ry = 0. 


The biharmonic equation and the boundary 
conditions are satisfied by: 


® = sinax|[D, coshay + D2 sinh ay 
+ Dsy cosh ay + Dy sinh ay]. (41) 


The constants in equation (41) are equal to: 








D,; = 0, 
D; Dyk; 
dD: == 
a Qa 
Dl sinhall +all coshall] . 
u— se. (42) 
at all sinh aH 
D; = —D,k; 
_ _p,| Sithall tall coshall] (4, 
aH sinh all 
i a 
aie 
1 
(44) 


~ alle +42, sinhall — aH coshall)’ 
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In example Ib, the lower boundary of an elastic layer 5 km thick and 15.7 km long undergoes half a 
wave length of sinusoidal verticai displacement. 


The general equations for 9, , ¢2, Tsy, v, and 


“are: 


A sin ax 


‘(ay + ky) sinh ay — ayk; cosh ay), 


A sin ax 
o, = — 


ke 


‘(ay — ky) sinh ay + (2 — ayki) cosh ayl, (46) 


A cos ax 


a 


Try 


(45) 


-[((1 — ayh)) sinh ay + ay cosh ay), 


(47) 
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aa [(1 — 2y) sinh ay 





= 


+ ayk; sinh ay — 2h,(1 — v) cosh ay 


— aycosh ay] + constant, (48) 
A COS ax 
ee thal ~ S5 ibes 
u QaGk> [ i( v) sinh ay 
+ ay sinh ay + 2(1 — v) cosh ay 
— ayk, cosh ay] + constant. (49) 


The elastic constants appear in the displace- 
ment equations, but not in the stress equations. 

Table 1 summarizes the results obtained from 
numerical example IIIa. In example IIIa, a 
sinusoidal vertical stress is applied to the lower 
boundary of an elastic layer. The elastic con- 
stants and fracture properties for the layer are 
the same as in the previous examples. The 
results of the calculations are diagrammed in 
Figure 10. In addition, calculations were made 
to determine the influence of layer shape on the 
relative magnitudes of maximum displacement 
at the boundaries (Fig. 11). Fracture point and 
critical displacement were determined by the 
Mohr fracture criterion (equations 38 and 39). 

A comparison of examples IIIa and Ib shows 
the maximum influence of the conditions at the 
lower boundary on the deformation resulting 
from a sinusoidal distribution of applied 
vertical displacement. The layers in these 
examples have identical dimensions, elastic 
constants, and fracture characteristics. Both 
undergo the same distributions of vertical 
displacement along the lower boundary. The 
only difference is the contact between the layer 
and underlying material. Example IIIa has a 
frictionless contact, and example Ib has a 
welded contact. 

The following are not greatly affected by the 
nature of the contact: 

(1) The point of initial fracture as de- 
termined by the Mohr criterion (com- 
pare position of F on Figs. 5B and 10B). 

(2) The distribution of vertical and _hori- 
zontal displacements at the upper bound- 
ary as a function of the dimensions of the 
layer (compare Figs. 7A and 11). 

The features which change when the lower 
boundary is frictionless are: 

(1) The amount of applied vertical displace- 
ment before fracture is slightly greater. 
Two additional areas are susceptible to 
yielding by plastic flow (compare Figs. 
5C and 10C). 

Horizontal displacements nearly equal to 
those on the upper boundary occur along 


(2 


a? 


(32 


— 


the lower boundary (compare Figs. 5A 
and 104A). 

If displacement does occur along the lower 
contact of a layer, it will always be less than the 
amount found in example IIIa. An absolutely 
frictionless contact between rock layers is not a 
reasonable boundary condition in geologic 
problems. For this reason, most features of the 
deformation in a homogeneous layer will tend to 
be closer to those occurring in the “welded” 
contact case. 


Prediction of Fractures 


According to the Mohr criterion, fracture 
occurs along lines which intersect the principal 
stresses at constant angles. The angle between 
the line of fracture and the stress trajectories 
is equal to (see Fig. 2): 


cH) 
6= 45°+-. 
~ 


(50) 
The extent to which the stress distributions in 
the analytical work can be used in predicting 
fractures with the Mohr criterion is difficult to 
determine. Any fracturing within the layer 
changes the stress distribution determined in 
the elastic analysis. If this change is confined 
to regions near the fracture, then formation and 
propagation of fractures are controlled by the 
original stress distribution. On the other hand, 
if the change extends a large distance from the 
fractures, formation and propagation of frac- 
tures are controlled by a stress distribution 
which changes continuously as the fractures are 
formed. One method of determining whether 
or not stress distributions from an elastic 
analysis can be used to predict fractures is to 
perform experiments. 


EXPERIMENTAL STUDY 
Problems Investigated 


Two problems similar to the analytical ex- 
amples were investigated in the scale-model 
experiments. The first dealt with the deforma- 
tion of a homogeneous layer resulting from the 
application of a broad curve in vertical displace- 
ment along its lower boundary. The second 
problem dealt with the deformation of a homo- 
geneous layer resulting from the application of 
a step in vertical displacement along its lower 
boundary. 


Theory 


Only the scale-model theory needed for this 
discussion is summarized here. A more detailed 
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account of scale-model theory and its applica- 
tion to geologic problems is given by Hubbert 
(1937). 

A geologic structure and a scale-model repre- 





/ 
/ 























Pa é 4 
3 sentation of the same structure must have 
n . . . on 6 ea . . 
—s Z certain similarities. The model must be similar 
se in shape to the geologic structure; that is, 
88 Z lengths, areas, and volumes must be _pro- 
5 portional: 
| 49 
! evi Dg XDn, 
2 28 
a — a = 
z ac Am = L,?, (51) 
2 a 
z aS 
& od sa Vn = AFL,,’. 
2 ais 
5 Ze — 
8 ££ (Subscripts m and x indicate model and nature, 
¥ * respectively.) The mass distributions must be 
% ‘= similar; that is, the mass for any element ot 
2 2 volume in the model must be proportiona! to the 
n e . . 
% 2 mass of the corresponding element in the 
zy E geologic structure: 
Sk 2 
eS dMm = pdMy . (52) 
a > 
5 ~ The time required for any given change of shape 
4 = or position in the model must be proportional 
x 3 to the time required in the geologic structure: 
2 Bo ‘ 
= — 
zg Tn = 1Tn. (53) 
me 
a i ; : ; ; 
mb The constants of proportionality in these 
<- & equations, A, w, and 7, are the model ratios of 
=e length, mass, and time, respectively. From 
= % these model ratios, it can be seen that dynamic 
& % — similarity between the model and the original 
= 2 exists if the ratio of each force acting on an 
S'S — element of mass in the model is proportional to 
a. - . . 
= £ the same force acting on the corresponding 
ZZ element of mass in the geologic structure: 
~_ 3S 
| 8 
= dM », Ln iy 3 
net - —_ = -2 = 5 
ee = par? = @. (54) 
we dM, Ln Tn? 
aso 
22 
& 2 In the formation of geologic structures and 
=+ F ; A 
5 = also in the modeling of these structures, ac- 
mm OS > e 
ta «= celerations are generally so small that forces due 
<3 = to inertia are negligible. If inertia forces are 
1g 8 € negligible, dynamic similarity is completely 
-+22 46 — satisfied (for independent values of the model 
Ew a . ’ ° ¢ 
So .& ratios A, wu, and 7) as long as all forces conform 
os to the model ratio of the gravity forces (Hub- 


Ficure 6. 


bert, 1937, p. 1489). The geologic structure and 
the model are subjected to the same gravita- 
tional field, therefore, the model ratio of gravity 
force is: 
F, @dMag 
FP, dMng 
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Figure 7.—Carecory I: Errecr or Layer DimENSIONS AND Poisson’s RATIO ON 


DisPLACEMENTS AT LAYER 


BOUNDARIES 


A shows relative change in maximum displacements at layer boundaries with change in layer dimen- 
sions; B shows relative change due to change in value of Poisson’s ratio. 


From equation (55) the ratio of strength be- 
tween the model and the original structure is 
derived: 


o = pr, (56) 
The model ratio of mass can be expressed as: 


a= 5X3, . (57) 


where 6 equals model ratio of density. Substitut- 
ing equation (57) into equation (56) gives for 
the model ratio of strength: 
ao = Ox. (58) 
Replacing A with the ratio of lengths gives: 
5L - 
L,=—. (59) 
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of the material of the layer. 


Equation (59) can be used to establish the size 
of the geologic structure being modeled if the 
model ratios of strength and density are known. 


Description of Experiments 


Properties of modeling materials.—The com- 
position of the materials used in the experi- 
ments was: 
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COHESION—- Z 
FicuRE 9.—ExampLe IIa: Fracture Point LOCATION AS A FUNCTION OF LAYER FRACTURE 
CHARACTERISTICS 
Location of fracture point in example Ila (F; or F: on Fig. 6B) is a function of fracture characteristics 


Material 1 Beach sand. 

Material 2 Beach sand *= ..2r cent, clay 15 per cent 
(OM-4 Ker..ucky-Tenn. Clay Co.). 

Material 3 Coarse St. Peter sand (Ottawa Silica 

Sands, flint shot). 

Fine St. Peter sand 

Sands, #102). 


Density values for the two degrees of com- 
paction used in the experiments were: 


Material 4 (Ottawa Silica 
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Ifa Diacrams: A—DIspLaceMENT Fretp; B—Stress DIstRiBuTION: 
C—DistTorTIONAL STRAIN-ENERGY DENSITY 

In example IIIa, the lower boundary of an elastic layer 5 km thick 
length of sinusoidal vertical displaceme nt and half a wave length o 
(90° out of phase with the vertical). 
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6 
HAL 


Errect OF LAYER DIMENSIONS ON DISPLACEMENTS AT LAYER 


BOUNDARIES 
Figure shows relative change in maximum displacements at layer boundaries with change in layer di- 


mensions. + 

Uncom pacted Com pacted 
Material 1 1.3 gms/cm* 1.6 gms/cm’ 
Material 2 1.3 14 
Material 3 KS Li 
Material 4 1.4 7 


The first value is an average density for the 
material when it is poured into a container. 
The second value is an average density for the 
material when moderately tamped (at approxi- 
mately 10,000 dynes/cm*). 

No water was added to any of the materials 
used. However, a small percentage of moisture 
was absorbed from the air. Water-content 
measurements showed that (1) 2.4 per cent of 
the total weight of the clay in Material 2 was 
water, and (2) 0.9 per cent of the total weight of 
Material 2 was water. 

The distribution of grain sizes in the ma- 
terials is shown in Figure 12. 

Numerical values were assigned for the 
degree of roundness and sphericity of the sand 
grains by comparing the grain shapes with 
published charts (Krumbein and Sloss, 1951, 
p. 81). Grains which approach a perfect sphere 
have the highest values of roundness and 
sphericity (maximum value equals 1.0). 


Roundness Sphericily 
Range Average Range Average 
Materials 1 0.1-0.9 0.45 0.3-0.9 0.65 
and 2 
Material 3 0.8-1.0 0.90 0.8-1.0 0.90 
Material 4 0.1-0.9 0.65 0.3-0.9 0.70 


Measurements of the strength of the modeling 
materials were made in two ways. In the first, 
the fracture characteristics for the material, 
T) and @, were determined by placing the 
material in a rectangular container and remov- 
ing the support for one vertical edge. The angle 
of fracture caused by removal of support at one 
vertical edge is related to the angle of internal 
friction by: 


(60) 


If the material possesses cohesive strength, an 
unsupported vertical edge of a thin layer will 
stand without fracturing. The maximum thick- 
ness or critical height (4.-) without fracture is 
related to the cohesive strength (Tschebo- 
tarioff, 1952, p. 169-172): 


pENter 


= . (61 ) 
4 tan @ 


T0 


(Equations (60) and (61) assume the special 
linear case of the Mohr fracture criterion.) 
Equation (61) does not consider the horizontal 
tensile stresses introduced by movement of 
material prior to fracture. These horizontal 
tensile stresses are located near the surface close 
to the unsupported edge. Inasmuch as_ the 
tensile strength of granular materials is low, 
vertical tensile fractures form at the surface 
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before shear fracture occurs lower in the layer. 
Formation of the tensile cracks reduces the 
critical height. The maximum depth of the 


100 


material fractures according to the special 
linear case of the Mohr criterion. Although the 
data are scattered at low normal stresses, the 
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FiGURE 12.—DISTRIBUTION OF GRAIN SIZES IN MODELING MATERIALS 


cracks is estimated to be one-half the critical 
height (Tschebotarioff, 1952, p. 171). Thus, the 
maximum effect of the tensile cracks can be 
considered qualitatively by multiplying the 
observed critical height by 34. 

This method was applied to all modeling 
materials with the following results: 


Uncom pacted 
ro(dynes/ 


Com pacted 
to(dynes/cm?) @ 


cm?) 
Material 1 0 34° 0 54° 
Material 2 140-210 42° 1100-1700 58° 
Material 3 0 22° 0 45° 
Material 4 0 27° 0 52° 


The strength measurements were taken under 
conditions similar to those used in the ex- 
periments. In the uncompacted tests, material 
was poured into the box. In the compacted 
tests, material was placed in the box by layers 
which were moderately tamped (approximately 
10,000 dynes/cm*) after each layer was added. 

The second procedure used for measurement 
of strength was the controlled-strain shear test 
(Tschebotarioff, 1952, p. 143-145). With these 
tests, the shear strength of Material 2 was 
determined under very low confining pressures. 
Tests were conducted for two degrees of com- 
paction, strong and moderate (Fig. 13). The 


best curves drawn by eye give the following 
values of vo and ¢ for Material 2: 


Strong Compaction Moderate Compaction 
TO 4000 dynes/cm? 1000 dynes/cm? 
¢ 53° 40° 


Also plotted on Figure 13 are data for con- 
trolled-strain shear tests on moderately com- 
pacted dry clay. 

Model Apparatus —The model experiments 
were performed in a rectangular box (inside 
dimensions 44.5 by 22.8 by 11.4 cm) with glass 
sides and wooden ends (Pl. 1). The bottom of 
the box was fitted with three felt-edged blocks 
for experiments that required steps in vertical 
displacement at the base of a layer. Each block 
could be moved independently in the vertical 
direction by four screws projecting through the 
bottom of the box. 

The bottom of the box was fitted with a felt- 
edged rubber pad for experiments that required 
gradual changes in vertical displacement at the 
base of a layer. The central and end portions of 
the rubber pad could be moved independently 
in the vertical direction by screws projecting 
through the bottom of the box. 

At the two ends of the box were vertical 
blocks which followed the vertical movement of 
the end blocks or end portions of the rubber 
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pad. This feature eliminated the drag which 
would occur if the ends of the box were fixed. 

One feaiure of the apparatus may influence 
the experimental results—the drag of the 
material along the glass face. Theoretically this 


compaction were used after each layer was 
added—zero and approximately 10,000 dynes/ 
cm?, 

Displacements were applied to the base of the 
layer in small increments. Direct and oblique 
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Ficure 13.—REsutts OF CONTROLLED-STRAIN SHEAR TESTS 


drag should be zero if the model is to represent 
accurately the deformation normal to the long 
axis of a geologic structure. Controlled-strain 
shear tests were performed to find the coef- 
ficient of friction between Material 2 and a 
glass surface (Fig. 13). The normal stress of 
Material 2 against the glass face in the box 
ranged from 0 to 10,000 dynes/cm?. The 
maximum coefficient of friction for this range of 
normal stress is 0.12. This agrees closely with 
the coefficient of friction for quartz grains 
moving over a solid quartz plate (Tschebo- 
tarioff, 1952, p. 122-124). Inasmuch as the 
measured coefficient for Material 2 moving 
over a glass surface was only one-fifth to one 
twenty-fifth the coefficient of friction of 
Material 2 moving over Material 2, drag along 
the glass face probably did not have a major 
influence on the experimental results. 
Experimental procedure.—Material was 
placed in the box in successive layers 0.6-1.2 
cm thick. These layers were separated by thin 
marker lines of flour or quartz sand dyed black. 
The marker lines extended a short distance into 
the box from the glass face. Two degrees of 


photographs (Pl. 1) of the model were taken 
between intervals of displacement. In some 
instances, several increments of displacement 
were recorded in a multiple exposure (PI. 2). 
This was used to record displacement fields. 
The total time required for a complete ex- 
periment was about 1 hour. 


Results of Experiments 


General statement.—Two types of experi- 
ments were performed in the model studies. 
In Type I, a broad curve in vertical displace- 
ment was applied to the lower boundary of a 
homogeneous layer. In Type II, a step in ver- 
tical displacement was applied to the lower 
boundary of a homogeneous layer. Layers of 
each kind of modeling material were used one or 
more times in both types of experiments. The 
experiments performed were: 





Experiments With Material Total 
1 2 3. 4 Experiments 
Type I 1 2 3 1 7 
Type Il 2 13 3 1 19 
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Several figures showing different stages of 
deformation in the model experiments have 
been used in the discussion of experimental 
results. These figures are tracings of photo- 
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Figure 14.—DIspLACEMENT FIELDS 
Tyre I ExPeRIMENTS 


graphs taken during the experiments. The co- 
ordinates for the figures are the same as those 
used in the numerical examples (see Fig. 1). 

Type I displacement fields —Figure 14 gives 
examples of the observed displacement fields. 
The important features of the displacement 
fields were: 

(1) No measurable change in the displace- 
ment field with a change in the compo- 
sition of the layer. 

(2) An increase in horizontal displacement 
from the lower to the upper boundary of 
the layer. 

(3) Maximum horizontal displacement at or 
very near the point of inflection in the 
fold (x = L/2 or —L/2). 

(4) A decrease in vertical displacement from 
the lower to the upper boundary of the 
layer. 

Type I fractures with Material 2.—Only 
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models in which layers of Material 2 were used 
produced fracturing at large displacements 
(see No. 9 and No. 19 of Fig. 15). The sequence 
of fracture formation with increase in applied 
displacement was as follows: 

(1) Vertical tensile cracks at the crest of the 
fold. 

A series of normal faults at or near the 
crest of the fold. 

Additional normal faults intersecting the 
surface at progressevely greater distances 
from the crest of the fold. 

The depth of the tensile cracks was not greater 
than 1 cm and in most cases was much less. The 
major normal faults were straight, had mazi- 
mum displacement at the surface, and dipped 
toward the axis of the fold. Average dip of the 
major normal faults was 65°. 

The angle of fracture and the depth of the 
tensile cracks (equivalent to the unsupported 
height in the strength measurements) can be 
used to determine the strength of the material 
in the layer at the time of fracture. The com- 
puted angle of internal friction from equation 
(60) is 40°. The computed cohesive strength 
from equation (61) is 100—200 dynes/cm’, 
These values are nearly identical to the values 
found in the strength measurements of un- 
compacted samples of Material 2. 

Type I deformation with Materials 1, 3, and 
4.—No fracturing was observed in the ex- 
periments with layers of Materials 1, 3, and 4 
However, a cup-shaped area over the crest of 
the fold was disturbed by the folding (see No. 
23 and No. 21 in Fig. 15). The main features of 
the disturbed zone were blurring of the marker 
lines, flattening of the upper boundary, and 
thinning of the layer. The thinning over the 
crest of the fold was about 6 per cent of the 
original thickness for an applied vertical dis- 
placement of 1 cm at x = 0, y = H. 

Type II displacement fields —Figure 16 gives 
examples of the observed displacement fields. 
The important features of the displacement 
fields were: 

(1) No measurable change in the displace- 
ment field with a change in the compo- 
sition of the layer. 

(2) An increase in horizontal displacement 
from the lower to the upper boundary of 
the layer. 

(3) Maximum horizontal displacement at or 
very near + = L/2 or —L/2 along any 
horizontal line through the layer. 

(4) A progressively less abrupt transition in 
vertical displacement from the lower to 
the upper boundary of the layer. 
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FicurE 15.—FRACTURING AND FoLpiInc—Type I EXPERIMENTS 


Type II fractures with Material 2.—Type IU 
experiments with layers of Material 2 produced 
distinct fractures at large applied vertical dis- 
placements (Fig. 17). The usual sequence of 
fracture formation with increasing applied 
vertical displacement was: 

(1) Simultaneous or nearly simultaneous 
formation of tensile cracks at the upper 
boundary and shear fractures at the 
lower boundary. 

(2) Propagation of one or more shear frac- 


tures (reverse faults) to the surface and 
deepening and widening of the tensile 
cracks (Fig. 17). 

Formation of normal faults in the tensile 
crack zone (see No. 12 of Fig. 18). 

(4) Formation of a major normal fault from 
the tensile crack zone to the edge of the 
uplifted block at the lower boundary of 
layer (see No. 12 of Fig. 18). 

In most of the experiments a series of reverse 

faults rather than a single fault formed at the 
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edge of the uplifted block. Formation of these 
fractures followed a definite sequence; pro- 
gressively younger fractures started at pro- 
gressively smaller angles to the vertical edge of 
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faults (Figs. 17 and 18) was the curvature of the 
line of fracture away from the uplifted block. 
Two factors influenced the curvature of the 
reverse faults, (1) the thickness of the layer, 
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FiGurE 16.—DIsPLACEMENT FreELps—Type II ExpEerRm™Ments 


the uplifted block. The end result was a series of 
diverging fractures, the oldest starting at 
about 8° to the vertical, and the youngest 
starting vertically. Most of the fractures that 
started at an angle to the vertical propagated 
only a short distance. Many of these fractures 
were probably undetected in the experiments. 
Most of the fractures that started vertically 
eventually extended all the way to the upper 
surface. The applied vertical displacement that 
was required to propagate a shear fracture all 
the way to the upper surface was equal to about 
one-twentieth of the total length of the frac- 
ture. The same ratio between applied displace- 
ment and length of fracture (1 to 20) appeared 
to hold for all intermediate stages in the de- 
velopment of the fracture as well. 

The most interesting feature of the reverse 


and (2) the number of short reverse faults 
formed before a particular reverse fault propa- 
gated to the surface. The influence of the 
second factor on curvature was eliminated by 
omitting the vertical portion of the fractures in 
the curvature measurements. These measure- 
ments were then correlated with an effective 
thickness (H’) for each layer which was equal to 
the vertical distance over which the fracture 
actually curved (Fig. 17). 

Curvature measurements for three intervals 
of effective thickness are summarized in 
Figure 19. Although curvature measurements 
in adjacent intervals of effective thickness 
overlap, the mean values in each interval show 
a small systematic change of curvature with 
the effective thickness of the layer. With a 
decrease in the effective thickness, the frac- 
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tures have (1) progressively greater curvature 
and (2) progressively lower dip at the upper 
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were thrust out over the upper surface of the 
stationary blocks. In one experiment (Fig. 2 
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The curvature of the reverse faults was inde- 
pendent of the relative movement of the blocks. 
In one experiment, the center section was 
allowed to drop while the two end sections 
remained stationary. The reverse faults pro- 
duced by this mode of displacement were nearly 
the same as the reverse faults produced by 
keeping the center section stationary and 
pushing up the two end sections (Nos. 3, 4, 5, 
and 8 of Fig. 17). 

In all Type II experiments with Material 2, 
displacements occurred along the reverse faults 
alter they had reached the upper surface. As a 
result of this movement, thin wedges of material 


of the wedge was 0.5 cm after a total of 2.0 cm 
vertical displacement of the lower boundary. 

In most of the experiments, vertical tensile 
cracks appeared at the upper surface at nearly 
the same time as the reverse faults were formed. 
The measured positions of the tensile cracks 
for three intervals of effective thickness are 
summarized in Figure 19. There is a systematic 
shift of the crack position with a change in the 
effective thickness of the layer. 

Fracturing in the region of the tensile cracks 
followed the same pattern as the fracturing at 
the crest of the folds in the Type I experiments. 
The tensile cracks became deeper and wider 
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with increased applied displacement. Even- 
tually the crack reached a critical depth at 
which the material on cach side of the crack 
was unable to support itself, and normal frac- 
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However, more displacement occurred along 
the normal fault than along the reverse faults, 

The depth of the tensile cracks and the angle 
at which the reverse faults intersected the 
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FicuRE 18.—FRACTURING AND FOLDING 


tures dipping toward the tension zone were 
formed. The maximum depth of the tensile 
cracks before normal faults formed 
was 2.5 cm. 

The final fracturing in the experiments was 
the propagation of a steep normal fault (aver- 
age dip 78°) from the tensile crack zone to the 
edge of the uplifted block. After formation of 
the normal fault, displacement occurred along 
both the normal fault and the reverse faults. 


were 
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upper surface can be used to estimate the 
strength of the material in the layer at the time 
of fracture. The computed cohesive strength 
from the depths of the tensile cracks is 200- 
500 dynes/cm*. The angle of internal friction 
from the average angle at which the fractures 
intersected the free upper surface (29°) is 32°. 
These values are fairly close to those obtained 
in the strength measurements on uncompacted 
samples of Material 2. 
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(GRID LINES ARE 0.5 CM APART) 


Ficure 1. — Type I Experiment 


PRAT 
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Figure 2. — Type II Experiment 


OBLIQUE PHOTOGRAPHS 
OF EXPERIMENTS 
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No. 23 Material 1 
(GRID LINES ARE 0.5 CM APART) 


FicureE 1. — Type I ExperRtmENT — MULTIPLE ExPposuRE 








No. 15 Material 2 
(GRID LINES ARE 0.5 CM APART) 


Figure 2. — Type II Experiment — MULTIPLE Exposure 
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Type IT fractures with Materials 1 and 4. 
The fracture behavior of layers of Materials 1 
and 4 differed in two respects from the fracture 
behavior of layers of Material 2; (1) tensile 


L/2 
| 


placement fields, at small applied vertical dis- 
placements, appeared to be independent of the 
size, shape, and sorting of the grains used for 
the layers. Fracture, however, was influenced 























75 <H'< 10.0 CM 
5.O0<H'< 75M 
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FiGuRE 19.—CuRVATURE OF REVERSE FAULTS IN TypPE II EXPERIMENTS 
Reverse faults in Type II experiments had (1) progressively greater curvature and (2) progressively 
lower dip at the upper surface with a decrease in the effective thickness of the layer. 


cracks were not formed (No. 8 of Fig. 17); and 
(2) the reverse faults formed were not sharp 
breaks, but narrow shear zones (No. 25 of Fig. 
18). The curvature of the shear zones, however, 
was nearly identical with the curvature of the 
sharp reverse faults produced in the experiments 
with layers of Material 2. Layers of Material 1, 
which has a greater range in grain sizes and 
more irregularity in grain shape and roundness 
than Material 4+, had narrower zones of frac- 
ture than layers of Material 4. 

Type II deformation with Material 3.—Fold- 
ing rather than fracturing was dominant in 
experiments in which layers of Material 3 were 
used. However, there were broad zones in the 
layers where the movement of the sand parti- 
cles described trajectories which curved away 
from the uplifted block. In addition, ridges 
were formed on the upper surface similar to 
those associated with the reverse faults in the 
other experiments. 

Summary of experimental results —The dis- 


by the properties of the modeling materials. 
The sharpest fractures were formed in layers 
of material which had cohesive strength 
(Material 2). The next sharpest fractures were 
formed in layers of material which lacked co- 
hesive strength but were composed of poorly 
sorted, irregular sand grains (Material 1). 
Layers composed of well-sorted, well-rounded, 
and spherical sand grains (Material 3) did not 
appear to fracture. 

The strength of the material in the layer at 
the time of fracture was nearly the same as 
measured strengths of uncompacted samples 
of the modeling materials. Movement of the 
layer prior to fracturing apparently decreased 
the initial compaction of the material and thus 
reduced its strength. 


Geologic Interpretation of Experimental Results 


Size relationship between model and geologic 
structure.—The size relationship between a 
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model and the geologic structure that it repre- 
sents is determined from equation (59). The 
model ratios of strength and density must be 
known to use this equation. Unfortunately, 
there is no means of estimating the strength of 
extensive layers of sedimentary rock in nature. 
The strength is probably determined to a large 
extent by the degree of jointing in the layer. 
Layers of rock which have undergone extensive 
jointing have little or no cohesive strength and 
probably react to applied vertical displace- 
ments like layers of Material 1 do in the model 
experiments. Layers of rock which have not 
undergone jointing have a cohesive strength 
and probably react to applied vertical displace- 
ments like layers of Material 2 do in the 
experiments. 

A rock layer with no cohesive strength can 
be modeled with a dry sand (similar to Ma- 
terial 1) which also has no cohesive strength. If 
the angles of internal friction in the sand and 
the rock layer are equal, the strengths of the 
layers are identical. However, the cohesive 
strength of the rock layer is zero only if one 
considers an extensive layer. The individual 
pieces between joints have considerable co- 
hesive strength. The strengths of the sand and 
rock layers will be identical only when the 
dimensions of individual blocks are related to 
the size of the rock layer in somewhat the same 
manner as the individual grains are related to 
the size of the sand layer. 

Equation (59) cannot be used to establish a 
size relationship between the model and the 
geologic structure if the rock and the sand 
layers have no cohesive strength. However, a 
crude size relationship can be obtained by 
comparing the size of sand grains with the 
size of joint blocks. The mean grain diameter 
for Material 1 is 0.2 mm. A reasonable diameter 
for a joint block is 10 m. The size relationship 
based on the above values is 1 cm of Material 1 
equals 0.5 km of jointed rock. 

A relationship between the size of the model 
and the geologic structure can be established 
with equation (59) if the layer of rock in the 
geologic structure has a cohesive strength. A 
rock layer with cohesive strength can be 
modeled with Material 2 which has a small 
cohesive strength of 1000 dynes/cm? and an 
angle of internal friction of about 40°. Repre- 
sentative values of cohesive strength and angle 
of internal friction for sedimentary rock are 
250 X 10° dynes/cm? and 40°. The size rela- 
tionship, on the basis of the above values, is 


1 cm of Material 2 equals 1.5 km of sedimentary 
rock. ; 

The 79 and @ values used to establish this 
size relationship are based on controlled-strain 
shear tests of moderately compacted samples 
of Material 2 and triaxial tests on small samples 
of sedimentary rock. In the model experiments, 
the strength of a moderately compacted layer 
of Material 2 is considerably reduced by de- 
formation of the layer prior to fracturing. A 
similar reduction in the strength of a layer of 
sedimentary rock probably occurs in nature, 
However, the strengths of the two layers may 
not be reduced in the same proportion. For this 
reason, the size relationship between the model 
and the geologic structure may change during 
the course of an experiment. 

Fractures—The characteristic fracture pat- 
tern for Type I experiments is 1 complex zone 
of normal faults at the crest of the fold (Fig. 15), 
The fractured zone tapers inward to the axis 
of the fold and dies out at depth. The Kettleman 
Hills anticline (Woodring et al., 1940) has this 
type of fracture pattern. 

The characteristic fracture pattern for 
Type II experiments is a series of curved re- 
verse faults intersecting the upper surface at 
low angles, and a series of normal faults in the 
uplifted block (Fig. 18). The reverse faults 
start as vertical faults at the lower boundary 
but become thrusts at the upper boundary. The 
low-angle or thrust portion of the fracture at 
the upper surface is due to horizontal compres- 
sive stresses. However, the horizontal compres- 
sion is the result of vertical movement at depth 
and not of horizontal compression of the entire 
layer. 

The characteristic fracture pattern formed 
for each of the two types of experiments is 
perhaps the most significant feature of the 
model experiments. Each pattern is relatively 
complex in view of the simplicity of the two 
distributions of applied displacement. The 
individual fractures in each pattern are related. 
Formation of one fracture triggers another in 
a definite sequence. The type and location of 
the primary fractures, which initiate the frac- 
ture sequence, are controlled by the distribu- 
tion of applied displacement. The type and 
location of secondary fractures are determined 
partly by the distribution of applied displace- 
ment and partly by conditions arising from the 
formation of the primary fractures. 

The manner in which the fracture pattern 
forms in the two types of experiments suggests 
that other distributions of applied displace- 
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ment at the base of a layer may also lead to 
distinctive fracture patterns. If the models are 
a true representation of fracturing in nature, 
fracture patterns observed in the models might 
be compared with fracture patterns observed 
on the surface of the earth, and in this manner 
be used to determine the distribution of dis- 
placement at depth. 


COMPARISON OF EXPERIMENTAL AND 
ANALYTICAL RESULTS 


Type I Experiments and Category I Numerical 
Examples 


Displacement fields.—The displacement fields 
in Type I experiments and Category I numeri- 
cal examples are very similar, although this is 
not immediately apparent from an examina- 
tion of the displacement-field diagrams. Dis- 
placement fields in the analytical work (Figs. 
44 and 5A) are produced by applied vertical 
displacements which change periodically in 
amplitude and direction. Displacement fields 
in the experimental work (Fig. 14) are produced 
by applied vertical displacements which change 
in amplitude but not direction. The only major 
difference between the two displacement fields 
isa uniform displacement in one direction. As 
previously noted, uniform movement of a layer 
does not change the stress distribution. There- 
fore, superposition of a uniform upward dis- 
placement (equal to one-half the total ampli- 
tude of the applied vertical displacement) on 
the analytical displacement field will bring the 
two displacement fields into approximate 
agreement. 

Although the general appearance of the dis- 
placement fields is similar, they differ consider- 
ably in detail. The maximum vertical and hori- 
zontal displacements at the boundaries of the 
model and the corresponding theoretical dis- 
placements from Figure 7A (for the average 
H/L ratio of 0.85 in the experiments) are: 


Experimental Theoretical 


Mmax. Lower boundary 1.00 B 1.00 B 
%max. Upper boundary 0.60 B 0.40 B 
Mmax. Lower boundary 0.30 B 0 

Umax. Upper boundary 0.70 B 0.30 B 


(8 equals the maximum applied vertical 
displacement.) 


Two factors are responsible for most of the 
difference between theoretical and experimental 
values above. First, the boundary conditions 


in the experiments differed slightly from the 
boundary conditions in the analytical examples. 
In the experiments, the applied vertical dis- 
placement did not vary exactly as the cosine 
of the horizontal distance, and horizontal dis- 
placement occurred along the lower boundary. 
Second, the displacements required in the 
experiments to record the displacement fields 
photographically were much larger than the 
critical displacements (B,) in the analytical 
work. In nearly all the experiments, there was 
evidence of failure in the layer while the dis- 
placement field was being photographed. 

Fractures —In the experiments with layers 
of Material 2, the position of the vertical ten- 
sile cracks is the same as the position of the 
initial tensile fractures predicted in the numeri- 
cal examples. 


Type II Experiments and Category II Numerical 
Example 


Displacement fields —The general appearance 
of the displacement fields in Type II experi- 
ments (Fig. 16) and Category II example 
(Fig. 6A) is very similar if a uniform upward 
displacement (equal to one-half the total ampli- 
tude of the applied displacement) is superposed 
on the analytical displacement field. However, 
the displacement fields differ slightly in detail. 
Listed below are the maximum vertical and 
horizontal displacements at the boundaries of 
the model and the corresponding. theoretical 
displacements from Figure 8 (for the average 
H/L ratio of 0.47 in the experiments): 


Experimental Theoretical 


Vmax. Lower boundary 1.00 B 1.00 B 
Vmax. Upper boundary 0.95 B 0.89 B 
tmax. Lower boundary 0 0 

Umax. Upper boundary 0.60 B 0.67 B 


(B equals the maximum applied 
vertical displacement.) 


The difference between the experimental and 
theoretical values is small. In this case, the 
boundary conditions in the experimental work 
were nearly equivalent to the boundary condi- 
tions in the theoretical example. Fracturing in 
the layer while the displacement fields were 
being photographed probably accounts for 
most of the difference in the experimental and 
theoretical values. 

Fractures —The predicted positions of the 
reverse faults in the Category II numerical 
example are based on the stress distribution in 
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Figure 6B and equation (50). The constant 
angle obtained from equation (50) is the orien- 
tation of the reverse faults with respect to the 
position of maximum compressive principal 


Apparently the formation of the reverse faults 
shifts the position of the large horizontal ten. 
sile stresses closer to the position of the reverse 
faults. 
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EXPERIMENTAL FRACTURE 
COMPARISON OF TITEORETICAL 


Fictre 20. 


AND EXPERIMENTAL FRACTURES 


The shear fracture determined from Example Ila stress distribution (Fig. 6B) and Mohr’s fracture cri 
terion (with @ = 28°) agrees closely with the shear fracture in the experiments which were geometrically 


similar to Example Ia. 


stresses throughout the layer. Inasmuch as the 
principal stress trajectories are curved, the 
reverse faults will also be curved. In Figure 20, 
three reverse faults based on the stress distribu- 
tion have been drawn for three different angles 
of internal friction. Also shown on Figure 20 
is a shear fracture based on measurements of 
reverse faults in model experiments which 
were geometrically similar to the numerical 
examples. The agreement between the experi- 
mental fracture and the theoretical fracture 
for a @ equal to 28° is very good. A value of 
28° is a reasonable average angle of internal 
friction for the entire layer. 

The tensile cracks in the experiments do not 
occur at the positions predicted in the analytical 
work (compare the position of tensile cracks on 
Fig. 19 with the position of F; on Fig. 6B). 
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HYDROTHERMAL ALTERATION IN THE FRONT RANGE MINERAL 
BELT, COLORADO 


By F£éi1x GonzALEz BONORINO 


ABSTRACT 


Hydrothermal alteration associated with ore veins in a large part of the Front Range 
mineral belt was studied in an attempt to correlate it with the mineralogic, petrologic, and 
structural features of the ore deposits. Six patterns of alteration were found; they range 
from the simplest type with a single zone of hydromica to the most complex type with 
five zones that are characterized by phlogopite carbonate (fringe zone), montmorillonite, 
kaolinite, orthoclase, and hydromica (innermost zone). A regular mixed-layer illite- 
montmorillonite mineral may substitute for hydromica either in the whole or in the outer 
part of the innermost zone, thus producing two subpatterns. Phlogopite and siderite repre- 
sent the alteration of biotite; they extend to and suddenly increase at the hydromica 
zone. The remaining zones are based on the alteration of plagioclase. Microcline and 
quartz are stable in all zones. 

Generally, the simplest alteration patterns are found in the southwestern part and 
along the western side of the belt. They are closely related to Laramide stocks, and most 
of the associated mineralization is of lead and zinc sulfides. The most complex patterns 
predominate along the eastern side and are concentrated in the middle to northeastern 
part of the belt; they are associated preferably with tungsten, pyritic gold, and telluride 
gold deposits. 

The origin of the complex patterns is explained in terms of chemical fronts developed 
as the altering fluids migrated away from the conduit. The altering fluids were probably 
gaseous and are believed to have been derived by fractional distillation from the main 
body of hydrothermal liquids. The remaining silica-rich, less volatile fraction filled the 
veins after wall-rock alteration. The fluids changed in composition through reaction with 
the wall rock and successively reached the stability fields of the various secondary miner- 
als. The critical changes responsible for the front (and zone) formation are increase of 
Si/Al ratio, decrease of K concentration, and increase of pH. No evidence of mutual 
replacement among hydrothermal minerals was found. The amount of matter transferred 
during alteration was small. The formation of other patterns is explained on the same 
principles, assuming slightly different compositions of the altering fluids. An attempt is 
made to correlate the patterns of alteration with the distribution of the intrusive bodies. 
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INTRODUCTION 
Scope of Problem 
Studies on wall-rock alteration associated 
with ore veins are important from the stand- 
point of the nature of the hydrothermal proc- 
esses, as well as of their possible use in the 
search for ore. Workers in this field disagree 
with regard to the mechanism of the alteration 
and its significance in the ore-forming processes. 
The present emphasis on the laboratory ap- 
proach to the study of hydrothermal phenomena 
looks fairly promising, but a program of sys- 
tematic field investigation seems necessary. 
Alteration should be studied and compared 
from district to district beginning with those 
with the simplest structural setting. 

The Front Range mineral belt, having various 
types of relatively simple vein deposits as- 
sociated with igneous bodies of diverse compo- 
sition, offers an excellent opportunity for the 
correlation between hydrothermal alteration, 
ore-forming processes, and magmatic differen- 
tiation. An excellent background for such a 
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study is provided by the detailed geologic and 
petrologic work done in this area by several 
workers, in particular Lovering and his as 
sociates (Lovering, 1941; Lovering and God- 
dard, 1938; 1950; Lovering and Tweto, 1953). 
Wright (1954) described the alteration related 
to the uranium-bearing veins of Caribou mine, 
west of Nederland, Colorado, and Tooker 
(1955; 1956) and his colleagues of the U. §$ 
Geological Survey studied the geology and 
hydrothermal alteration in the Central City- 
Idaho Springs area. The results of the latter 
work were not available when this was written 

The exploitation of all the possibilities of- 
fered by the Front Range mineral belt for the 
study of hydrothermal alteration would require 
many years. Only a few of the problems have 
been studied here. Future work calls for inte- 
gration of the results here presented with those 
of others working elsewhere in the same area, 
for further studies of critical features related to 
age and magmatic relationships of the veins 
and intrusives, and for extension of the studies 
to other districts in the Front Range. 
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INTRODUCTION 


Present Study 


Field studies were carried out during parts 
of two summers (1955-1956). Although an 
effort was made to cover the entire mineral 
belt (with the exception of the Central City- 
Idaho Springs district), the major part of the 
work was focused on the northwest part of the 
belt, where alteration is most complex and 
workings most accessible. Only a few mines were 
active throughout the entire belt at the time of 
this study, and these only in part. The majority 
of the veins examined were in idle or abandone:i 
tunnels where most of the ore had been mined. 
The difficulties presented by this area for this 
kind of work were aptly summarized by Tooker 
(1955). 

Samples were studied by optical methods, 
supplemented by X-ray and thermal analyses. 
Most of the petrographic analyses were made by 
immersion techniques and thin-section studies. 
Using properly spaced X-ray analyses as sup- 
porting data a large number of profiles can be 
studied in a short time. For estimating the rela- 
tive proportions of clay minerals, thin sections 
were also used, but if the minerals were too 
intimately mixed, X-ray studies were necessary. 
The mineralogical study of clay minerals and 
other alteration products was handicapped by 
the scarcity of material, which generally had to 
be dug painstakingly out of the pseudomorphs. 
Fortunately, the writer had access for a short 
time to an ultrasonic disperser at the X-Ray 
Laboratory of the U. S. Geological Survey in 
Denver, through the courtesy of A. J. Gude, 
3rd. This apparatus proved to be very effective 
in the quick separation of clays from pseudo- 
morphs, using the minimum amount of crushing 
and therefore reducing contamination. 
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SUMMARY OF REGIONAL GEOLOGY AND ORE 
Deposits 


Lovering and Goddard (1950) gave a detailed 
account of the geology of the Front Range 
mineral belt and its ore deposits. The following 
summary is taken from their monograph. 

The belt extends from Jamestown to Brecken- 
ridge (Fig. 1). The veins are located in igneous, 
metamorphic, and mixed rocks of different 
ages. The oldest rocks are the micaschists and 
hornblende schists of the Idaho Springs forma- 
tion (early Precambrian), with an approximate 
thickness of 20,000 feet. The schistosity strikes 
generally between north and east and dips very 
steeply. Migmatization by granitic material is 
widespread, especially on the eastern end, and 
various forms of gneiss and arterites are found. 
Hornblende schists are predominant in the 
southwest end of the belt. 

Irregular, largely interconnected bodies of 
Precambrian granitic rock cut the Idaho 
Springs formation and form the host rock of the 
majority of the veins. They are most abundant 
in the northeast and central parts of the belt 
and least abundant in the southwest end where 
metamorphic rocks predominate. Two main 
granite types are distinguished in the mineral 
belt: The Boulder Creek granite, in the north- 
eastern part of the belt, is a coarse-grained, 
foliated rock with considerable amount of 
biotite and numerous schist inclusions and 
schlieren. The other type is represented by the 
Silver Plume granite, which predominates 
around and west of Georgetown-Empire. The 
latter is a medium-grained, slightly porphyroid 
biotite granite, which lacks the strong foliation 
of the Boulder Creek granite. Numerous 
pegmatite dikes cut both the granites and the 
metamorphic rocks. 

Outward dipping, Paleozoic and Mesozoic 
sedimentary formations flank the Front Range 
at both ends of the mineral belt and enclose the 
stock to which the Breckenridge veins are re- 
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lated. Tertiary rocks are also included in the 
sedimentary series. 

A series of small igneous bodies of Laramide 
age extends parallel to the mineral belt (‘por- 


the veins are generally richer closer to their 
intersection with the reefs. The mineral belt js 
almost coincident with the line of Laramide 
porphyry stocks, although somewhat to the 
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Figure 1.—Inpex Map or THE Front RANGE MINERAL BELT, CoLorapo, SHOWING 
DISTRIBUTION OF ALTERATION PATTERNS 


Base map after Lovering and Goddard, 1950. 


phyry belt”) on the northwestern side. In the 
mineral belt proper dikes and small irregular 
bodies are abundant. The rock types range 
generally from diorite (and andesite) porphyry 
to quartz monzonite (and latite) porphyry, but 
gabbroic and syenitic bodies are also found. In 
general, the age of the Laramide intrusions 
increases toward the southwest. 

The mineral deposits in the belt are dis- 
tributed in several districts, which, northeast to 
southwest, are: Jamestown, Gold Hill, Ward, 
Caribou, Nederland, North Gilpin County, 
Central City, Idaho Springs, Georgetown, Silver 
Plume, Montezuma, and Breckenridge. They 
are represented by veins of predominantly 
northeast to east trend, formed in Laramide 
times. Many persistent, northwest-trending 
barren breccia zones (‘‘reefs’’) cross the belt: 


southeast. The related origin of igneous bodies 
and veins seems evident. The veins are younger 
than the porphyry and in places cut the Lara- 
mide igneous rocks. 

The Jamestown district is characterized by 
pyritic gold telluride veins, more or less clustered 
around a composite stock of granodiorite-sodic 
granite. Fluorspar veins and lenses are closely 
associated with gold telluride veins. Silver-lead 
veins are also present. Most of the veins are in 
Boulder Creek granite. 

The Gold Hill district (Fig. 1) has pyritic 
gold telluride mineralization. Most of the veins 
are in Boulder Creek granite, some in gneissic 
schists. The district is separated from the Ward 
district by an irregular body of Laramide age. 

The deposits in the Ward area are varied, 
with pyritic gold-silver, copper-gold-silver, 
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lead-silver, and tungsten veins. Only a silver- 
lead deposit (White Raven) was accessible to 
the author. 

The Nederland tungsten district (Lovering 
and Tweto, 1953) is characterized by simple 
quartz-ferberite veins. Most of the veins are 
in Boulder Creek granite and biotite aplite. No 
Laramide intrusions are found here except 
some hornblende monzonite and hornblende 
diorite dikes. 

The Caribou vein system, west of Nederland, 
is partly in a composite, gabbro-syenite stock 
with concentration of titaniferous magnetite. 
The veins contain lead-silver-zinc ore and in 
places also pitchblende. 

The districts in the middle part of the mineral 
belt (Blackhawk, Central City, Idaho Springs, 
and Empire), which have not been included in 
this study, contain a variety of pyritic gold, 
lead-silver-zinc, gold telluride, and uranium 
veins. 

The Georgetown-Silver Plume district con- 
sists of silver-lead-zinc, as well as some pyritic 
gold veins cutting Idaho Springs schists and 
gneisses and Silver Plume granite. The near-by 
Argentine district has similar characteristics. 

The Montezuma district is separated from 
the Argentine district by a body of quartz 
monzonite porphyry. Lead-silver-zinc veins cut 
Idaho Springs hornblende-rich schists. The 
Breckenridge district contains also lead-silver- 
zinc veins, but they cut extensive bodies of 
quartz monzonite and quartz latite porphyries 
which in turn intrude Paleozoic and Mesozoic 
sediments. 

The general sequence of mineralization was as 
follows (Lovering and Goddard, 1950): (1) py- 
rite, (2) sphalerite, (3) chalcopyrite, (4) galena 
and chalcopyrite, (5) silver-bearing sulfoanti- 
monides, sulfarsenides, and _bismuthinides, 
(6) pyrite and subordinate chalcopyrite, (7) free 
gold, (8) minor amounts of sphalerite, galena, 
and silver minerals, (9) gold tellurides and 
sparse pyrite, gold, sphalerite, and galena, and 
(10) ferberite and sparse sulfides. Most of the 
ores between Breckenridge and Idaho Springs 
belong in groups 1-5. The ores to the north 
belong in groups 6-10. This is at least partly in 
agreement with the fact that the latest members 
of the porphyry series are found in the north- 
eastern part of the belt. Southwest of George- 
town, where complex lead-silver-zinc ores pre- 
dominate, the igneous rocks are diorite, mon- 
zonite, quartz monzonite, and sodic quartz 
monzonite. Northeast of Silver Plume, as the 
pyritic gold deposits appear, the igneous bodies 
change to more potassic types (bostonite, 


syenite). The gold telluride veins between 
Idaho Springs and Jamestown are coexistent 
with dikes of biotite latite, biotite monzonite, 
and latite. The pitchblende veins at Central 
City represent a local modification of the 
pyritic gold veins. 


PRE-LARAMIDE ALTERATION 


The igneous rocks of the Front Range 
mineral belt show alteration that is unrelated 
to the veins of the Laramide cycle. The chief 
alteration is sericite in plagioclase. Penninite 
and epidote replace biotite in a few places, and 
a little carbonate of doubtful origin is present. 
The amount of alteration varies within a given 
district, but no attempt was made to map such 
changes. Sericite is moderately abundant, 
amounting in places to about a quarter of the 
plagioclase host. Compared with the hydromica 
associated with ore veins, sericite shows much 
better crystallinity; in places, coalescence of 
subparallel plates has given rise to large indi- 
viduals. The plates are commonly arranged 
along the crystallographic planes of plagioclase. 
Measurements of a few samples gave refractive- 
index values well above 1.590 for gamma, ap- 
proaching the index of well-crystallized musco- 
vite. 

The nature of sericite alteration of this type 
is a problem of general significance but is beyond 
the scope of this paper. Significant features of 
this alteration are the widespread distribution, 
the coarse grain, and the simplicity of the 
mineral assemblage. The first two features 
point toward a highly pervasive, hot fluid; the 
third characteristic indicates lack of dif- 
ferentiation. These conditions are best given 
during the deuteric stage of the consolidation of 
igneous bodies. At this stage, the residues of the 
magmatic crystallization, containing little else 
than water, COs, silica, and alkalies, are dis- 
persed in the rock under considerable vapor 
pressure. In the presence of this fluid, plagio- 
clase is unstable, whereas potassium feldspar is 
generally insoluble. Therefore, sericite replaces 
the former but leaves generally untouched the 
potassium feldspar. This is an example of what 
the writer believes to be a general principle in 
alteration: the localization of a secondary 
mineral is often a question of available space 
rather than of direct chemical contribution from 
the host. 

The interpretation of the sericite replacement 
as deuteric presupposes that most of the 
granitic rocks of the mineral belt are of igneous 
origin. They are thus explained by Lovering 
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and Goddard (1950), and the writer sees no 
evidence to the contrary. 

The recognition of this alteration is im- 
portant for a correct interpretation of the 
later cycles. Failure to recognize and separate 
the pre-Laramide alteration would result in a 
much more confused picture of the later 
phenomena. 


LARAMIDE ALTERATION 
General Statement 


In contrast with the older alteration, the 
Laramide hydrothermal alteration in the 
mineral belt is closely related to the ore veins. 
The rock decomposition is more intense as well 
as more concentrated than in the pre-Laramide 
sericitization. The intensity and, generally, the 
mineralogical nature of the alteration vary 
laterally from the veins. These changes manifest 
themselves in the color and hardness of the 
wall rock, but some mineral variations are not 
detectable by simple observation. All external 
changes correspond to mineral variations, but 
the reverse is not true. 

The Laramide alteration is not everywhere 
restricted to the individual veins. In some de- 
posits (e.g., in the Georgetown area) a wide- 
spread alteration comprises several veins or 
perhaps 4 whole vein system. No attempt was 
made to determine the spatial relationships in 
the deposits. The present study deals with the 
alteration that can be correlated with individual 
veins, which give a better opportunity to study 
the critical genetic features. 

The zoning, i.e., the systematic lateral 
distribution of the mineral assemblages—varies 
from simple (one zone) to complex (up to five 
zones, excluding the fringe zone), each zone 
characterized by a given mineral assemblage. 
Each combination of zones is here called a 
“pattern”. 


Zone Patterns 


According to the lateral distribution of the 
alteration minerals, the veins of the Front 
Range mineral belt studied by the writer can 
be divided into six patterns of zoning (Fig. 2): 

PATTERN I: The simplest pattern is a single 
zone of bleached “sericitized”? wall rock, in 
places silicified. It is found at the Wellington 
(Breckenridge) and Silver King (Montezuma) 
mines (Pb-Ag-Zn) in the southwestern part of 
the belt, and at the Oregon-Quaker City mine, 
Nederland (tungsten). 
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PATTERN I1: A single zone characterized bya 
mixture of kaolinite, hydromica, and 
montmorillonite in variable proportions, 
Georgetown and Ward districts (Pb-Ag-Zn- 
Au-pyrite). 

PATTERN 111: A single zone with predominant 
hydromica and potassium feldspar. Jamestown 
district (Au-Te-pyrite). 

PATTERN IV: Two separate zones, character. 
ized respectively by kaolinite (inner) and 
montmorillonite. Caribou mine, Nederland 
(Pb-Ag-Zn). 

PATTERN V: Three zones, characterized te. 
spectively by hydromica (inner), kaolinite, and 
montmorillonite. Found sporadically within the 
area of pattern VI. 

PATTERN VI: The most complex pattem 
consists of four main zones: hydromica (inner), 
orthoclase, kaolinite, and montmorillonite. It 
is typical of Nederland tungsten-Gold Hill 
pyritic gold districts. Two subpatterns are 
recognized in a few places, as the regular mixed- 
layer clay mineral described above substitutes, 
either partly or totally, for the hydromica in the 
inner zone. 

The above classification is based on the 
alternation of plagioclase. Further subdivisions 
of zones and/or patterns may be made by 
considering also the alteration of mafic con- 
stituents. For instance, a fringe zone, char- 
acterized by the partial alteration of biotite 
and (in places) hornblende is distinguished in 
most veins. A zone of complete alteration of 
biotite is superimposed on the hydromica zone 
in nearly all veins of the most complex pat- 
terns. 


Mineralogy 


The prominent primary constituents of the 
rocks of the mineral belt are plagioclase, quartz, 
potassium feldspar in most places, and biotite. 
Hornblende and muscovite average less than 
1 or 2 per cent each in the mines examined, 
whereas pyroxene was found at only one locality 
(Caribou mine). 

Plagioclase is unstable in all but the fringe 
zone. The degree of alteration is low in the 
montmorillonite zone of patterns V and VI 
but reaches 100 per cent between there and the 
vein. Quartz is essentially stable in all zones; 
microscopic textures suggesting peripheral re- 
placement by clay aggregates are seen in the 
kaolinite zones at a few mines, but modal and 
chemical analyses do not support the idea of 
quartz replacement to any measurable extent. 
Potassium feldspar is stable in all zones. Biotite 
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» 
is generally unstable, but the degree of altera- 
tion varies considerably. Its replacement is 
complete only in the innermost zone, where it 


linite is more abundant than second and is 
strictly related to the alteration zoning; it 
appears in white, soft aggregates replacing 
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FictrE 2.—PATTERNS OF HYDROTHERMAL ALTERATION IN THE FRONT RANGE MINERAL 
BELT, COLORADO 


has the character of a cationic substitution 
within the mica framework (phlogopitization). 
Hornblende and pyroxene are also unstable in 
all zones; however, in the fringe zone the 
alteration of the amphibole may be incomplete. 
Muscovite, as well as apatite, is stable every- 
where. ~ 

The main alteration minerals, representing 
more than 99 per cent of most of the hydro- 
thermal products, are described below. 

KAOLINITE: Kaolinite is the only representa- 
tive of its group that the writer found in associa- 
tion with the veins of the mineral belt. Two 
generations are present. First-generation kao- 


plagioclase and, to a lesser extent, biotite, 
hornblende, and microcline. Most individuals 
are tiny booklets, but a few show hexagonal 
outlines under the highest power. Vermicular 
development is uncommon. Most of the 
kaolinite is mixed with montmorillonite 
stringers that separate the booklets. 
Second-generation kaolinite is found in 
veinlets in which quartz lines the walls and 
kaolinite fills the middle part. Opal may be 
found intimately mixed with the clay. This 
kaolinite is in general coarser-grained than the 
first generation and free from montmorillonite. 
Kaolinite that replaces biotite and horn- 
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blende differs in habit from the type that 
replaces plagioclase: it forms white, fibrous 
aggregates. Montmorillonite (or nontronite) 
may be mixed in the aggregates. In some ag- 


TABLE 1.—(ool) SPACINGS IN MONTMORILLONITES 
AND NONTRONITE 


1 2 4 
Poorman | _— Caribou | — 
| Raven | 


Spring 


a(A)| ¥ | a(A) | 1 {a(A)| 1 | aA)| I 


17.1 |10 /16.7 110 16.4 10 





001 |16.9 10 | | 

002 | 8.42) 1 | 86 | 1/83) 1 | 8.4) 3 
003 | 5.6211 | 5.72} 1| 5.641 | 5.56 1 
004 | 4.23) 1 | 4.261) 4.231 ? 
005 | 3.38) 1 | 3.44) 1) 3.34) 1 |? 
006 | 2.87) 0.5, 2.86 | 1 | 2.84 0.5! 2.69 1 





(1-3) Montmorillonite; (4) Nontronite; all 
samples glycolated. 


gregates of this type (Chance mine) these 
fibers are seen to be subhedral plates 10-20 
times longer than wide, elongated in the direc- 
tion of one of the edges of the pseudohexagonal 
basal pinacoids. Thin montmorillonite films 
separate the kaolinite laths. 
One sample of kaolinite gave the following 
indices of refraction: a = 1.561, y = 1.567. 
No dickite, halloysite, or any other member 
of the group was found in spite of extensive 
microscopic and X-ray investigations. 
MONTMORILLONITE: Montmorillonite is nearly 
everywhere mixed with kaolinite to form 
pseudomorphs chiefly after plagioclase. Where 
montmorillonite predominates, the aggregates 
are light gray; if kaolinite is predominant they 
are white. Veinlets of pure montmorillonite are 
found in largely unaltered plagioclase in the 
marginal zone of alteration of some veins (PI. 2, 
fig. 6). They form a network of stringers and 
threads enclosing the kaolinite individuals. 
The gamma index of refraction measured in 
several samples has a value around 1.530, 
although values as low as 1.513 were found in 
samples from the Prospect tunnel. Montmoril- 
lonite about 95 per cent pure, with a little 
kaolinite, was recovered by ultrasonic dispersion 
of samples from Poorman, White Raven, and 
Caribou mines. X-ray-difiraction spacings of 
these clays are shown in Table 1. The dry and 
glycolated dg; spacings of the H-, Na-, NH4-, 
and Ca-saturated montmorillonite from Poor- 
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man and *iite Raven mines are shown on 
Figure 7. The diffraction patterns suggest that 
montmorillonite varies little throughout the 
mineral belt. 

NONTRONITE: Nontronite is found, together 
with kaolinite and siderite, replacing horm- 
blende. It is pale greenish brown, and the 
indices of refraction are: a = 1.545; y = 1.580, 
variable. Although comb texture is found in 
some aggregates, nontronite generally forms 
films or screens separating kaolinite aggregates 
along planes parallel to the original cleavage of 
the hornblende. Table 1 shows the X-ray 
basal reflections of an impure sample obtained 
from a pseudomorph after hornblende. Very 
small amounts of this clay mineral were avail- 
able for study. 

HYDROMICA: The term hydromica is used 
here, following Brown and McEwan (1951), to 
include illites as well as interstratified mixtures 
of illite and other clay minerals in which the 
former predominates. 

Hydromica is the most abundant alteration 


mineral in the Front Range mineral belt. It | 


appears in aggregates of tiny, ill-defined, leafy 
booklets and in more elongated individuals ar- 
ranged approximately parallel to the twinning 
planes of the plagioclase they usually replace 
(Pl. 2, fig. 1). Most hydromica can be dis- 


tinguished from the earlier sericite, which shows | 


better defined plates. There are, however, many 
indistinguishable border cases. Hydromica 
forms grayish or creamy-white, locally iron- 
stained aggregates, pseudomorphic after plagio- 
clase. Creamy or yellowish tinges indicate 
phosphate. 

The argillaceous (as opposed to micaceous) 
nature of the mineral forming the bulk of the 
secondary minerals of the inner alteration zone 
in most of the veins is substantiated by optical, 
X-ray, thermal, and chemical analyses. The 
results prove that the mineral is a hydromica, 
not a sericite. Sericite of an earlier period of 
alteration is present in many aggregates. 

Optically, the hydromica shows relatively 
low refraction indices; the majority have 
a = 1.545 + 3,y = 1.580 + 3. The range of 
variation for gamma is 1.585 to 1.570, and the 
alpha varies accordingly. Associated sericite 
has a gamma index of about 1.590. The above 
measurements were made on untreated flakes. 

X-ray-diffraction patterns (Fig. 3) indicate 
some degree of disorder in the stacking of 
layers. Moreover, the (001) reflections show 4 
relatively low spacing, generally between 
9.9 A and 9.8 A, and in one case 9.7 A was 
recorded (Fig. 3). This low spacing is attributed 
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to the random interstratification of montmoril- 
lonite within the illite, in a proportion of 10 to 
15 per cent in the common case, and as much as 
25 per cent in the case of the lowest spacing, 
9.7 A. These proportions are obtained by means 
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The curve of Figure 4 may be used to ascer- 
tain the approximate proportions of illite and 
montmorillonite in the mixed-layer clays of the 
mineral belt. It could presumably be used also 
in other areas, but the nature of components 
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Figure 3.—X-RaAy-DiFFRACTION CURVE OF HypROMICA CONTAINING 25 PER CENT OF 
RANDOMLY INTERLAYERED MONTMORILLONITE 


; of the curves calculated by Brown and McEwan 


(1951); see Weaver, 1956). The slight but 
significant displacement of the (001) reflection is 
due to the influence of the (002) reflection of 
montmorillonite, and it is best observed after 
glycolation. As the amount of randomly inter- 
layered montmorillonite increases, the (001) 
illite peak is displaced, and the (001) mont- 
morillonite peak answers (Fig. 3). 

The comparatively low indices of refraction 
of the hydromica of the Front Range mineral 
belt are also attributed to the presence of 
interlayered montmorillonite. Figure 4 shows 
the y indices of three samples plotted against 
the percentage of montmorillonite as deter- 
mined by means of Brown and McEwan’s 
curves. Sample 1 represents the common 
hydromica in the area. Sample 2 is the hydro- 
mica with the highest content of interlayered 
montmorillonite (Cold Spring mine). Sample 3 
corresponds to a regularly interlayered, 1:1 
illite-montmorillonite clay mineral described 
below. A straight line traced through these 
three points intersects the ordinates of pure 
illite and montmorillonite at values of 1.590 
and 1.513, respectively. Such indices of re- 
fraction fit well with the standard values for 
these minerals. 


must be established first. The presence of 
interlayered montmorillonite in illite may be 
recognized or suspected by a decrease in the 
indices without perceptible change in bire- 
fringence. Interlayered kaolinite would diminish 
the birefringence with relatively little effect on 
refractivity, except where present in large 
amounts. Chlorite would be indicated by re- 
duced birefringence and by coloration. Py- 
rophyllite would be impossible to detect. 

The displacement of the (001) peak of illite 
by montmorillonite interlayering is enhanced 
by glycolation; this is due not to the expansion 
but to the intensification of the second-order 
peak of montmorillonite, which is very weak in 
the air-dried state. 

Differential thermal analysis of an Ill;Mtis 
clay mineral (Fig. 5) shows the mid-tempera- 
ture endothermic peaks of the two components 
to be about 600°C. and 700°C., respectively. 
The X-ray pattern of the same clay indicates 
almost pure mixed-layer hydromica. The 
thermal curve is practically identical to that of 
the clay from Sarospatak, Hungary, which is 
also a mixed-layer hydromica (Grim and 
Rowland, 1942). The name “bravaisite’” has 
been suggested for this type of clay mineral. 

The chemical analysis of a typical hydromica 
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(Ingram mine) is shown in Table 2,a. The 
sample was purified by sedimentation but still 
contains minor amounts of quartz (5 per cent), 
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FIGURE 4.—RELATION BETWEEN GAMMA RE- 
FRACTION INDEX AND MONTMORILLONITE CONTENT 
IN Mrxep-Layer ILLITE-MONTMORILLONITE CLAY 
MINERALS 

1, “theoretical’’ illite; 2, hydromica, Hummer 
vein; 3, hydromica, Cold Spring mine; 4, 1:1 














illite-montmorillonite clay, Chance mine; 5, 
‘theoretical’? montmorillonite. 
T T T 
Hummer Hydromica 
(1) 
Clark Regular mixed-layer 
(2) 
Chance Regular mixed-layer 
200° 400° 600° 800°C. 
1 1 rn i rn l n 





FiGuRE 5.—DIFFERENTIAL THERMAL CURVES 
OF Mixep-LAvYER, ILLITE-MONTMORILLONITE 
CLay MINERALS 

1, hydromica with 13 per cent montmorillonite; 
2 and 3, 1:1 illite-montmorillonite. The presence of 
about 10 per cent kaolinite enhances the 600°C. 
peak in the last two samples. 


pyrite (3.2 per cent), phosphates (2 per cent), 
and gypsum (.8 per cent). With the exception 
of quartz, these impurities can be reckoned out 
of the analysis. The percentage of quartz was 
found by X-ray methods, using the following 
procedure: The relative intensities of the 4.45 A 
(clay) and 4.25 A (quartz) peaks were deter- 
mined for samples containing increasing 


amounts of added quartz. A curve was thus 
obtained which by extrapolation gave a quartz 
content of 5 (+1) per cent for the original 
sample. The structural formula shown below, 


TABLE 2.—CuHEMICAL ANALYSES OF HyprRomic, 
AND MIxEp-LAYER CLAy MINERAL 
Analyst: Doris Thaemlitz, Rock Analysis 
Laboratory, University of Minnesota 


(a) (b) 
Ingram FR112 | Chance FR156 


53.13 


SiO» 53.74 
AlO; 23.78 27.02 
TiO, .76 .03 
Fe.O;* 3.87 1.01 
MnO .02 .00 
MgO 2.02 1.30 
CaO ao li 
NasO «85 06 
K:0 6.96 4.96 
BaO 02 .02 
H.O+ 3.94 6.83 
H.0— Ruae 3.66 
SO; .89 n.d. 
S 1.70 .08 
P.O: .34 08 
99.97 98.29 


* Total Fe as Fe.O; 

(a) Hydromica, including about 10 per cent of 
randomly interstratified montmorillonite 

(b) Regular mixed-layer 1:1. illite-montmoril- 
lonite mineral 


caiculated from the corrected analysis, shows 
that the charges are distributed between the 
tetrahedral and octahedral layers: 


Ky. 32Na_.os(Siz.2g Ti. o7Al 67) (Als 37FeisMg 45) Ox9(OH), 


This formula does not reflect the real structure 
because, according to X-ray data, this mineral 
contains about 10 per cent interlayed mont- 
morillonite. 

REGULAR MIXED-LAYER CLAY MINERALS: The 
existence of clay minerals formed by regularly 
interstratified unit layers of two different 
species has been recognized only recently. 
Alternating vermiculite and pyrophyllite (recto- 
rite) wat described by Bradley (1950). Vermic- 
ulite or montmorillonite interstratified with 
chlorite (corrensite) has been found by Stephen 
and McEwan (in Bradley and Weaver, 1956), 
Honeyborne (in Bradley and Weaver, 1956), 
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Lippmann (1954), Sudo et al. (1954), Early et 
il, (1956), and Bradley and Weaver (1956). 
[llitemontmorillonite has been found by 
Heystek (1954). In all instances the ratio of 
the two constituents is 1:1. 
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vealed by its X-ray-diffraction pattern. Figure 
6 shows the diffraction curves of three mixed- 
layer clays from different localities, in air-dried 
and glycolated conditions. The mounts are 
oriented aggregates, so that only the basal 


TABLE 3.—TABULATION OF LATTICE BASAL SPACINGS OF REGULAR MIXED-LAYER 
ILLITE-MONTMORILLONITE CLAY MINERALS 
See Figure 6 for relative intensities. 
(In A Units) 











3.01 


Examples of the illite-montmorillonite mixed- 
layer mineral were found in the course of the 
present study in the Chance (Tungsten), Clark, 
Crackerjack, Prospect, and Lily mines. The 
mixed-layer clay occupies a definite position 
within the alteration halo of the veins and plays 
essentially the same role as hydromica (see 
below). Where reasonably pure, it forms gray- 
white aggregates pseudomorphic after plagio- 
clase. Under the microscope the aggregates 
resemble those of hydromica although they are 
somewhat finer-grained. These aggregates are 
much more easily dispersed than those of hy- 
dromica. In some aggregates the ill-defined, 
tiny (5-10 microns) booklets of mixed-layer 
clay show interstratification of a low-index 
substance. X-ray analysis failed to reveal the 
nature of these intercalations. The refraction 
indices are fairly consistent: 


@=1520+2, y = 155242, B = 0.032. 


The angle 2V(—), as measured on oriented ag- 
gregates, is somewhat smaller than that of 
muscovite. Possibly randomness in the a, 0 
orientation of the plates has reduced the ap- 
parent optic angle. The indices of refraction 
are the arithmetic mean of the two interlayered 
clays, as shown in Figure 4. 

The interlayered nature of the clay is re- 





Chance | Crackerjack Clark | Lily Prospect 
| 
Onder Air- Glyco- | Air- Glyco- Air- | Glyco- Air- | Glyco- Air- Glyco- 
dried lated | dried lated dried lated dried | lated dried lated 
- = — = = mnt en = = 
001 | 23.10] 26.8 | 23.5 26.3 24.5 | 26.8 22.3 | 26.0 24.5 26.0 
002 11.55 13.4 11.9 13.2 11.5 13.2 11.5 | 13.2 12.0 13.3 
003 7.87 8.94 | 7.5 8.93 | 9.0? — | 8.4 — — 
04 5.46, 6.73) 5.95 6.6 6.6 — | - 
005 4.70 5.37 | 4.97 5.30 4.97 5.2 5.3 4:95) 3:2 
006 3.93 4.49) — 4.44 4.40 
007 3.36 3.88 3.30 g.49 3.78 | 
008 _- 3. | 3.13) 3.33 3.30. 
09 == (2.62 2.98 2.99 





reflections are recorded. The corresponding 
spacings are shown in Table 2. The significant 
feature of all patterns is the presence of a 001 
periodicity between 23 A and 24 A in the air- 
dried state and between 26 A and 28 A in the 
glycolated state. This large spacing is due to the 
alternation of illite and montmorillonite layers 
in a 1:1 ratio. The glycolation increases the 
basal spacing by about 3 A. Furthermore, 
glycolation reinforces the basal reilection (es- 
pecially the weak higher orders) and makes the 
series more rational. In the Chance and Cracker- 
jack clays the 003 (8.94 A) reflection appears 
after glycolation. This is explained by its partial 
reinforcement by the 001 reflection of illite 
(10 A) and the 002 reflection of montmorillonite 
(8.2 A-8.4 A). The presence of free hydromica 
as impurity, as in the samples from Clark and 
Lily tunnels, masks this reflection. 

The differences in the 001 spacings depend 
largely upon the exchangeable cations. The 
result of leaching in Chance clay with different 
cations can be seen in Figure 7, where the 001 
spacings of the air-dried and glycolated H-, 
Na-, NH4-, and Ca-saturated clays are plotted 
together with those of two montmorillonite 
clays. To facilitate the comparison the spacing 
corresponding to the illite layer (10 A) has been 
subtracted from the 001 spacing of the mixed- 
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layer clay. The coincidence in the plotting is 
further proof of the presence of montmorillonite 


layers in the structure. 
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temperature, exothermic peak shown by the 
Chance sample is probably reinforced br 
kaolinite. ; 
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FiGuRE 6.—X-Ray-DIFFRACTION CURVES OF REGULAR Mrxep-Layer, 1:1 
ILLITE-MONTMORILLONITE CLAY MINERALS 
Some peaks have been reduced in a ratio indicated in each case. Glycolated (top) and 
air-dried (bottom) samples used in each pair. 


The thermal behavior of the mixed-layer 
clay is demonstrated by the D.T.A. curves oi 
Figure 5. They show a low-temperature endo- 
thermic peak with a shoulder (indicating that 
Ca as the predominant exchangeable cation) 
and two middle endothermic peaks near 600°C. 
and 700°C. The last two peaks correspond, 
respectively, to the illite and montmorillonite 
lavers. The 600°C. peak appears to be rein- 
jorced by the kaolinite (Chance) and hydro- 
mica (Clark) present as impurities. The high- 


Table 2, > presents a chemical analysis of the 
mixed-layer clay from the Chance mine. The 
sample, which contains about 2 per cent kaoli- 
nite, has been HCl-leached. The structural 
iormula, after correcting for kaolinite and a 
small amount of pyrite and woodhouseite, is as 
iollows: 


K s9(Al.zoSiz 30) (Fe. 3} Mg soAls, ¢2)020(OH): 


Ii the formula is split into a 1:1 illite and 
montmorillonite layering it becomes: 
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Illite: Ky.7s( Al. s0Sis. 60) (Mg.ssAls.¢2)O20(OH)s 
M..22 


Montmorillonite: Sig,oo(Fe®s2 Mg.22 Als.62)O20(OH4) 


The above formulas probably do not repre- 
sent the actual distribution of atoms in the 
layers, Which is possibly as follows: 


Sirs + Als 
Or 

| Kis 

| Os 

| Sis-x + Als 
|O,; + (OH) 


Aly + (Mg, Fe*, Fe*)y (y = 0-4) 


| Exchangeable cations + water 
p 

| Sig 

|O; + (OH)2 

| Aly-y + (Mg, Fe*, Fe?*)y 

|O, + (OH)2 

| Six + Als 

Os 

Ki-2 


PHLOGOPITE: The magnesian mica phlogopite 
results from the alteration of biotite. It appears 
as monocrystalline pseudomorphs generally in- 
cluding abundant siderite or ankerite (Pl. 2, 
fig. 2). Unaltered remnants of biotite parallel 
to the cleavage planes are common. The pseudo- 
morphs are light green, darker where the iron- 
tich carbonate is abundant, and yellowish where 
jarosite is also present. Under the microscope 
phlogopite appears pale green to colorless and 
shows numerous inclusions of rutile needles, 
tiny euhedral tablets of epidote, and iron oxides. 
The rutile needles lie in hexagonal pattern on 
the cleavage planes. The gamma index ranges 
irom 1.560 to 1.585, with the majority between 
1.565 and 1.570. The birefringence is about 
0,035, and the optic negative angle is about 
10°. Some single plates show changes in 2V 
values from one area to another, from 0° to 
45°, indicating different degrees of phlogopitiza- 
tion. X-ray analyses show a 060 reflection at 
1.54 A, thus confirming the trioctahedral 
character of this mica. The above data indi- 
cate a very low Fe: Mg ratio. 

BIOTITE: Secondary biotite is found only in 
the Jamestown area. It appears in fine-grained 
clusters of relatively well-formed, sharp indi- 
viduals with good pleochroism. In one case 
a bluish tint (indicating Cr?) was detected on 


Z. Pale biotite with weak pleochroism is also 
found. 

Secondary biotite is not related to the vein 
alteration but belongs to an earlier Laramide, 
greisenlike replacement (see below). 
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Figure 7.—RELATION BETWEEN ool SPACING 
AND NATURE OF THE EXCHANGE CATION 
IN MONTMORILLONITE 
Triangles, Poorman mine; circles, White Raven 
mine; squares, Chance mine. The Chance sample 
is a 1:1 illite-montmorillonite; plotted value is 
doo: minus 10 A. 


POTASSIUM FELDSPAR: Orthoclase is found as 
replacement of plagioclase, generally associated 
with kaolinite and hydromica (Pl. 2, figs. 3 
and 4). Whether or not feldspar is the pre- 
dominant secondary mineral in the pseudo- 
morphs after plagioclase, it preserves the 
crystallographic orientation of the host, as 
indicated by external form or ghost twinning. 
Orthoclase is also found lining microscopic 
fissures, especially where these cut microcline 
crystals. Tiny, euhedral crystals are found 
within a few kaolinite aggregates. Universal 
stage measurements show that the optic plane 
is parallel to (101), and gamma is normal to 
(010). The indices of refraction are somewhat 
variable: a = 1.516 — 1.519; y = 1.523 — 
1.527; 2V(—) = 40° — 42°. 

WOODHOUSEITE: A phosphate of the svan- 
bergite-goyazite group was mentioned by 
Lovering in his work on the Boulder County 
tungsten veins (Lovering, 1941, p. 249). It 
appears as pale-yellow, tiny cuboid crystals 
disseminated especially in hydromica but also 
in the vein-filling microcrystalline (“horn’’) 
quartz. Clay aggregates containing the alunite- 
like phosphate may show yellowish tint. The 
refraction indices on some of the largest rhombo- 
hedrons were: a = 1.646;y = 1.654: B = 0.008. 
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These values fit well with those of the mineral 
woodhouseite. 

JAROSITE: Jarosite is found in small quanti- 
ties, chiefly in the innermost (hydromica) zone 
of alteration, both as an alteration product of 
biotite and as a superficial coating. In the first 
case it is probably hydrothermal, since it shows 
the same type of occurrence as the iron car- 
bonates; the superficial jarosite was probably 
reprecipitated at a supergene stage. 

QUARTZ AND OPAL: Outside the vein filling, 
quartz is present in small to moderate amounts 
in the wall rock as replacing mineral. At the 
Wellington mine, Breckenridge, vein no. 2 
shows considerable silicification accompanying 
the hydromica replacement. Small quantities of 
quartz and opal associate with the kaolinite 
replacing plagioclase in the Cold Spring mine; 
quartz occurs in small, equant grains with 
marginal intergrowth with the clay, whereas 
opal forms tiny vesicles. 


Description of Individual Patterns 


PATTERN I: Pattern I was found in the south- 
western end of the belt, in the only mines active 
at the time of this study—the Wellington mine 
at Breckenridge and the Silver King mine at 
Montezuma. It also extends to the Argentine 
district, where it has been recognized in the 
Grizzly mine. The ores are essentially lead and 
zinc sulfides in carbonate gangue. At Welling- 
ton, pyrite and quartz are also abundant. The 
country rock in the Wellington mine is a quartz 
latite porphyry with abundant, small andesine 
and biotite phenocrysts in a_ fine-grained 
groundmass of quartz, potassium feldspar, and 
little orthoclase. Pre-Laramide (?) alteration 
is represented by some sericite and chlorite re- 
placement of, respectively, plagioclase and 
biotite. At Silver King, the country rock is 
quartz-oligoandesine amphibolite schist with 
layers of fine-grained tonalite and granodiorite, 
probably metasomatic; coarser, more granitic 
areas are also present. 

The Laramide alteration is basically the same 
in both mines. There is a single, symmetrical 
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bleached zone, where plagioclase has been 
replaced by hydromica to a distance of about 
15 times the width of the vein. Iron-rich carbo. 
nate is a common associate of hydromica, 
Biotite is altered to kaolinite, quartz, and sider. 
ite in the porphyry, and to phlogopite and 
siderite in the schist. These assemblages are 
typical of the entire mineral belt. Silicification 
is common but irregularly distributed at Wel. 
lington; it attacked only the groundmass, which 
otherwise is generally less altered than the 
phenocrysts. The bleached zone grades into the 
fresh rock as hydromica disappears. Kaolinite 
is present erratically in the outer edge of the 
hydromica zone at Wellington and is also mixed 
in small quantities with the hydromica in the 
pseudomorphs after plagioclase. Nontronite 
occurs in the border zone at Silver King. Pyrite 
is widely though unevenly disseminated in the 
bleached zone in both mines, and iron carbon- 
ates replaced hornblende in the amphibolite 
schist. 

Outside the southwest end of the belt, pattem 
I was found only at the Quaker City vein in 
Gordon Gulch, near Nederland. A_ profile 
studied just west of the caved crosscut tunnel 
showed a single zone, extending as much as 
100 feet from the vein, where the predominant 
alteration product is hydromica. The plagio- 
clase is not wholly altered throughout the zone, 
and the granite texture is well preserved, but 
in a zone about 2 feet wide enclosing the vein 
the wall rock has been replaced by white-gray, 
fine-grained, moderately soft material composed 
of quartz, hydromica, and sericite. The forma- 
tion of this zone seems to have been favored by 
shearing and brecciation of the wall rock. 

PATTERN I1: Pattern II zoning is predominant 
in the Georgetown-Silver Plume-Argentine and 
the Ward lead-zinc-silver districts. The wall 
rock is injected schist, gneiss, and amphibolite 
of the Idaho Springs formation (Georgetown 
and Argentine), Silver Plume granite (Ward), 
and Laramide quartz monzonite (Ward). The 
alteration here is not so clearly related to the 
veins as in the other patterns and shows as an 
intensive replacement of plagioclase in wide 


Pirate 1.—VEINS SHOWING PATTERN VI OF HYDROTHERMAL ALTERATION 


Ficure 1.—Pyritic gold vein. Gardner tunnel, I 


mancipation mine, Gold Hill district. V, vein; hm 


hydromica zone; or, orthoclase zone; ka, kaolinite zone. Iron staining between hydromica and _ orthoclase 
zones. mf, montmorillonite and fringe zones (not distinguishable on the picture). Kaolinite zone is partly 
covered by dust coating. Width of alteration profile, about 2 feet. 

FicurE 2.—Rock specimen showing alteration zoning. Vasco 6 mine, Nederland (tungsten). V, veil 
(quartz); HM, hydromica zone; Or, orthoclase zone; Ka, kaolinite zone; Ph, zone of total phlogopitization 
of biotite, occupying the inner part of hydromica zone. Note absence of ferruginous bands. Height of speci 


men, 4 inches. 
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VEINS SHOWING PATTERN VI OF HYDROTHERMAL ALTERATION 
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areas enclosing one or more veins. The re- 
placing material is an intimate mixture of 
hydromica, kaolinite, and montmorillonite. 
Biotite appears only slightly bleached, and 
carbonate is abundant. In the White Raven 
mine, near Ward, the sulfide veins show no con- 
spicuous alteration halo, whereas one un- 
mineralized fissure was found with well-de- 
veloped alteration zoning: a main, inner zone 
with the clay mixture, and a less altered, grayer 
outer zone with kaolinite replacing plagioclase 
phenocrysts of the quartz monzonite porphyry, 
and montmorillonite predominating in the 
groundmass. 

East of Ward the vein alteration changes 
gradually into that of the Gold Hill type. About 
2 miles southeast of Ward, on Lefthand Canyon 
road, an abandoned tunnel shows both types 
of veins—that is, sulfide veins surrounded by 
an indefinite, weak hydromica-kaolinite altera- 
tion halo, and pyrite-quartz veins showing 
pattern V or VI, recognized by the ferruginous 
bands. 

PATTERN II: Pattern III, found in the 
Jamestown district, includes four types of 
veins (Lovering and Goddard, 1950, p. 260- 
265): lead-silver, fluorspar, pyritic gold, and 
telluride, in order of decreasing age. No lead- 
silver veins were sampled by the writer. Gen- 
erally, the alteration consists of a single zone 
of hydromica, quartz, and orthoclase, in places 
superimposed on a mild greisenlike replacement 
of the granite, characterized by poikilitic 
muscovite and abundant secondary orthoclase 
and quartz. Fine-grained, well-crystallized 
biotite is common, generally in veinlets asso- 
ciated with muscovite, quartz, and other 
secondary minerals. Secondary is distinguished 
from primary orthoclase—which it partly re- 
places—by its finer grain and association with 
other secondary minerals and by its smaller 


optic angle (40° as against 60°-70°). Fluorite 
is a common associate, as well as pyrite; topaz 
is uncommon. The greisenlike alteration is 
found all over the Jamestown area and is ap- 
parently related to the fluoride mineralization. 
It is not restricted to the area of productive 
fluorspar veins, having been found as far east 
as Springdale (Gladiator mine). The quartz 
monzonite porphyry of the Jamestown stock 
also shows replacement by orthoclase (as 
patches within the plagioclase), fluorite, hydro- 
mica or sericite, and fine-grained biotite. 
Secondary biotite has replaced primary biotite 
and hornblende, or forms veinlike aggregates. 

The wall rock of the fluorspar veins shows 
little apparent alteration, though the micro- 
scope reveals more secondary orthoclase and 
hydromica than the average. Hydromica and 
sericite or hydromica are abundant in the fine- 
grained fluorite veins. 

The alteration in the Jamestown district can 
be summarized as follows: (1) a widespread, 
greisenlike replacement of the granite, repre- 
sented by poikilitic white mica, milky quartz, 
orthoclase, a little fluorite, fine-grained biotite, 
and in a few places topaz; (2) a restricted al- 
teration zone of hydromica and orthoclase 
replacing plagioclase, related mainly to the 
gold telluride pyritic veins; orthoclase and 
fine-grained fluorite aggregates replaces the 
rock almost indiscriminately in places. Veins 
with a zone of almost exclusively hydromica 
(or sericite) alteration are also found (Grand 
Central telluride mine). 

The Jamestown district deserves more de- 
tailed studies to establish more precisely the 
structural and petrographic relations between 
the different phases of alteration. 

PATTERN Iv: The alteration associated with 
the galena-sphalerite-pyrite veins in the Idaho 
tunnel of the Caribou mine (Smith, 1938) near 





PitaTE 2.—MICROSCOPIC FEATURES OF HYDROTHERMAL ALTERATION 
FicurE 1.—Hydromica pseudomorphs after plagioclase. Note twinning reflected in the arrangement 
of clay individuals. Unaltered individuals are quartz. Snowbound mine. Crossed nicols, < 110. 
FicurE 2.—Phlogopite (Ph) and iron carbonate (Si) replacement of biotite. Adjacent plagioclase is 


replaced by hydromica. Logan mine. X80. 


FicurE 3.—Orthoclase (light gray, Or) and kaolinite (dark-gray aggregates) in pseudomorph after 
plagioclase. Some hydromica marks the original plagioclase twinning. Q, quartz. Eureka mine (tungsten). 


Crossed nicols, X 100. 


Ficure 4.—Orthoclase (light gray, Or) and kaolinite (dark-gray aggregates) pseudomorphic after plagio- 





clase. Orthoclase veins, containing some hydromica, are optically iso-oriented. Crevices in the middle of 
veins suggest that replacement progressed from cracks in the host. Eureka mine (tungsten). Crossed nicols, 
X110. 


Ficure 5.—Kaolinite (dark-gray aggregate) and hydromica (light gray) replacement of plagioclase. 
Note relict twinning in hydromica. Q, quartz inclusion. Clay-rich area in the orthoclase zone. Eureka 
mine. Crossed nicols, X 100. 

oo 6.—Montmorillonite veinlets in plagioclase. Montmorillonite zone, Cold Spring mine. Crossed 
nicols, XK 130. 
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Nederland, Boulder County, represents pattern 
IV. This was the only working to which the 
writer had access, since the mine was closed 
down at the time of his visit. The alteration 
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lets. The kaolinite pseudomorphs contain 
subordinate hydromica and sericite; some ¢p. 
close also thin veinlets of potassium feldspar 
The thickness of the zone averages 8 inches, 
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Ficure 8.—Cross SECTION OF ALTERATION ZONES AT CARIBOU MINE, IDAHO TUNNEL 
Pattern IV. HM, hydromica; Ka, kaolinite; Mt, montmorillonite; Q, quartz; V, vein; FR, fringe zone. 
Total width of profile, 3 feet. Curves in upper part of figure show approximate proportions of secondary 


minerals. 


here consists of a bleached, moderately co- 
herent, kaolinite inner zone, and a narrower, 
crumbly montmorillonite-kaolinite outer zone 
(Fig. 8). 

Two veins, the Caribou and the Poorman, 
were sampled at several places. The wall rock 
(Wright, 1954) is a fine-grained monzonite to 
quartz monzonite with variable proportions 
of biotite, hornblende, and—in lesser amounts— 
augite. In the area sampled, the total propor- 
tion of mafic minerals varies widely within a 
few feet or inches. In the monzonite stock, the 
rock outside the range of the Laramide altera- 
tion is very fresh. 

The veins consist of brecciated rock re- 
placed by microcrystalline quartz and iron- 
rich carbonate. Sulfides are discontinuous along 
the veins. Banded texture is common, generally 
with microcrystalline quartz and carbonate in 
the middle and coarse-grained ankerite on the 
sides. 

The bleached, inner zone is characterized by 
kaolinite pseudomorphs after _ plagioclase. 
Biotite is replaced by kaolinite and siderite, 
but areas of slightly altered biotite are present. 
The other mafic minerals are replaced by iron- 
rich carbonate, which is also common in vein- 


The montmorillonite-kaolinite outer zone is 
1-2 inches wide and is softer and grayer than 
the inner zone. All plagioclase crystals have 
been altered into montmorillonite and kaolinite 
mixed in similar amounts. Some biotite plates 
have been replaced by a mixture of phlogopite, 
siderite, and quartz. As in the inner zone, 
potassium feldspar, quartz, and apatite are 
fresh, but magnetite is oxidized. Hornblende 
and augite are less altered in the outer zone 
The zone limits are sharp. 

Wright (1954) described in detail the wall- 
rock alteration related to a uranium-bearing 
vein in the 1040-foot level of the same mine, 
540 feet below the Idaho tunnel level. His find- 
ings coincide only partially with the writer's 
Outward from the vein, he finds on one side 
three zones: (1) sericite-kaolinite zone, bleached, 
hard; (2) montmorillonite-kaolinite zone, gray 
green, crumbly; (3) sericite (kaolinite) zone, 
white, compact. On the other side there is only 
a sericite zone analogous to (1) and (3). In 
the kaolinite pseudomorphs of the writer's 
inner zone hydromica or sericite reach in places 
20 per cent, thus approaching the composition 
of Wright’s inner zone. Zone 3 is not represented 
in the veins here described. Rock fragments 
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found in the mine dump show, however, sulfide 
yeinlets surrounded by sericitized wall rock. 
Apparently, two separate phases of alteration 
are partially superimposed in this deposit. 

In the Idaho tunnel, the granite-injected 
schists of the Idaho formation show a wide- 
spread argillic alteration similar to that found 
in the Georgetown area (pattern II). 

PATTERN V: Pattern V may be considered as 
an underdeveloped pattern VI, from which it 
departs only by the absence of the orthoclase 
zone. It has been found only in the Corning 
tunnel near Gold Hill and the Lily tunnel near 
Nederland. Since only one vein was accessible 
in each mine, pattern V may be typical of the 
whole mine or of only the vein sampled. It may 
also be present in deposits where pattern VI 
is predominant. Pattern V corresponds to the 
zoning described by Lovering (1941) and Lover- 
ing and Tweto (1953) from the same area, and 
by Sales and Meyer (1948) from Butte, Mon- 
tana. Its megascopic and microscopic charac- 
teristics are identical, with the above exception, 
to those of pattern VI. 

The Corning tunnel, located half a mile 
northeast of Gold Hill, cuts the Iowa vein 200 
feet from the portal. The wall rock is a medium- 
grained biotite granodiorite (Boulder Creek 
granite). The vein belongs to the pyritic gold 
type and consists here of fine-grained quartz 
stringers, each of which shows its own altera- 
tion zoning where separation is sufficient. Thin 
fissures filled by dolomite show the same 
pattern. The width of the zones is roughly pro- 
portional to that of the veins. Ferruginous 
bands parallel the veins at a distance within 
the inner zone; the inner part of this zone has 
been bleached by probably supergene solutions 
that dissolved the iron-rich carbonates and 
precipitated the iron hydroxides near the outer 
boundary of the zone, generally leaving a light- 
green-gray margin of the unbleached hydro- 
mica zone. The kaolinite zone is well developed 
and characterized by white clayey pseudo- 
morphs that contrast with the relatively un- 
altered biotite. 

In the Lily tunnel (Lovering and Tweto, 
1953) there is access only to a thin quartz ferbe- 
tite vein about 50 feet from the portal. The 
zoning, that is obscured by ferruginous staining 
and pre-Laramide sericitization, is asymmetri- 
cal; the zones are wider and better defined in 
the hanging wall and contain in the outer part 
of the kaolinite zone abundant illite-mont- 
morillonite regular mixed-layer clay, which is 
absent in the foot wall. 

PATTERN VI: Pattern VI represents the most 


complex zoning in the mineral belt and is found 
in the Nederland tungsten and Gold Hill 
pyritic gold districts, Boulder County. More 
than 12 mines! were studied; the alteration 
zoning is consistent, and a description of each 
therefore will not be given here. The only dif- 
ference in the alteration of tungsten and pyritic 
gold veins is a ferruginous staining (Pl. 1, fig. 
1) that is either absent or weak in the tungsten 
veins, in accordance with the relative absence 
of pyrite. The mode of occurrence of the dif- 
ferent secondary minerals can be taken as repre- 
sentative of the entire mineral belt and will 
therefore be described in some detail. 

Most of the veins cut Boulder Creek granite 
(Lovering and Tweto, 1953), a highly migma- 
titic rock ranging from medium-grained granite 
to biotite-hornblende tonalite gneiss. The most 
common facies is a medium-grained gneissoid 
biotite-adamellite. Masses of fine-grained and 
aplitic facies are abundant (e.g., Hurricane 
Hill), as are pegmatite dikes. The plagioclase 
(calcic oligoclase) is generally replaced by 10-20 
per cent of sericite of pre-Laramide age. 

There are three clearly visible alteration zones 
(Pl. 1, fig. 1). The innermost zone is texture- 
less, grayish white with yellowish jarosite stain- 
ing, and moderately friable; biotite is discolored. 
The middle zone, less conspicuous and narrower 
than the other two, is greenish gray and of 
normal hardness; it grades into the inner zone, 
whereas its outer boundary is marked by a sharp 
decrease in hardness. The third zone is friable 
and characterized by white argillic pseudo- 
morphs after plagioclase, contrasting with the 
moderately fresh biotite; it grades into the 
fresh rock. This zone corresponds to Lovering’s 
argillic zone, and the middle zone is roughly 
what he called “hard casing” (Lovering, 1941). 
To each zone corresponds a different mineral 
assemblage: inner zone, hydromica-phlogopite- 
iron carbonate; middle zone, orthoclase- 
kaolinite-hydromica; outer zone, kaolinite- 
montmorillonite. In many veins a closer study 
reveals a fourth zone in which montmorillonite, 





1Sample localities include: tungsten: Chance 
mine, Nederland, first level; Clark tunnel, Neder- 
land, Barker 1 vein; Clyde tunnel, Nederland, 
Footwall vein; Cold Spring mine, Nederland, 
fifth and seventh level; Crackerjack tunnel, Mil- 
lionaire Gulch; Eureka tunnel, Middle Boulder 
Creek; Hoosier-Hummer mine, 225-foot level; 
Prospect tunnel, Nederland. Pyritic gold: Emanci- 
pation (Gardner) tunnel, Gold Hill; Ingram 
tunnel, 120-foot level, Gold Hill; Logan and 
Grand Republic mines, Gold Hill, tunnel level; 
Poorman mine, Gold Hill, tunnel level; Snowbound 
tunnel, Gold Hill. 
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lying between the kaolinitic zone and the fresh 
wall rock, is the only clay mineral. A fifth, 
outermost zone may be distinguished in which 
the only alteration is a partial replacement of 
the mafic minerals by phlogepite and carbonate. 


Ka Or 





orthoclase, and barite. Cavities filled with opal 
and kaolinite are common. The “dark horn” 
(Lovering, 1941) contains abundant interstitial] 
ferberite. 

The hydromica zone coincides partially with 








FiGURE 9.—IDEALIZED Cross SECTION OF A VEIN SHOWING ALTERATION PATTERN VI 


See Figure 8 for definition of symbols. Ph, phlogopite-siderite alteration of biotite. Iron-staining 
bands between hydromica and orthoclase zones. 


Figure 9 shows an idealized section of a vein 
and its alteration zoning. The curves indicate 
the relative abundance ef the main secondary 
minerals. The zoning was studied, wherever 
possible, on both sides of the veins and was 
found to be essentially symmetrical in most 
veins. In the following description, each zone 
is designated by its critical secondary mineral. 
Thus, the successive zones, outward from the 
vein, are: hydromica zone, orthoclase zone, 
kaolinite zone, montmorillonite zone, and fringe 
zone. Figure 10 shows the type and degree of 
replacement for each of the main primary con- 
stituents of the wall rock. Alteration zones are 
generally described in the veinward sequence, 
a procedure that is based on the belief that the 
zones have developed successively from the 
fresh rock toward the vein. The writer believes 
that this assumption is unjustified because the 
zones may develop simultaneously, each under 
a particular physicochemical environment. This 
is thought to be the case of the zone patterns 
described in this paper. 

The veins are formed by microcrystalline 
quartz (“horn quartz”, Lovering, 1941); in 
places they have drusy cavities lined by quartz, 
kaolinite, and, in tungsten veins, ferberite. 
Pyrite is fine-grained and disseminated. The 
microscope shows interstitial carbonate, wood- 
houseite, kaolinite, and hydromica, scanty 


shear zones in many veins; cataclasis, in addi- 
tion to the discoloration of biotite, causes the 
texturelessness of this zone (PI. 1, fig. 1). In 
a few veins, completely altered biotite does not 
occupy the entire width of the zone (Pl. 1, 
fig. 2). Yellowish jarosite staining is generally 
present. In the pyritic veins, staining by iron 
(and manganese) hydroxides hides the transi- 
tion with the orthoclase zone. This staining, 
very conspicuous in most places, is banded and 
forms a veneer on old exposed surfaces; on 
fresh surfaces the hydroxides appear in the 
intergranular interstices of the rock. Gypsum 
also coats the wall rock in this zone. 
Mineralogically, the inner zone is charac- 
terized by the replacement of plagioclase by 
hydromica (Pl. 2, fig. 1), and—generally—the 
discoloration of biotite to pale-green phlogopite 
(Pl. 2, fig. 2). Iron-rich carbonate (both siderite 
and ankerite have been detected) appears as 
lenses in the phlogopite and also forms irregular 
aggregates and veinlets; these are scarce in the 
bleached part of the zone. The hydromica 
pseudomorphs enclose granules of iron carbon- 
ate and woodhouseite; the phosphate may 
produce a yellowish tint. If pre-Laramide seri- 
cite is present, its well-defined lamellae con- 
trast with the smaller, fluffy individual crystals 
of hydromica. Quartz and opal are present in 
the pseudomorphs in certain veins (e.g., Cold 
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Spring mine). X rays reveal small quantities 

of kaolinite intimately mixed with the potassic 

day. | | | 
The width ratio of hydromica zone to vein 


ZONES 
Or Ka 


Hydromica Kaolinite 


kaolinite and quartz take partially the pace 
of phlogopite and carbonate (Fig. 11). The 
orthoclase zone is from 2 to 6 inches wide in 
most of the veins studied by the writer. 


Mt Fringe 


Unaltered Plagioclase 


Unaltered Quartz 


Unaltered Microcline 


> Ko,Mt 








DISTANCE FROM VEIN 


Kaolinite, Nontronite 





Ficure 10.—DIAGRAM SHOWING DEGREE AND TYPE OF ALTERATION OF MAIN PRIMARY 
MINERALS IN THE PATTERN VI VEINS 
The relative abundance of the primary minerals is indicated by the width of the blocks. a, range of 
rapid expansion of the hydromica front during the initial stage; b, c, expansion of hydromica and ortho- 
clase fronts, respectively, during the main stage of alteration. 


is about 5:1. Maximum thickness was found 
in the Clyde tunnel, where the zone extends 
almost 25 feet (Fig. 16). 

The orthoclase zone retains in part the hard- 
ness of the fresh rock because of the presence 
of the feldspar in pseudomorphs after plagio- 
clase, in contrast with the exclusively argillic 
alteration of the adjacent zones. Secondary 
arthoclase preserves the crystallographic orien- 
tation of the host, whether it replaces the whole 
or only veinlike parts of it (Pl. 2, figs. 3, 4). 
The pseudomorphs contain also kaolinite (PI. 
2, fig. 5) and hydromica. Orthoclase is easily 
distinguished from the microcline by the lack 
of twinning and abundance of inclusions and 
appears in places to have replaced plagioclase 
along cracks (PI. 2, fig. 4). Hydromica follows 
cleavage and twinning planes, whereas kaolinite 
forms compact aggregates which do not reflect 
the structure of the host. Most of the biotite is 
preserved in the orthoclase zone, though its 
plates are paler and less straight than normal 
and lack their characteristic luster. Many mica 
plates show partial phlogopite carbonate altera- 
tion that appears to the naked eye as pale- 
green-gray flat lenses intercalated in the 
Cleavage planes. Near the kaolinite zone, 


The kaolinite zone has a sharp inner bound- 
ary and transitional outer boundary. The white 
argillic pseudomorphs have almost exclusively 
kaolinite, in places with some orthoclase 
stringers, near the inner boundary. In most 
of the zone, however, montmorillonite is also 
present, amounting to about 50 per cent of the 
pseudomorphs in some veins. Pseudomorphs 
with abundant montmorillonite are grayer and 
more crumbly (e.g., Logan group; Snowbound 
mine). Montmorillonite forms interstitial films 
between kaolinite individuals. In this zone 
biotite shows the type of alteration present in 
the orthoclase zone, but the quartz-kaolinite- 
montmorillonite assemblage (Fig. 11) is here 
more abundant. In the Chance and Prospect 
mines, the kaolinite in the biotite is fibrous. 
The kaolinite zone is about as thick as the hy- 
dromica zone (average, 1 foot) and roughly 4-5 
times thicker than the orthoclase zone. 

The montmorillonite zone, as defined in this 
paper, begins where kaolinite becomes negli- 
gible; in this zone montmorillonite is the only 
clay mineral replacing plagioclase. Since this 
replacement is partial, the total amount of 
montmorillonite may be, and commonly is, 
less in this zone than in most of the kaolinite 








~I 
bdo 


zone (Fig. 9). The transition between the two 
zones begins with the appearance of unreplaced 
areas in plagioclase, separated by argillic vein- 
lets of decreasing width. The veinlets show a 
kaolinite-rich argillic mixture in the middle 
and a wall zone of montmorillonite. The dis- 





FIGURE 11.—ALTERATION OF BIOTITE IN 
THE KAOLINITE ZONE 
Ka, kaolinite, with intercalated montmorillonite 
films. Q, quartz; S, siderite; J, jarosite. 


appearance of the kaolinitic middle part marks 
the beginning of the montmorillonite zone. 
The alteration in this zone is shown in Figure 6 
of Plate 2. The only megascopic sign of this 
alteration is a dull luster in the cleavage sur- 
faces of the plagioclase. Biotite lacks the quartz- 
kaolinite-montmorillonite replacement typical 
of the kaolinite zone; otherwise it shows the 
dull luster and interleaved phlogopite carbonate 
layers that characterize the latter zone. Horn- 
blende (Cold Spring mine) is partially or totally 
altered into an intimate mixture of iron carbon- 
ate, kaolinite, and nontronite. 

In about half of the veins studied (Ingram, 


Logan, Poorman, Hoosier, Chance, and 
Crackerjack) montmorillonite and kaolinite 


disappear at nearly the same distance from the 
vein; the montmorillonite zone is therefore 
lacking here. This zone is everywhere much 
narrower than the kaolinite zone. 

The fringe zone lies beyond the last trace of 
montmorillonite in plagioclase, where the only 
hydrothermal alteration is the partial phlogo- 
pitization of biotite and (at Cold Spring) the 
iron carbonate kaolinite-nontronite replace- 
ment of hornblende. This zone is equivalent to 
the propylitic zone of other areas (Sales and 
Meyer, 1948). 

SUBPATTERNS Via AND Vib: In variations of 
pattern VI montmorillonite enters as a regu- 
larly interstratified constituent of hydromica 
either in the outer part or the whole of the 
innermost zone (Fig. 2). The first, pattern VIa, 
was found in the Crackerjack (Fig. 12) and 
Prospect tunnels; pattern VIb is represented 
by the Chance mine and the Clark tunnel. 

The presence of the mixed-layer clay mineral 
cannot be established without laboratory exam- 


F. GONZALEZ BONORINO—HYDROTHERMAL ALTERATION, FRONT RANGE 


ination, but in places it can be suspected from 
some of the physical properties of the clay 
Rocks containing the mixed-layer clay dis. 
aggregate more easily than those containing 
hydromica; where the two clay minerals occy; 








Hydromica (mi | Kaolinite Fringe 





Ficture 12.—Cross SECTION OF THE ALTERATION 
ZONES AT CRACKERJACK TUNNEL 
Subpattern VIa. See Figure 8 for definition of 


symbols. Ph, phlogopite-siderite alteration 
biotite; ML, mixed-layer 1:1  illite-montmoril- 
lonite. 


in adjacent zones, as in the Crackerjack tunnel, 
this difference may show up very clearly. In 
the Crackerjack mine, in a little more than ? 
inches the alteration of the plagioclase passes 
successively from hydromica to mixed-layer 
clay, orthoclase, and kaolinite. 

Some veins with pattern VI zoning may show 
an incipient development of subpattern VIa, 
as the randomly interlayered montmorillonite 
in the hydromica becomes more abundant 
toward the outer part of the hydromica zone. 
In the Cold Spring mine (fifth level) for in- 
stance, the content of montmorillonite in the 
hydromica increases from about 13 per cent 
in the main part of the zone to about 20 per 
cent in the outer part. This is accompanied by 
a decrease in the refraction indices from about 
1.582 to about 1.570 for gamma. 


CHEMICAL CHANGES 


Table 4 shows the chemical composition of 
five samples from the vein and four successive 
zones of the Chance mine. This mine was chosen 
because the fine grain ahd homogeneous com- 
position of the wall rock ensures more nearly 
representative samples than in most other 
mines. The variation curves (Fig. 13) have been 
calculated correcting for the bulk density of 
the rock as given in Table 4. It may be pointed 
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SiO, in the hydromica zone are due to replace- 
ment of plagioclase by the clay mineral, which 
in this mine is a mixed-layer species with 50-50 
illite-montmorillonite. The gain of MgO must 


Kaolinite Fresh 
zone rock 
| | 





A\2O3 














with atmosphere. Expansion occurs at the face 
of the workings by relief of pressure and hy- 
dration. Bulk-density values close to 2.0, 2° 
reported in some papers, must be abnormal and 
probably do not represent the density of the 
rock immediately after alteration. 

Bulk chemical changes in the successive zones 
are moderate and probably indicate a limited 
transfer of matter. Only plagioclase and, to a 
lesser degree, the mafic minerals were affected 
by the alteration. In Figure 14 the gains and 
losses have been calculated with respect to the 
iresh rock. Solid blocks represent gains, ruled 
blocks represent losses; the height of the blocks 
is proportional to the percentage of the con- 
stituent in the fresh rock divided by its per- 
centage in the altered rock, plotted on a 
logarithmic scale. 

The gains of K,O and Al,O3 and the loss of 





FicurE 13.—CHEMICAL CHANGES IN THE ALTERATION ZONES, CHANCE MINE 


be attributed to both phlogopitization and 
formation of the clay mineral. The orthoclase 
zone shows a gain of KO, which has been 
limited by the presence of kaolinite in the 
pseudomorphs after plagioclase. The kaolinite 
zone lost SiO, and gained Al,O3; and K,O; the 
alkali may be accounted for by the presence of 
fine veinlets of orthoclase in the clay pseudo- 
morphs. With the alteration of plagioclase 
Na,O and CaO have been strongly leached from 
all zones. The proportion of H,O follows closely 
the clay content of the zones. 

The broken lines in Figure 14 represent quali- 
tatively the variations in the oxide concentra- 
tion in the hydrothermal liquids as they mi- 
grated through the various zones. A gain in the 
rock corresponds to a loss in the solution and 
vice versa. The slope of the curve is qualita- 
tively proportional to the magnitude of the 
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gain or loss in each zone, assuming that the 
solutions travel away from the vein. 
The chemical changes in the Chance mine 


TABLE 4+.—CHEMICAL ANALYSES OF THE ALTERED 
WaALL-Rock ZonEs, CHANCE MINE, NEDERLAND, 
BouLpER County, COLORADO 
Doris Thaemlitz, analyst for FR155a, FR155b, 
FR158, FR159; Eileen H. Oslund, analyst for 
FR157; Rock Analysis Laboratory, University of 
Minnesota. 








Vein | mica | clase | linite viene 

ee igh ae : rock 
‘F R15Sa, zone | zone | ZOne |ipisg 
[FR155b, FR157 | FR158 ‘ 
SiO. | 88.97 | 74.96 | 80.85 | 78.29 | 79.55 
Al.O; $.23 | 12.38 8.02 9.09; 8.45 
Fe.0; 0.28} 0.68! 2.07| 1.50] 2.48 
FeO 035; 1.26, 1.89) 2.94) 2.57 
MgO 0.43 1.45 1.02| 1.29! 1.25 
CaO 0.09 | 0.20; 0.11 0.15 1.21 
Na:O 0.13 | 0.14) 0.07 0.06 1.36 
K,O0 4.29 | 5.06; 3.13 | 2:4 1.37 
H,0+ 1.07 ee 1.35 2.17 | 0.69 
H,O— 0.82 1.18 | 0.45; 0.57 | 0.12 
TiO, 0.22 | 0.28; 0.49; 0.56! 0.62 
PO; 0.08 | 0.03) 0.03 | 0.03 | 0.02 
MnO 0.01 0.03 | 0.03 | 0.05 | 0.03 
99.47 99.38 | 99.76*| 99.12 | 99.72 
Bulk 2.48 2.45 2.49 | 2.48 | 2.62 

density 


are very close to those found by Lovering in 
the Cold Spring vein (Lovering, 1941, p. 242, 
Table I). Unfortunately, analyses of the ortho- 
clase zone of the latter vein are not available. 
The hydromica zones show a gain in MgO in 
the Chance mine and a loss in the Cold Spring 
mine. This may be attributed to the composi- 
tion of the clay mineral, that at Cold Spring 
has only little interlayered montmorillonite. 


ORIGIN OF HyDROTHERMAL ALTERATION 
General Statement 


The following discussion refers—unless other- 
wise stated—to pattern VI and its variations, 
although it can in general be applied to the 
other patterns. 


Alteration and Vein Formation 


Study of the veins in the Front Range 
mineral belt indicates that wall-rock alteration 
and vein filling took place in two separate 
phases. This conclusion is based on the fol. 
lowing: 

(1) Fissures containing no vein filling byt 
showing a well-developed alteration zoning 
are common. These fissures were reached by 
the thin, gaseous altering fluids but not by the 
more viscous, less strongly propelled veip. 
filling solutions. 

(2) Contemporaneity of the phases would 
mean that the altering fluids escaped into the 
wall rocks directly from the silica-rich, vein. 
forming liquids that filled the veins. Under 
such conditions the escaping phase would be 
expected to have carried substantial amounts 
of silica into the wall rocks. In general, no such 
silicification is found. If fractionation was in- 
complete, and the vein liquids still carried 
volatile constituents, some degree of silicifica- 
tion might have taken place (e.g., at Brecken- 
ridge). 

(3) Vein filling tends to seal off the fissures. 
Therefore, the altering phase must have ar- 
rived before the vein-forming liquids. 

The contrasting physicochemical behavior 
of the two phases becomes apparent if the per- 
vasive character of the wall-rock alteration is 
compared with the restricted, sharply bounded 
vein filling. 

Lovering (1941) discussed in detail the forma- 
tion of the tungsten and related veins of Boulder 
County, Colorado, and concluded that the 
argillic (kaolinite-montmorillonite) alteration 
preceded ore deposition, but that sericite (the 
writer’s hydromica) and adularia were the 
product of a post-ore, alkaline phase (Lovering, 
1941, p. 267, Fig. 17). He finds that the veins 
are genetically related to biotite latite dikes. 
Hydrothermal solutions would have escaped 
from the latite magmas as an acid, gaseous 
phase which later condensed and became alka- 
line by reaction with the wall rocks (Lovering, 
1941, p. 274-279). 

The present writer believes that the above- 
mentioned secondary minerals were formed in 
the same stage, before the formation of the ore 
veins. 


Hydrothermal Differentiation 


On separation of a hydrothermal phase in- 
cluding compounds of wide range of solubility 
and volatility, further differentiation may have 
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taken place with separation of a gaseous phase 
out of the liquid phase. The gaseous phase, 
containing a large part of the water and most 
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fluorine would escape as SiFs, the boiling point 
of which is 178.3°K at 1 atm. (Kelley, 1935) 
rather than as WF,, the boiling point of which 
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Figure 14.—D1acRAm OF GAINS AND LossEs, CHANCE MINE 


HM, hydromica zone; Or, orthoclase zone; Ka, kaolinite zone. Dashed curves represent 
qualitative concentration changes in the fluids. 


of the acid radicals including CO, and ap- 
preciable alkali chlorides, would have advanced 
ahead of the liquid phase. The liquid phase con- 
tained most of the silica and carbonate, as well 
as the ore metals. 

The relatively high volatility of the ore-metal 
halides seems to indicate that the metals be- 
longed to the gaseous phase. In the veins here 
studied, however, it is evident that most were 
deposited with the late silica-rich, vein-form- 
ing liquids. (See Lovering, 1941, p. 247-279.) 
Ferberite, for instance, is contemporaneous 
with at least much of the “horn” quartz, par- 
ticularly with the coarser, subhedral variety. 
Drusy ferberite is common. Thus, the ore metals 
did not in any considerable extent escape with 
the vapor phase upon fractionation. The reason 
ior this may be that in a silica-rich solution 


is 290.4°K at 1 atm. (Kelley, 1935). The same 
would hold for the chlorine compounds. 

The degree of fractionation of the ‘‘mother 
liquor” would depend upon the pressure 
gradient and the distance traveled. Lovering 
(1941, p. 277) states that “relatively open 
channels were available and, therefore, the 
pressures driving the solutions upward were 
not abnormally large’. These conditions would 
have favored the separation of the two phases. 

On the way upward the gaseous phase would 
have diffused into the wall rocks and caused 
alteration. Alteration could have lasted only 
until arrival of the vein-forming liquid which 
sealed off the fissures. If the liquids arrived at a 
given point at an early stage of fractionation 
the lateration halo would be narrow. Veins 
formed under these conditions would show no 
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constant relation between thickness of the vein 
and that of the alteration halo. The fact that 
such a relation exists in most veins of pattern 
VI indicates that fractionation has taken place, 
since the volume of the gaseous phase must be 
roughly proportional to that of the vein-filling 
column. A time gap probably separated the 
wall-rock alteration from the vein formation; 
because of differences in mobility, this gap 
would widen upward. 

Probably no important change took place 
in the hydrothermal fluids as they traveled 
along the fissures. The portion not penetrating 
into the walls moved along too fast for any 
considerable reaction with the wall rocks. Many 
geologists have postulated that acid solutions 
ascending along fissures would gradually turn 
neutral or alkaline by reaction with the walls, 
on the assumption that the dissolved matter 
diffuses inward into the fissures. Assuming that 
the alteration progresses outward from the con- 
duits it seems more logical to expect a pH de- 
crease caused by gradual insulation of the wall 
rock by the alteration sheath. According to the 
hypothesis presented in this paper, no sig- 
nificant change of pH has taken place in the 
ascending solutions during the development of 
zoning in the Front Range, because the dis- 
solved matter was carried away from the fissures 
through the successive zones and into the fresh 
wall rock by the bodily diffusing fluids. The 
writer believes that alteration progressed 
simultaneously in all zones, and that no totally 
altered, insulating sheath was therefore formed 
until the end of the process. 


Mechanism of Zone Development 


Sequential versus simultaneous origin of 
zones.—The zoning of the hydrothermal altera- 
tion has been attributed to one or more of the 
following mechanisms: 

(1) Waves of hydrothermal solutions in- 
vading the rocks in separate phases or stages. 
Most workers using this concept (Lovering, 
1949; Kerr, 1950; Kerr ef a/., 1950; Peterson 
et al., 1946; Schwartz, 1947; 1953; Stringham, 
1953) imply definite time intervals between the 
waves or at least between their peaks. Chemical 
and/or physical properties of the fluids would 
vary from wave to wave. The stage concept 
has been applied mainly to porphyry copper 
deposits, where the spatial relationships of the 
zones are somewhat obscure. 

(2) A single wave, the composition of which 
changed as the solution circulated along the 
channels (Lovering, 1941; 1950). Except for the 


continuity of the fluid invasion, this proces 
may be theoretically considered similar to the 
one above. In both processes the zones ar 
formed in a time sequence. 

(3) Zones formed simultaneously by physical 
and/or chemical fronts as the solutions pj. 
grated away from the conduits into the wal 
rocks (Sales and Meyer, 1948). 

Mechanism (2), as generally applied to veins 
(Lovering, 1941; 1950), provides for the forma. 
tion of each zone at the expense of the inne 
part of the preceding one, so that the zon 
closest to the vein is the youngest. 

The case for mechanism (3), simultaneoys 
formation of the zones, has been convincingly 
stated by Sales and Meyer (1948, p. 32) with 
regard to the veins at Butte, Montana. Their 
chief argument is that the same arrangement 
of zones is found throughout the district 
“around every ore-bearing fracture, regardless 
of size, attitude or relative age. . . .” A separate 
origin of the zones would require ‘‘that the in- 
vading zone of sericitization never advanced 
completely through the argillized zone to attack 
fresh rock directly.” A zone pattern produced 
by successive, encroaching waves of alteration 
must depart considerably from the regular 
arrangement invariably found at Butte. 

Zoning at Butte is almost identical to the 
writer’s pattern V, and Sales and Meyer’ 
reasoning can just as properly be applied to the 
veins studied by the present writer. The strong- 
est argument in favor of a simultaneous origin 
seems to be the striking persistence of the 
relative thickness of the zones and the lack 
of transgressive boundaries. There are clear 
indications in the veins of the Front Range 
mineral belt, that the formation of each zone 
was the direct consequence of the processes 
taking place during the formation of the next 
innermost zone. The changes in composition 
must have taken place as the fluids migrated 
through the wall rocks away from the veins, 
rather than during their ascent along the 
fissures. 

The theory of simultaneous, parallel fronts 
seems to be the most logical explanation for the 
structural features of the multiple zoning in 
the Front Range mineral belt. The behavior of 
the fronts during the alteration presents several 
possibilities. Sales and Meyer (1948, p. 33) be- 
lieve that at Butte “each zone migrated away 
from the fissure, that is,: it grew at its outer edge 
and simultaneously receded at its inner edge 
because of encroachment by the next innermost 
zone.” This means that when the sericite front 
reached a given plagioclase crystal the latter 
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had already been replaced by kaolinite. As far 
4s the Boulder County deposits are concerned, 
the evidence is contrary to the idea of front 
migration. 

Figures 1-5 of Plate 2 show the typical 
textural relationships found in the hydromica 
and the orthoclase zones. Hydromica seems to 
have replaced only plagioclase; the aggregates 
generally reflect the polysynthetic twinning. 
The inherited optical orientation of orthoclase 
indicates the same for this mineral; in one 
instance, relict albite twinning was observed 
in the secondary feldspar. No textural indica- 
tion of replacement of kaolinite by either of the 
other two secondary minerals has been found. 
Kaolinite forms massive aggregates which do 
not reflect the structure of the host and con- 
trast with the more ragged, structurally con- 
trolled hydromica or sericite, thus, giving in 
places a false impression of mutual replacement 
(See Sales and Meyer, 1948, Fig. 4.) The inter- 
mingled kaolinite and montmorillonite show 
strong indications of having crystallized to- 
gether. These facts oppose the idea of displace- 
ment of zone fronts as expressed by Sales and 
Meyer and indicate that the zones developed 
simultaneously and the alteration fronts became 
nearly stabilized early in the process. As the 
origin of the alteration is discussed in more de- 
tail, it will be evident, however, that a limited 
expansion of the fronts did take place. 

Nature of alteration fronts——The formation 
of the zone pattern VI can be described in terms 
of six fronts. Four represent the alteration of 
plagioclase; the other two correspond to the 
alteration of biotite: 


Fronts for plagioclase: 

(1) Hydromica front: outer boundary of the zone 
of stability of hydromica; inner boundary of the 
orthoclase stability zone 

(2) Orthoclase front: outer boundary of the 
stability zone of orthoclase; inner boundary of the 
kaolinite zone 

(3) Kaolinite front: outer boundary of the 
stability zone of kaolinite 

(4) Montmorillonite front: outer boundary of 
the stability zone of montmorillonite (and of the 
instability zone of plagioclase), in places coincides 
with the kaolinite front 
Fronts for biotite: 

(1) Fringe front: outer boundary of the partial 
alteration of biotite (to phlogopite-siderite) 

(2) Phlogopite front: outer limit of the complete 
decomposition of biotite. This is much sharper than 
the fringe front and generally coincides approxi- 
mately with the hydromica front 


~I 
~s 


There is also an instability front for horn- 
blende, roughly coincident with that of biotite, 
or fringe front. 

Fronts may be determined by one or more 
of the following factors: (1) temperature, (2) 
pH, and (3) chemical composition of the fluids. 
The temperature gradient in the wall rocks has 
been investigated by Lovering (1950) who con- 
cluded that, for a period of 10,000 years, the 
alteration takes place in ‘an almost isothermal 
environment” (p. 251). Although nothing is 
known of the duration of the process in this 
instance, the gradient probably could not have 
been steep enough to account for the formation 
of as many as four or five fronts within a 
distance of a few feet or even inches, even 
though the differences in the temperature 
fields of the minerals involved were larger than 
they seem to be, as well as in right order for 
the zone sequence—a doubtful condition. In 
a region of steep geothermal gradient, as in 
the mineral belt during Laramide times, the 
difference in temperature between hydro- 
thermal fluids and wall rock must have been 
slight . The writer concludes, therefore, that 
temperature was a subordinate factor in the 
formation of the alteration zoning. 

Changes in pH cannot explain the zoning. 
If other factors are constant, the passage from 
hydromica to orthoclase to kaolinite requires 
first an increase of alkalinity, followed by a 
sharp turn into acidity. Such changes would be 
difficult to explain. 

Chemical gradients, determined by changes 
in composition of the hydrothermal fluids dur- 
ing their course through the wall rocks, remain 
the only major factor in the formation of zones. 
Temperature may have affected the position 
of the fronts or the width of the alteration halo; 
pH, a significant factor, must be considered in 
connection with the chemical gradient. 

Development of alteration fronts ——The fluids 
permeating the wall rocks contained, as far as 
can be ascertained from the alteration products, 
H,0, CO2, S=, PO=, Cl-, F-, K+, Mgt". Their 
temperature was well above the lower limit of 
stability of sericite in acid solutions, which 
Gruner’s (1944) experiments put at about 
350°C., and was probably between 400° and 
500°C. within the conduits. The fluids were 
probably supercritical and moderately acid. 
In this environment, plagioclase, biotite, and 
hornblende were unstable; other minerals, such 
as quartz, potassium feldspar, and most ac- 
cessories, were either stable or metastable. As 
the concentration of the different chemical 
constituents changed by reaction with the wall 
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rock, chemical fronts were established, and 
the alteration zones began to form. 

Figure 15 shows the stability fields of the 
main hydrothermal minerals involved, plotted 


In the system represented by Figure 15, the 
hydromica field is enlarged by an increase oj 
Si concentration with respect to Al. The fie 
boundaries of kaolinite with respect to both 





Orthoclase Fi 





ocid—'—neutral 





vi 


TO 





————> distance from vein 


nas 
~ 
a 











G Cee 





No, Ca 


Figure 15.—Stasitity Fietp oF HypDROTHERMAL MINERALS IN THE WALL-Rock 
ENVIRONMENT DuRING ALTERATION 
Short-dashed lines represent courses of differentiation of the fluids migrating through the wall 
rocks. Numbers IV, V, and VI indicate patterns. 


against the Si/Al ratio as the ordinate and the 
change of Kt and Nat + Ca** concentrations 
as the abcissa. The latter also corresponds to an 
increasing distance from the conduit; therefore, 
the course of differentiation of the fluids as 
they migrated away from the conduits can be 
followed through the successive stability fields. 
Field boundaries and variation curves have 
purely qualitative significance and have been 
traced according to the following premises: 

According to Gruner (1944), hydromica (or 
sericite) is formed out of feldspar by potassium- 
bearing solutions in acid environment above 
350°C. Below this temperature kaolinite is the 
stable phase, except in neutral or alkaline en- 
vironment where, as shown by Noll (1936), 
either hydromica or potassium feldspar will 
form. Noll also demonstrated that an increase 
in the Si/Al ratio favors the formation of 
potassium feldspar instead of hydromica from 
solutions within the system (Ca,Mg)O— 
(K, Na)»,O—AI,0;—SiO.—H.0. Although kaoli- 
nite forms preferably in acid medium, Noll’s 
investigations show that with very low K con- 
centration it will also form from silica-alumina- 
water solutions under mildly alkaline condi- 
tions, in which circumstances it is accompanied 
by montmorillonite. 


orthoclase and hydromica are determined 
chiefly by the K content, but an increase in Al 
enlarges the field of kaolinite—more with 
respect to orthoclase than with respect to 
hydromica. Kaolinite and montmorillonite 
fields overlap along a wide area, for which the 
boundaries depend principally upon the pH 
(or the Na + Ca concentration) and second- 
arily upon the Si/Al ratio. The pH variation 
curve in Figure 15 corresponds to pattern VI. 

HYDROMICA ZONE: The fluids entering the 
wall rock from the fissures reacted with 
plagioclase and formed hydromica in its place; 
they also attacked biotite, in which Mg was 
substituted for Fe without destroying the 
mica. The volume-to-volume replacement of 
oligoandesine by hydromica represented an 
addition of Al,O; of about 15 per cent, and a 
subtraction of about 10 per cent of the SiO, con- 
tained in the plagioclase. Therefore, an increase 
in the ratio of Si to Al took place in the fluids 
as they moved away from the conduits into the 
wall rocks (Figs. 14, 15). 

Most of the iron released by the alteration of 
biotite was fixed as carbonate without sig- 
nificant transport, because of the diminished 
acidity. The iron and magnesium released by 
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hornblende became also fixed as carbonate, and 
Si and Ca were carried away. 

ORTHOCLASE ZONE: After traveling a certain 
distance the Si/Al ratio in the fluids decreased 
to a value below which orthoclase instead of 
hydromica was formed at the expense of plagio- 
dase. At this point, the pH had risen to about 
neutrality by incorporation of Na and Ca from 
plagioclase. This incorporation overbalances 
the loss of K gone into the hydromica. As indi- 
cated above, neutralization must have had a 
strong effect in bringing about the formation 
of orthoclase. The orthoclase ceased to form 
at a distance from the conduit where the con- 
centration of potassium descended below a 
critical value; that is, at the orthoclase front. 

KAOLINITE ZONE: As they crossed the ortho- 
clase front, the fluids began to form kaolinite 
instead of feldspar. In agreement with experi- 
mental data (see above), montmorillonite was 
formed with kaolinite in most of the kaolinite 
zone of patterns V and VI. Montmorillonite 
began to appear a little beyond the orthoclase 
front, increased as the fluids crossed the zone, 
and decreased at a rate either smaller than or 
equal to that of kaolinite. In a few veins the 
expanding clay equals or outweighs kaolinite 
in part of the zone. The amount of montmoril- 
lonite probably depended on the Na and Ca 
concentration in the fluids. Kaolinite ceased to 
form where the pH and base content became too 
high. 

MONTMORILLONITE ZONE: As Na + Ca con- 
centration increased and the pH rose, plagio- 
clase became less susceptible to the attack by 
the fluids, which affected only the walls of 
cracks in the plagioclase, and deposited mont- 
morillonite. The influence of microgradients 
is shown by the banding of many veinlets, in 
which kaolinite occupies the middle and the 
expanding clay lines the walls. The place where 
the fluids ceased to attack the plagioclase at all 
represented the montmorillonite front. 

PHLOGOPITE ZONE: The zone of intense 
phlogopite-carbonate replacement of biotite 
coincides roughly with the hydromica zone. The 
exceptions (Pl. 1, fig. 2) show that the two 
fronts were partly independent of one another. 
The front was the place where the migrating 
fluids ceased to attack biotite on a large scale. 
Beyond it the attack decreased suddenly from 
100 to 5-10 per cent of the mica; in the kaolinite 
zone, this limited attack was more radical be- 
cause it produced in parts a mixture of kaolinite, 
quartz, and montmorillonite (Fig. 11) that 
replaced the mica and destroyed its crystalline 
framework. The stability field of phlogopite 


extended to at least the outermost, or fringe, 
front, but the ability of the hydrothermal fluids 
to decompose biotite dropped abruptly at a 
place near the hydromica front. This drop was 
apparently due to neutralization of the fluids. 
In the presence of Mg-bearing fluids, the re- 
action biotite — phlogopite may take place 
much faster in acid than in neutral or alkaline 
conditions. Neutralization was also probably a 
contributing factor in the formation of ortho- 
clase and therefore in the setting of the hydro- 
mica front. Consequently, the hydromica and 
the phlogopite fronts must have been partially 
related to one another. 

Problem of front migration—Hydromica and 
kaolinite generally extend into the orthoclase 
zone and may even overlap one another slightly. 
The coexistence of the three secondary minerals 
(Pl. 2, figs. 2, 3) may be tentatively explained 
by one of the three following processes: (1) 
simultaneous crystallization, (2) successive re- 
placement of the secondary minerals by one 
another, and (3) successive crystallization with- 
out mutual replacement. 

Simultaneous formation of the three minerals 
out of the same solution seems unlikely on 
phase-rule considerations. As to process (2), it 
has been shown above that no signs of mutual 
replacement are found in the aggregates and 
that all three minerals appear to have formed by 
replacement of the same host, i.e., plagioclase. 
Process (3) seems to be the logical conclusion, 
but it requires a shift of the fronts to produce 
the overlapping. An explanation of the shift 
may be as follows: Early in the process, the 
hydromica front was established near its final 
position, between the outer boundaries of hy- 
dromica and orthoclase zones (see Fig. 10) 
as the flow across the wall rock became stabi- 
lized. After this, the gas pressure continued to 
increase at a slower rate, and the increasing 
flow pushed the fronts outward, over-riding 
the next outer zone but replacing only the 
remaining areas of plagioclase. As the hydro- 
mica front reached the former position of the 
kaolinite front, the three secondary minerals— 
kaolinite, orthoclase, and hydromica, in this 
order—would have formed side by side. 

The zone boundaries, as arbitrarily defined 
by the writer, do not wholly coincide with the 
final position of the fronts. For instance, 
whereas the outer boundaries of the montmoril- 
lonite, kaolinite, and orthoclase zones represent 
the fronts at their maximum expansion, the 
outer limit of the hydromica zone corresponds 
to the initial front. 

The stages in the evolution of the fronts can 
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be summarized as follows: (1) Initial stage. 
Irruption of gases in the fissures, followed by 
a short period of rapid increase in the rate of 
flow through the wall rock, and then rough 
stabilization. During this period, the fronts 
moved rapidly outward from the conduits to 
assume their respective positions in the wall 
rock. This sweeping displacement was too fast 
to leave any marked effect on the rocks; the 
sort of effect to be expected is exemplified by the 
kaolinite present in small amounts within the 
hydromica zone (Fig. 10), which record the 
path of the kaolinite front. 

(2) Main stage. After partial stabilization, 
the fronts continued to expand slowly, as the 
rate of flow increased to culmination. Most of 
the alteration took place at this stage. 

(3) Waning stage. Decline in the rate of 
flow may have occurred by simple exhaustion of 
the channel fluids, by decrease of permeability, 
or by irruption in the conduits of the vein- 
filling liquids. If, as believed, a time gap existed 
between the alteration and the vein formation, 
the last cause should be ruled out. During the 
decline the fluids had little, if any, effect upon 
the previously formed alteration zones, as 
every alteration mineral was either stable or 
metastable in the field of the others. 

During the main stage, the alteration of 
plagioclase took place simultaneously and at 
approximately the same rate in the three inner 
zones. Therefore, it was completed about the 
same time in all three zones. Assuming that this 
occurred before the flow began to decline, there 
would have taken place a sudden outward shift 
of the fronts to new positions into the fresh 
rock, and a new set of zones could have begun 
to form adjacent to the outermost zone. No 
such double zoning, however, was found in the 
veins studied by the writer. A reason for this 
may be a decrease in the flow rate caused by a 
drop in permeability as a result of the presence 
of an argillic zone which may have acted as a 
barrier for the fluids. 

Transfer of matter during alteration—The 
problem of the transport of chemical con- 
stituents during the alteration has been dis- 
cussed by several authors. Two opposite mecha- 
nisms may have active in a wall rock next to 
a fissure: (1) diffusion of ions through a static 
pore-filling solute, and (2) bodily movement 
of the solutions through pores and cracks, away 
from the channel. Diffusion is conceived as a 
two-way transport and responds to various 
concentration gradients built up during the 
chemical reaction between wall rock and alter- 
ing fluids. According to Lovering (1950, p. 


249-260) diffusion predominates in wall-tog 
alteration. Sales and Meyer (1948) attribut 
the alteration at Butte to the same proces, 
Other authors (Burbank, 1950; Schmitt, 1959) 
favor the bodily movement of solutions. 

If the pore solvent is static, ionic diffusion 
may accomplish the required transfer to bring 
about the zoning. However, concentration 
gradients built up by chemical exchange be. 
tween rocks and solvent would have been to 
weak for any considerable matter transport, 
because of the low solubility of most of the 
constituents involved (Si, Al). The rise in con. 
centration of, for example, silica in the kaolinite 
zone would produce gradients sloping in oppo- 
site directions—toward the vein as well as 
toward the fresh rock. Apparently even a 
modest rate of outward fluid migration should 
overcome ionic diffusion. Actually, the strongest 
concentration gradient established in the wall 
rock during the circulation of fluids in the 
channel is that of the solvent itself. 

The only advantage of ionic diffusion over 
bodily movement of fluids is that it apparently 
provides a simple mechanism for disposal of the 
excess matter. In the writer’s opinion, however, 
the problem of matter disposal is not a serious 
one in the Front Range veins. First, the total 
transfer during alteration was not large, as can 
be seen from Table 4 and Figure 14. The total 
SiO, gained by the solutions from both the 
hydromica and kaolinite zones (disregarding 
the orthoclase zone, where there was a slight 
loss) in the Chance mine amounts to 0.18 gram 
per cubic cm, or 5.4 grams per square cm 
through an integrated thickness of 30 cm, mak- 
ing due allowance for the bulk densities (Table 
4). Secondly, as the fluids traveled away from 
the fissure their silica-transporting power in- 
creased substantially with the pH and the 
concentration of sodium. When the fluids 
reached the fresh rock they contained mainly 
a sodium silicate solute which must have been 
disposed of by bodily diffusion in the wall 
rocks. 


Nature of the Replacement Mechanism 


As far as can be ascertained from microscopic 
observations, the replacement of the plagio- 
clase by the various secondary minerals in the 
veins studied by the writer was generally a 
volume-to-volume process. This generalization 
is wholly valid only for plagioclase, because the 
replacement of biotite by the phlogopite-iron 
carbonate assemblage involved some volume 
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increase, indicated by the growth of the carbon- 
ate between cleavage planes. 

In a general sense the mineralogical nature 
of the alteration depends largely on the com- 
position of the host rock, but some qualifica- 
tions are necessary for the veins discussed in 
this paper. The different alteration of plagio- 
case and biotite in most zones expresses the 
influence oi the host. The difference is, however, 
not only chemical but also structural; hydro- 
mica represents a complete rearrangement and 
substitution of plagioclase; phlogopite is an 
exchange of one type of cation for another 
within the same crystalline framework. Where 
complete reworking of the biotite has taken 
place, as it has partly in the kaolinite zone, 
biotite and plagioclase were replaced by nearly 
the same minerals. Direct precipitation of hy- 
dromica and kaolinite in open fissures also 
occurred in places, showing that‘the solutions 
contained all the necessary ingredients. 

The writer believes that the composition 
of the host was essential for localizing the al- 
teration, but most of its effect on the type of 
alteration was indirect. Secondary minerals 
formed where solution of a primary constituent 
had occurred, mainly because free space was 
there available and because crystallization of 
the secondary mineral was triggered and sus- 
tained by the loss of energy spent in the attack. 
The nature of the mineral was determined pri- 
marily by the over-all composition of the fluids, 
secondarily by the composition of the replaced 
mineral. The dissolution of the host affected the 
chemical nature of the fluids and determined 
the formation of fronts as described above. The 
writer believes that the above conclusion can 
be applied to other areas. The basic idea is that 
the altering solutions are more important in 
determining the nature of the secondary product 
than the composition of the host. Sericite re- 
placing both plagioclase and potassium feldspar 
in the inner zone at Butte (Sales and Meyer, 
1948) is an example of the limitations in the in- 
fluence of the host mineral. Where the com- 
position of the replaced mineral contrasts 
strongly with that of the altering fluids, local, 
steep chemical gradients may originate crystal 
phases which—though in disequilibrium with 
the solution as a whole—can persist throughout 
the process if the rate of solution is sufficiently 
slow. It is, therefore, largely a question of de- 
gree of equilibrium between the fluids and the 
wall rock. The writer believes that the degree 
of equilibrium was relatively high in the Front 
Range mineral belt. 


Origin of Other Patterns of Alteration 


The above discussion refers to pattern VI. 
However, the processes outlined should also 
be valid for other patterns in the Front Range 
mineral belt. The discussion of the remaining 
patterns will follow approximately the order of 
increasing departure from pattern VI. 

SUBPATTERNS Via AND Vib: The chemical 
analyses in Table 4 correspond to subpattern 
VIb, in which the clay of the innermost zone 
is a mixed-layer minera]. A straight pattern 
VI would show a somewhat greater gain of 
potassium and a smaller gain of magnesium in 
the innermost zone. With hornblende present, 
magnesium might have shown a loss (Cold 
Spring). The formation of interstratified illite- 
montmorillonite clays is explained in terms of 
a lower K/Mg ratio in the fluids. This ratio 
decreased (see Fig. 15) away from the fissures 
and reached the critical value before arriving 
at the orthoclase zone; hydromica then gave 
way to the mixed-layer clay, thus forming 
pattern VIa. 

PATTERN V: Pattern V differs from pattern 
VI by the lack of orthoclase zone. The fluids 
that gave rise to pattern V probably had a 
slightly lower potassium content. Thus the 
depletion of the alkali would have come before 
the Si/Al ratio reached the adequate value for 
the formation of feldspar. Such fluids should 
have also been more acid, thus reducing the 
orthoclase field. (See Fig. 15, curve V.) In- 
creased free-acid content would also have had 
the effect of eliminating the feldspar zone. In 
the diagram this would be indicated by shrink- 
age of the orthoclase field, in which case curve 
VI also fails to enter it. The ratio of acid radi- 
cals to potassium may have increased, for 
instance, by a greater-than-ordinary distilla- 
tion of the hydrothermal solutions. 

PATTERN IV: The formation of predominant 
kaolinite in the innermost zone may be due to 
low temperature, high acidity, or low potassium 
content, or a combination of these factors—re- 
ferring to pattern VI solutions as standards. 
An evaluation of these factors will not be at- 
tempted here. Curve IV in Figure 15 shows the 
interpretation based on either one or both of the 
two last factors. 

PATTERN I: Alteration at Jamestown seems 
to have started with a pervasive, high-tempera- 
ture replacement of the country rock (greisen- 
like alteration). As this stage declined, the 
lagging less penetrative part of the hydro- 
thermal column could advance through better 
defined channels; some fissures tapped fluorine- 








82 F. GONZALEZ BONORINO—HYDROTHERMAL ALTERATION, FRONT RANGE 


rich segregations, whereas others, perhaps at a 
slightly more advanced stage of differentiation, 
served as channelways for sulfide and telluride- 
bearing liquids. Within these liquids further 
differentiation seems to have taken place; a 
light, probably gaseous phase separated and 
caused most of the replacement of plagioclase 
by hydromica, whereas the development of 
orthoclase and silicification occurred at a some- 
what later stage, upon arrival of the heavier, 
ore-bearing portion of the differentiated solu- 
tion. Apparently the wall-rock altering phase 
was rich in potassium and probably of neutral- 
to-alkaline character, as indicated by the 
absence of kaolinite. 

PATTERN II: Pattern II shows best the separa- 
tion between the altering phase and the vein- 
forming phase. The sulfide veins are not gen- 
erally accompanied by conspicuous alteration 
zones, although argillic replacement is wide- 
spread in the surrounding rocks. Some non- 
mineralized fissures show intense argillization. 
The interpretation of this pattern is handi- 
capped by the lack of information on the 
horizontal and vertical distribution of the 
alteration. The hydrothermal phase probably 
was potassium-poor. The presence of intimately 
mixed hydromica, kaolinite, and montmoril- 
lonite cannot be due to simultaneous crystal- 
lization but is probably due to successive 
replacement, first by hydromica, then by 
kaolinite-montmorillonite. These stages were 
not separate waves of alteration but probably 
represented differentiation of the fluids as they 
migrated in different directions. Apparently 
each altered crystal has been affected by solu- 
tions coming from different local sources (fis- 
sures) and at various stages of differentiation. 

PATTERN I: The problem of pattern I is what 
prevented the altering phase from following the 
normal course of differentiatign after it had 
formed the hydromica zone. The writer has 
stressed the point that wall-rock alteration 
depended upon two conditions: (1) chemical 
composition of the fluids, in so far as it deter- 
mined the minerals with which they were 
saturated; (2) the dissolving power, essentially 
related to the free-acid content. A solution may 
be saturated with a given clay mineral, but 
precipitation may not occur if the solution does 
not react with the wall rock. Precipitation 
would also fail to occur if the solution, though 
able to dissolve the primary minerals, is under- 
saturated of the ions needed to form the second- 
ary products. This solution will leach the rock 
and leave behind open spaces. In pattern I, the 
writer believes that the reacting capacity of the 


fluids became neutralized before leaving the 
hydromica field, and the saturated phase went 
on without further crystallization of secondary 
minerals. 


Silicification 


Silicification of the wall rocks was a restricted 
phenomenon in the veins of the Front Range 
mineral belt. Excluding brecciated areas, 
where silica is related to the vein-filling phase, 
and a little quartz and opal present in argillic 
pseudomorphs in some veins (e.g., Cold Spring), 
the only appreciable silicification was found at 
Breckenridge (Wellington). Silica is here 
somewhat irregularly distributed but con. 
centrates near the veins. The association of 
silicification with pattern I may perhaps be 
explained in terms of degree of differentiation 
of the hydrothermal column. If a high pro- 
portion of the volatile compounds remained in 
the siliceous, vein-forming liquid as it arrived 
at a given place, the vapor pressure would 
greatly facilitate the introduction of silica into 
the wall rock. This interpretation fits with that 
proposed above with regard to the formation of 
a single hydromica zone in pattern I. In effect, 
poor reacting ability (low acidity) in the alter- 
ing phase may be attributed also to lack of 
differentiation in the hydrothermal column. 


Alteration Pattern and Wall-rock Type 


The range of petrographic variation of the 
wall rock in the Front Range mineral belt is 
wide, but the mineralogy is monotonous: 
plagioclase, quartz, potassium feldspar, and 
mica (mainly biotite) represent more than 95 
per cent of the total composition. This is valid 
for the igneous as well as for the metamorphic 
rocks. Some _ hornblende-bearing facies of 
granite gneisses and schists are also present. 
Plagioclase is the most affected by the altera- 
tion, followed by biotite and hornblende; 
potassium feldspar and quartz are practically 
inert to hydrothermal attack. The proportion 
of plagioclase is therefore a main factor in the 
development of zoning. 

Patterns V and VI are within the area of the 
Boulder Creek batholith. Most facies in this 
body are rich in plagioclase, and the zoning 
shows exceptionally well. Where the veins cut 
zones of plagioclase-poor, schistose rocks the 
zones are weakly developed, as in the George- 
town-Argentine area. Since most of these rocks, 
however, are granite-injected schists with ap- 
preciable amounts of plagioclase, the pattern 
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liffierences cannot be attributed solely to wall- 
rock influence. 

Pattern I was studied in three different wall- 
rock types. In the Wellington mine the rock is a 
quartz latite porphyry; at Silver King it is a 
handed quartz-biotite-hornblende rock with 
abundant metasomatic tonalite, and at Oregon- 
Quaker City it is a granite. The zoning is formed 
essentially at the expense of plagioclase, and 
the three places show no fundamental dif- 
ferences. 

The effect of wall rock on the width of the 
ones Was noticed in a vein at Cold Spring, 
fith level, where the wall rock on one side is a 
biotite- and plagioclase-poor leucocratic granite 
and on the other side a biotite-rich granodiorite. 
The zones are nearly twice as wide on the 
leucocratic side as they are on the granodiorite 
side, whick is explained by the lower rate of 
chemical change of the fluids traveling through 
the plagioclase-poor rock. 


Alteration Pattern and Ore Tvpe 


The zone pattern in the mineral belt bears a 
definite relationship to the type of ore in the 
veins. The most complex zoning (patterns V, 
\I) is related to the pyritic gold and ferberite 
veins of Boulder County. Pattern ILI is found 
in the pyrite telluride gold veins of the north- 
eastern end of the belt, which are spatially and 
genetically related to fluorite deposits. Pat- 
tern IV, also associated with lead-silver-zine 
mineralization seems to be localized in a small 
area west of the tungsten district. The simplest 
zoning (patterns I and II) is generally re- 
stricted to the lead-silver-zinc veins pre- 
dominating in the southwestern part of the belt 
but is also present in scattered deposits along 
the western side of the belt up to Ward and 
Jamestown. 

As a broad generalization, then, the simple 
patterns (I-III) are associated with sulfide 
ores along the western border of the mineral 
belt, whereas the more complex patterns (V- 
VI) are related to the sub-belt situated to the 
southeast of the former where pyritic gold or 
tungsten predominates over sulfides. The latter 
includes the ores of the Central City area, 
which belong to the complex type (Tooker, 
1955; 1956). 


Alteration Pattern and Laramide Intrusions 


The petrography and sequence of the Lara- 
mide igneous intrusions of the mineral belt have 
been thoroughly studied by Lovering and 
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Goddard (1938; 1950). The intrusions pro- 
gressed in general from southeast to north- 
east; simultaneously, the magmas became 
slightly more alkalic. Correspondingly, the 
pyritic gold telluride veins as well as the 
tungsten veins are related genetically to the 
more alkalic types (latite porphyries) (Lovering 
and Goddard, 1938). The sulfide ores seem to 
be located preferably within or close to the 
stocks, whereas most of the pyritic gold, 
telluride, and tungsten veins are in the zone of 
dikes along the southeastern side of the mineral 
belt. 


Summary and Conclusions of Relation of Wall- 
rock Alteration to Laramide Intrusions 
and Ores 


Alteration patterns I and II are generally 
associated with lead-silver-zinc sulfides and 
older igneous Laramide bodies. Patterns V and 
VI are related preferably to tungsten and gold 
pyrite telluride veins and to relatively alkalic, 
Laramide dikes. Pattern IIT is related to gold 
telluride and fluorspar veins and partly alkaline 
composite stock. The alkalic (potassic) nature 
of the alteration increases northeastward and 
is at a maximum in the Jamestown district. 

The significance of these relationships is un- 
known to the writer. Perhaps the increasing 
potassium content in the alteration of the 
northeastern part of the belt is related to the 
more alkaline character of the intrusives. This 
correlation, however, is not clear. The most 
intense potassium replacement (Jamestown) is 
associated with calc-alkalic and sodic igneous 
bodies. (See Lovering and Goddard, 1950, 
p. 256-257.) The more complex patterns were 
probably related to acidic hydrothermal fluids. 
This may have some connection with younger 
age and greater distance from the igneous 
stocks, which in turn may be favorable condi- 
tions for a greater hydrothermal differentiation. 
Pattern VI seems to have been formed by solu- 
tions with high free-acid and potassium content. 

The ore type seems to have no direct con- 
nection with the wall-rock alteration. This 
conclusion is supported by the fact that in other 
districts (e.g., Butte) a zoning similar to that of 
pattern V is found in association with sulfide 
ores (Sales and Meyer, 1948). In the Front 
Range mineral belt the alteration is fairly 
homogeneous in the individual districts, al- 
though the mineralization may change con- 
siderably, as in the Boulder County tungsten 
and pyritic gold deposits. The above conclusion 
is in accordance with the idea that alteration 
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and ore formation belong to different phases of 
the process. 

No simple, direct relationship between wall- 
rock alteration and Laramide igneous bodies 
may be derived irom the available data. The 
writer believes that the character of the hydro- 
thermal alteration is only broadly related to 
the nature of the magmas and that the dif- 
ferences found within the mineral belt depend 
principally upon the evolution of the fluids 
after they leit the source areas. Structural 
conditions would, therefore, have considerable 
influence on the nature of the alteration. 

Before these and other problems related to 
the origin of the alteration in the Front Range 
can be satisfactorily answered, much more 
study is necessary. Specifically, the mapping of 
the alteration patterns should be carried to a 
much greater detail than shown in this paper. 
Pattern II should be investigated more thor- 
oughly, with the purpose of establishing the 
spatial relationship between the diffuse, argillic 
alteration and the addition directly associated 
with the veins. 


Comparison with Alteration in Other 
Areas 


Although potassium feldspar is an alteration 
product in many districts, it had never been 
found in definite zones around veins as in 
pattern VI. Among the few vein deposits where 
alteration studies have been made, the deposit 
at Butte, Montana (Sales and Meyer, 1948), 
offers the closest similarity with the deposits 
described in this paper, and in particular with 
pattern V. Minor differences are: (1) silica is an 
important constituent of the inner zone at 
Butte, whereas it is subordinate or nonexistent 
in Boulder County; (2) primary potassium 
feldspar was replaced by sericite in the same 
zone as Butte; (3) At Butte, biotite 
slightly chloritized in the montmorillonite 
zone, recrystallized into fine-grained biotite in 
the kaolinite zone, and replaced by sericite in 
the sericite zone. Fine-grained, secondary bio- 
tite is not found in patterns V or VI but in the 
area of pattern III (Jamestown). Despite these 
differences, the alteration at Butte and in 
Boulder County are similar enough to warrant 
the assertion that they must have had es- 
sentially the same origin. 

Most alteration studies in the literature refer 
to deposits with complex structural histories 
(Lovering, 1949; Kerr ef al., 1950; Kerr, 1951; 
Peterson et al., 1946; Gonzalez Bonorino, 1949; 
Schwartz, 1947; 1953; Leroy, 1954; Stringham, 
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1953). Schwartz (1956) presents an excellen; 
summary of the alteration in various areas, The 
avenues for the altering solutions are ill-define 
at best. Most authors believe that the alteration, 
took place in all instances in several separat, 
stages, each stage corresponding to a ney 
inflow of solutions. Although the number oj 
stages may be as many as seven (Stringham, 
1953) and the general picture of the alteration 
is complex, a sericite-quartz zone or area sup. 
rounded by an argillic (kaolinite, hydromica, 
montmorillonoid) zone is present almost every- 
where (Schwartz, 1956). A propylitic zone 
characterized by the alteration of the mafic 
minerals alone and equivalent to the fringe 
zone, is also found around the argillic zone. The 
close analogy of this sequence with that of the 
veins at Butte and Boulder County suggests 
common basic mechanism of alteration. 

The main grounds for the “‘stage’”’ hypothesis 
are: (1) the irregular contours of the zones or 
areas covered by each “‘stage’’, and (2) the 
replacement relationships between the minerals 
of the adjacent zones. 

(1) The front mechanism may produce zones 
that are nonparallel at the surface if the fronts 
were concentric with respect to irregular chan- 
nelways located below the present surface. 
Furthermore, solutions traveling in wide areas 
are likely to find differences in composition and 
permeability (7.e., fissures) that affect the 
extension of the zones. 

(2) Encroachment of one stage or zone by 
the adjacent one is not precluded by the front 
hypothesis, since it generally occurs during 
expansion of fronts. This encroachment does 
not necessarily entail mutual replacement of 
the secondary minerals but rather a successive 
replacement of the common host. The criteria 
generally used to recognize replacement ate 
based on textural relationships of the aggre- 
gates and are highly unreliable. Even criteria 
such as veining must be used with caution. The 
various clay minerals form aggregates with 
peculiar habits, from which no “mutuai bound- 
ary” relationships can be expected in general. 
For instance, if kaolinite and hydromica (or 
sericite) are present in a pseudomorph the 
former will make, whatever their relative age, 
compact aggregates surrounded by an I- 
regular, unevenly textured hydromica mass that 
appears to penetrate into the kaolinite. Many 
clear-cut replacement textures could be relict or 
inherited textures, resulting from the replace- 
ment of a host that has already been replaced 
by another. This happens when an alteration 
front overrides the adjacent zone. 
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excellent \lthough the writer believes that the front — to-high-temperature hydrothermal processes, 
reas. The iypothesis can also be applied to porphyry — however, the effect of supergene solutions should 
ll-define opper bodies, it is also probable that repeated — be readily distinguished. The difference lies 
alteration ipsurges of hydrothermal fluids were operative. | mainly in their CO, and H.O pressures, which 
Separate in excellent example of multistage alteration— are the main factors in their dissolving power. 
© 4 nev | hough not a porphyry copper deposit—is the In the mineral belt, cold, low-pressure vadose 
umber oj fast Tintic district, Utah (Lovering, 1949). liquids dissolved and precipitated interstitial 
tringham, However, complexity in the alteration can be matter but could not attack the primary or 
alteration brought about even by a single phase of altera- secondary silicates. 

area sur | on ii adequate structural conditions are The most noticeable supergene phenomenon 
‘dromica present. For instance, solutions radiating from in the mineral belt is the formation of fer- 
wines dispersion center may find fissures that allow — ruginous_ bands on both sides of the veins. 
<< Zone, | em to reach outer zones without further These bands are especially well developed in 
he Matic | ange in composition. A zoning will develop _ the pyritic veins of the Gold Hill district (Pl. 1, 
wile dong the fissures in apparent independence of — fig. 1). Iron and manganese hydroxides are 
one. The the diffuse alteration that surrounds it and that concentrated in fissures and intergranular 
at of the is forming more or less simultaneously. Such a spaces in the outer part of the hydromica zone; 
iggests a F tation could be interpreted as produced in exposed surfaces are covered by the hydroxides 
i . [two stages. . that migrated out of the interstices and formed 
‘pothesis In conclusion, the front theory of alteration conspicuous red-brown to bluish-black bands 
Zones of F an be applied, in greater or lesser extent, to with concretionary and mammilary structure 
(2) the almost any hydrothermally altered body, The iron staining is attributed to iron-bearing, 
minerals especially alter the tridimensional picture and acid solutions formed by oxidation of pyrite at 
the stricture of the country rock have been _ high levels. Because of its sheared texture, the 
Ce zones | ken into account. After the alteration pattern zone next to the vein is favorable for the circu- 
1€ Eronts has been analyzed with the front theory in lation of the supergene solutions, which 
ar chan. mind, the presence of more than one stage of — bleached the inner part of the zone as they dis- 
surface. | iteration can be ascertained. solved part of the siderite and ankerite out of 
de _— the altered mafic minerals. The hydroxides were 
wees precipitated as the solutions sought to migrate 
ect the laterally through the less permeable outer part 
The relative importance of hydrothermal of the zone and into the orthoclase zone. 
zone by | and supergene alteration in mineral deposits Jarosite. gypsum, and barite were also formed 
he iront | as been a matter of disagreement among out of the sulfate-bearing supergene solutions. 
during investigators. Lindgren (1915) considered that Gypsum is in places abundant and relatively 
nt does supergene waters carrying sulfuric acid from widespread (e.g., Clyde tunnel). 

nent 0! | the oxidation of pyrite were the main agent of Supergene transportation of opal and 
cessive | \aolinization, but later work (Burbank, 1932; kaolinite from hydrothermal veinlets into 
criteria } Graton and Bowditch, 1936; Lovering, 1941) younger fissures seems to have occurred in 


Supergene Alteration 


nt aft | tablished the predominantly hydrothermal many places. 
oe origin of argillization. The problem remains, 
wm ~ however, of determining to what extent, if any, Alteration as a Guide to Ore 
‘ th supergene processes affect the hydrothermal 
S Wit) | alteration products in ore deposits Kerr (1950) The use of the wall-rock alteration in 


bound: | concluded that in porphyry copper deposits the prospecting for ore has been discussed _princi- 
eneral. “argillic and sericitic minerals formed by _ pally by Lovering (1949) and Schwartz (1955). 
ica (ot hydrothermal alteration suffer little change No systematic technique has been developed, 
oh the luring later supergene processes.’’ The present although miners have long used rock alteration 
ve age; | writer agrees with this conclusion with regard as an indication of a vein. The only contribution 
an ™ | to the Front Range mineral belt. that alteration studies can give to prospecting is 
ss that The distinction between the two types of — the recognition of hidden alteration features 
Many ilteration is not always apparent. In the last — bearing a systematic relation to ore or at least 
lict ot phases of hydrothermal processes the solutions _ to the vein. If, as the writer believes, the wall- 
eplace- must behave like heat-activated vadose waters. rock alteration and the vein formation belong 
placed Intermingled vadose and juvenile solutions to separate phases of the hydrothermal process, 
TatON | vould also give rise to a series of transitional — the chances of finding such systematic relation- 
alteration phases. In environments of moderate- ships are indeed few. 
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In the Boulder County tungsten and pyritic 
gold deposits, the recognition of the complex 
zoning and, in particular, of the orthoclase 
zone might have been of some use in guiding the 
exploration work. In the widest veins, where 





FiGuRE 16.—MaAp oF ALTERATION ZONES, 
FootwaLL VEIN, CLYDE TUNNEL 


the zones may be several feet thick, the rock on 
either side of the kaolinite zone seems similarly 
fresh. Microscopic examination discloses, how- 
ever, that on one side plagioclase is preserved, 
whereas on the other side it has been replaced 
by orthoclase. The vein would be on the latter 
side. A good example of a wide orthoclase zone 
that is not outwardly different from the fresh 
rock is found in the Clyde tunnel, a mile north 
of Nederland (Fig. 16). The country rock is the 
fine-grained, biotite granodiorite or adamellite 
that characterizes Hurricane Hill and vicinity 
(Hoosier-Hummer, Clyde, and Lily mines). 
The kaolinite zone is recognized by the “salt 
and pepper” texture that results from the 
contrast between biotite and the white argillic 
pseudomorphs. The hydromica zone shows 
similar pseudomorphs, but the alteration of 
biotite to phlogopite makes the texture less 
contrasting. The orthoclase zone looks much 
like the fresh rock. The identification of this 
zone could have helped during the exploration 
work. 


GENERAL CONCLUSIONS 


The writer attempts in this paper to develop a 
“fron” theory that explains the alteration 
zoning found in the veins of the Front Range 
mineral belt. The theory is developed with re- 
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gard to the most complex zoning, pattern VI. 
and extended therefrom to the other patterns, 
The conclusions can be summarized as follows: 

(1) Wall-rock alteration and vein formation 
are two separate phases of the hydrothermal 
process. The alteration was accomplished by a 
moderately acidic, gaseous phase that separated 
by distillation from the hydrothermal solution. 
The remaining less volatile, silica-rich, alkaline. 
liquid portion arrived later, filling the fissures 
with the gangue and most of the ore minerals. 
This conclusion is based mainly on the presence 
of wall-rock alteration around fissures that 
have not been reached by the vein-filling solu- 
tions, and of veins with no alteration halo. 

(2) Zoning like that found in patterns IV to 
VI is the result of a one-stage process in which 
the fluids changed in composition as they 
migrated through the wall rocks away from 
the fissures. This is strongly supported by the 
parallelism and lack of transgression of the 
zone boundaries, as well as by the constant 
relative thickness of the zones. 

(3) The change in chemical composition of 
the migrating fluids as a result of the effects of 
the wall rock gave rise to several fronts which 
represented the stability limits of the various 
secondary minerals. In the setting of these 
fronts the hydrogen-ion concentration was a 
contributing factor whereas temperature was of 
subordinate importance. The maximum number 
of fronts (pattern VI) is six: hydromica, ortho- 
clase, kaolinite, montmorillonite, phlogopite, 
and fringe. The first four are’ defined on the 
basis of the alteration of plagioclase, the last 
two on that of the biotite and other mafic 
minerals. The hydromica and the phlogopite 
fronts are coincident in many places. 

(4) Chemical analyses indicate a significant 
increase of the Si:Al ratio in the migrating 
fluids within the innermost zone. Hydromica 
ceased to form where the above ratio reached a 
critical value (hydromica front). The field oi 
orthoclase was thus reached, with the help of a 
rise in pH. Neutralization also reduced radically 
the attacking power of the solutions with 
respect to biotite, which thereafter is only 
partially altered. The phlogopite front repre- 
sents this abrupt change in the alteration of 
biotite. 

As the fluids traveled away from the fissure 
dissolving plagioclase and forming orthoclase 
in its place, the potassium content diminished 
to a critical concentration (orthoclase front) 
beyond which the feldspar ceased to crystallize 
and gave way to kaolinite, joined later by 
montmorillonite. The introduction of mont- 
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morillonite was forced by the alkalinization of 
the fluids, mainly as a result of the incorporation 
of Na and Ca. 

Kaolinite and montmorillonite continued to 
crystallize together until the solutions became 
too alkaline for the former (kaolinite front). In 
some veins the power of the fluids to attack 
plagioclase was also lost at this point, and the 
two clay minerals ceased to form here. In 
others, montmorillonite continued to form fora 
short distance, and a montmorillonite zone was 
thus originated. In these veins the respective 
ironts were separated. 

After the alteration of plagioclase stops, the 
unly effect of the fluids on the wall rock is the 
partial phlogopitization of biotite and replace- 
ment of hornblende by siderite, kaolinite, and 
nontronite, until the partial pressure of COs 
becomes too low even for this alteration (fringe 
front). 

(5) Petrographic evidence shows that no 
mutual replacement has taken place to any 
significant extent among the various secondary 
minerals. These repiaced only the primary 
constituents of the rock, chiefly plagioclase. 
More than one alteration mineral can be present 
in the same pseudomorph because of simul- 
taneous or successive replacement of the 
common host. For instance, the association of 
hydromica, orthoclase, and kaolinite in the 
orthoclase zone is due to expansion of fronts 
and encroachment of zones before the complete 
alteration of plagioclase. Simultaneous crystal- 
lization is exemplified by kaolinite and mont- 
morillonite. 

(6) The process of zone formation is based 
upon the hypothesis of outward migration of 
the fluids, as opposed to ionic diffusion in a 
static solute. The pressure gradient that drove 
the fluids away from the conduits was stronger 
than the concentration gradients of the indi- 
vidual ions. 

Siliceous, alkaline residual solutions leaving 
the alteration halo could travel long distances 
in the country rock and present no problem of 
disposal. 

(7) Formation of other patterns of zoning 
can be explained by differences in the composi- 
tion of the original fluids, which depended 
largely on the degree of hydrothermal fractiona- 
tion. Low potassium content may be the cause 
of the absence of the orthoclase zone (pattern 
V) or of the orthoclase and hydromica zones 
(pattern IV). The simpler patterns, where the 
kaolinite-montmorillonite alteration is poorly 
developed or missing, may be the result of 
lesser fractionation of the hydrothermal 


column. Pattern II requires further study. The 
hypothesis of limited fractionation agrees with 
the fact that the simpler patterns are more 
closely related to Laramide igneous stocks. A 
phase of greisenlike alteration was found in the 
Jamestown district, associated with pattern ITI. 

(8) No clear-cut relationship has been recog- 
nized between hydrothermal alteration and ore 
type and Laramide intrusions. In _ general, 
patterns I and II are associated with lead- 
silver-zinc deposits and also with somewhat 
older intrusives. Patterns V and VI are related 
to pyritic gold telluride and tungsten veins and 
to slightly alkalic dikes. Pattern III is related 
to gold telluride and fluorspar deposits and to a 
partly alkaline (saturated) composite stock. The 
alkalic (not alkaline) character of the alteration 
increases toward the northeast; this roughly 
parallels a similar change in the Laramide 
intrusions (Lovering and Goddard, 1950). 

No direct connection was detected between 
ore type and wall-rock alteration. This is in 
agreement with the conclusion that the two 
represent different phases of the hydrothermal 
processes. The type of ore may be largely a 
matter of availability of a particular metal at 
the moment of the process. No difference in 
pattern of alteration was found in districts 
where more than one ore type was present, as 
shown in the Boulder County tungsten and 
gold district. 

(9) The writer believes the explanation here 
offered for the origin of the zoning in the Front 
Range mineral belt can be applied to other 
districts, not only to those with similar char- 
acteristics, such as the Butte district, but also 
perhaps to others with more complex alteration 
such as the porphyry copper deposits. In these 
deposits, however, complicating structural 
features have to be taken into account. 

(10) The only supergene phenomena detected 
in relation to the hydrothermal alteration are 
transport and precipitation of iron and man- 
ganese hydroxides in the pyritic veins, pre- 
cipitation of sulfates and carbonates, and clay 
transportation. No supergene clay formation 
was found. 


ADDENDUM 


This paper had been submitted for publica- 
tion before the writer became aware of Tooker’s 
recent publication (1956). Since Tooker’s work 
deals with the Central City—Idaho Springs 
area within the Front Range mineral belt, a 
brief review of its conclusions and a comparison 
with those given in this paper follow: 
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The host rocks in the Central City—Idaho 
Springs area are quartz feldspar gneisses and 
granites of Precambrian Tooker dis 
tinguishes four alteration least 
altered (fresh) rock; (2) weakly altered; plagio- 
clase partly or totally replaced by montmoril- 
lonite; (3) moderately altered, soft clay-min- 
eral rock with primary textures preserved or 
indistinct; subdivided according to the pre- 
dominant clay mineral in (a) montmorillonite- 
rich rock, (b) kaolinite-rich rock, and (c) illite- 
rich rock; (4) intensely altered wall rock, 
original textures obliterated, sheared; silicified, 
sericitized, and pyritized. 

This zone sequence corresponds fairly well 
with that of the Nederland-Gold Hill area. Zone 
2 resembles the present writer’s montmorillonite 
zone, Whereas subzones 3a and 3b correspond to 
the outer and inner parts, respectively, of the 
kaolinite zone, which in many veins of the 
Boulder County area contains abundant 
montmorillonite in the outer half. Zone 3c is the 
hydromica zone. The silicified, innermost zone 4 
is lacking in the Nederland-Gold Hill district 
(exception, Quaker City vein), unless the wall 
rock crisscrossed by siliceous veinlets next to 
the main vein may be considered as such. 
Zone 4+ may perhaps be interpreted as a late- 
stage feature related to vein filling rather than 
to the argillic alteration stage. The orthoclase 
zone has not been recognized in the Central 
City-Idaho Springs area, the zone sequence thus 
corresponding to the present writer’s pattern V. 
In partial agreement with Sales and Meyer’s 
(1948) findings at Butte but contrary to the 
rule in the portion of the Front Range mineral 
belt studied by the present writer, potassium 
feldspar “is altered to kaolinite and perhaps 
other clay minerals” and “quartz is in part 
corroded and recrystallized” in zone 3 (Tooker, 
1956, p. 351). 

Tooker concludes that in the Central City- 
Idaho Springs district most wall-rock alteration 
preceded the sulfide phase of hydrothermal 
activity. 


age. 


zones: (1) 
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of an earlier generation. 


them dealing with the Scottish Highlands, have 
drawn attention to the widespread complex 
and superposed folds in tectonites (see, for 
instance, Sutton and Watson, 1954; King, 1955; 
King and Rast, 1955; Ramsay, 1956; Rast and 





Two simple kinematic models of folding in tectonites have been 
grams prepared to show ideal geometrical transformations of a set of parallel surfaces 
(S) as a result of superposition of one generation of folds on another. 

It was found for both models of folding that a second generation does not necessarily 
form folds with constant trend and plunge. These folds depend for their axial orientation 
on the attitude of S in which they form. It was found also that cylindroidal folds of a 
second generation are necessarily smaller and less persistent in axial direction than those 
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adopted and dia- 


Geometry of rotation of early formed B lineations by later folding is different for the 
two models of folding adopted. Where the second folding is by flexural slip of S, the 
earlier lineations follow small-circle paths in projection, as described by Sander. Where S 
is slip-folded by the second deformation, earlier lineations ideally follow great-circle paths 
in projection. Some departures from this ideal geometry to be expected in nature are 


listed. 
CONTENTS 
TEXT Figure Page 
sie 2. Geometry of acylindroidal foldsystem.... 94 
, fs 3. Relationship between S’ and 7S maxima in 
Introduction Beeb ee deca ses enww nsec esenee. 91 Sede OP ee ee 95 
Pete MOONEE, 55555 5 oe tans csacsst es 92 4. Structural geometry of a field of cylin- 
Models of folding............... ++. 9 re hc 96 
Kinematic geometry of the two models..... 92 5. Geometry of superposition of later folds on 
Geometry of cylindroidal fold systems...... 92 earlier cylindroidal folds. . . . 97 
Superposed deformations. . . . . 93 6. Geometry of 5 eg ga of p’ folds on 
Interfering or superposed sy stems of folds. .. & the 8 folds of Fi igure errr 99 
Geometry of superposition of folds......... 95 7. Geometry of rotation of a lidar ateantues 
Geometry of rotation of B structures by by later folding Pee ae 100 
flexural slip folding. . 98 8. Geometry of superposition of slip folds on 
Geometry of rotation of B structures by ‘slip a cylindroidal fold (axis 8 = B) by non- 
felling... ..-.-. 100 affine slip on S’ normal to B’. 101 
Departures from ideal geometry which may 9. Diagrammatic projection of field in Figure 
heespected if MAUS... 6. 5 osc ee 102 102 
Conclusions... ......-+0000ee eee e cece eens 105 10, Synoptic diagrams for data from Figure 9... 103 
References cited. ........ 2.66660 0eesee scene 106 11. Geometry of intersection of a group of 
tautozonal S’ surfaces with S on a cy- 
ILLUSTRATIONS lindroidal 6’ fold Asie . 104 
12. Rotation of plane of rotation of an earlier 
Figure Page linear structure (B) where nonaffine slip 
1. Kinematic geometry of flexural slip and on S’ is accompanied by elongation 
slip folding. . : : along either B’ or a’.. 105 
INTRODUCTION Platt, 1957; Clifford, 1957; Clifford et al., 
. ‘ . OS7- Jeiec , Os “hic 
Some recently published papers, many of — 1957; Weiss and McIntyre, 1957). This paper 


examines the special problems arising where 
a single set of parallel S surfaces (bedding or 
foliation) contains more than one generation of 
folds. 

Successive generations of folds may result 
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from different phases of one orogenic movement 
or from unrelated deformations separated 
widely in time. In the following discussion of 
geometry, no distinction is drawn between 
these. Examples are considered in which two 
fold systems are formed in two successive acts; 
however no genetic relationship between the 
acts is implied. Excluded from the following 
discussion of ideal geometry are complex fold 
systems oi the type described by Clifford ef al. 
(1957), who considered them products of 
essentially inhomogeneous deformation caused 
by mutual interference, during formation, of 
large folds. 

Geometrical transformations shown in the 
figures of this paper are ideally simplified and 
therefore may not correspond in detail to 
natural deformations. However, the close 
resemblance between ideal and natural geom- 
etry observed by the writer has encouraged 
him in the belief that a detailed study of the 
ideal cases will increase understanding of less 
regular natural geometry. In particular, such 
a study will help separate the relationships 
between axes of superposed folds which are 
geometrically possible from those which are 
geometrically impossible. All projections illus- 
trated are from the lower hemisphere of an 
equal-area projection. 

The writer gratefully acknowledges the 
helpful criticism and advice received from 
Dr. Iris Borg of Princeton and Dr. F. J. Turner 
of the University of California. 


Foips 1N TECTONITES 
Models of Folding 


Most geologists recognize two mechanisms 
of folding in tectonites and believe that most. 
folds probably are formed by a combination 
of these mechanisms together with mechanisms 
less easily defined. The two principal models of 
folding are as follows: 

(1) Folds in a body with layering S are 
formed by flexure or buckling of individual 
competent layers combined with simultaneous 
slip of layer over layer. This type of folding 
corresponds to the generally accepted model of 
flexural slip folding (see Knopf and Ingerson, 
1938, p. 155; Turner, 1948, p. 173-174). 

(2) Folds in a body with layering S are 
formed by nonaffine laminar flow or slip on a 
second set of surfaces inclined to S. Such a 
mechanism implies that, under certain condi- 
tions of deformation, layers in rocks lack co- 
hesion and become incapable of transmitting 


stresses, and that surfaces between individual 
layers do not act as surfaces of preferential 
slip. Carey (1954) has outlined features oj 
“rheid” behavior in rocks which resemble 
features of the generally accepted model oj 
slip folding, as this second type of folding js 
sometimes termed (see Knopf and Ingerson, 
1938, p. 157-158; Turner, 1948, p. 165-166). 
Ideally, slip folding, like rheid folding, has 
geometrical features like those of laminar flow 
in a fluid, and therefore it is sometimes termed 
flow folding or flowage folding; but these terms 
can denote a less homogeneous type of deforma- 
tion—more akin, perhaps, to turbulent flow in 
a fluid—with a correspondingly less regular 
geometry (Bain, 1931). 

It is not suggested that all folds in tectonites 
are of one of the two kinds just described. Most 
iolds, where considered in successive layers of 
different competence, combine features of both 
kinds and are formed by the two mechanisms 
acting together. In particular, neither mech- 
anism by itself accounts for the elongation 
along fold axes that is so common in tectonites. 


Kinematic Geometry of the Two Models 


The geometry of flexural slip folding is pro- 
jected in Figure 1a. ‘he slip suztace is S, the 
axis of slip is B (kinematic B axis), and the 
direction of slip is a (kinematic a axis). For 
opposing senses of external rotation of S about 
the fold axis (8), S; and Se are attitudes oi 
S. In S; and S»2 the directions: of slip are re- 
spectively a; and ae. For flexural slip folds, 8 is 
parallel to B (8 = B). 

The geometry of slip folding of S by non- 
affine slip on S’ is projected in Figure 1b. The 
internal rotations of S for opposing senses of 
slip on S’ are given by S; and Se. S is folded 
about its line of intersection with S’ (@’) irre- 
spective of the axis of slip in S’ (B’) or the 
direction of slip in S’ (a’). In slip folding, for 
all but special cases of deformation (for in- 
stance, where S intersects S’ either parallel to 
B’, so that 8’ = B’, or parallel to a’, so that 
no folding occurs), 8’ is not parallel to B' 


(8’ ~ B’). 


Geometry of Cylindroidal Fold Systems 


The ideal concept of a fold axis as the recti- 
linear generator of a fold (Clark and McIntyre, 
1951, p. 594) is only statistically fulfilled, 
because even regular folds do not persist un- 
changed in profile for a great distance along 
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FOLDS IN TECTONITES 


their axes. But, whether it is a flexural slip 
jold or a slip fold, and however complicated it 
may be in profile, a fold may be described as 
cylindroidal within a given field if the folded 





Figure 1.—KiINeEMATIC GEOMETRY OF FLI 


set of S surfaces is statistically tautozonal about 
a rectilinear axis. Such a statistically defined 
jold axis is termed a @ axis by Sander (1948, 
p. 132-146). The geometry of such a fold and 
its associated B lineations is summarized in 
Figure 2. 

In most fields of cylindroidal folding, axial 
planes of individual folds are subparallel. 
Where present, this parallelism is demonstrated 
in projection by a symmetrical disposition of 
poles of S (7S) within the wS diagram (Fig. 2). 
In the plane of the girdle of wS there may be 
one elongated maximum (Fig. 3a, upper) denot- 
ing a considerable appression of the folds 
(Fig. 3a, lower). The center of gravity of this 
maximum defines the pole of the axial plane 
(rS’) for the field as a whole. In most tectonites 
where folding is advanced, geometry approxi- 
mates this example. But, where folds are less 
appressed (Fig. 2), limbs of folds may {all into 
two sets, so that mS for each set forms a dis- 
tinct maximum in the 7S girdle (Fig. 36). The 
axial plane for the whole field (S’) is one of the 
symmetry planes of the mS diagram. Most folds 
found in tectonites have hinges of much smaller 
extent than limbs (the folds are similar rather 
than parallel or concentric), so that fewer 
measurements attitudes in the field are 
gathered from hinges than from planar limbs 
But the geometry of aS girdles can vary greatly 
in nature, especially where folds are consist- 


of 


Figure 2 
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ently overturned and have limbs of unequal 
length. 

Extension of a field such as that shown in 
may result in patterns of preferred 





EXURAL SLip FOLDING AND SLIP FOLDING 


orientation of S and B (B lineations, mullions, 
small folds, and so on) of decreased regularity. 
Where, for instance, the poles of S do not lie 
in a clearly defined girdle, or B structures do 
not lie in a single maximum, the field can be 
considered noncylindroidal, and the field must 
be divided into subfields in which folding is 
cylindroidal. Only by such a procedure can the 
geometry of noncylindroidal fields be estab- 
lished. Mapping of contacts alone is inadequate. 

The geometrical features of most systems of 
cvlindroidal folds can be summarized as follows: 
(1) The axes of folds in the system are statisti- 
cally rectilinear (this is implicit in the term 
cylindroidal). (2) Less common, but of more 
fundamental significance, is a statistical paral- 
lelism of axial planes of folds in the system. In 
slip folds this plane is the flow or slip surface. 


SUPERPOSED DEFORMATIONS 
Interfering or Super posed Systems of Folds 


Mutually interfering fold systems have long 
been recognized in some greatly deformed 
bodies of rock. They have been attributed to 
superposition of folds either in independent 
orogenies or in successive phases of one deforma- 
tion. (For recent discussion of the geometry of 
layered bodies so folded, see Ramsay, 1956; 
Clifford et al., 1957; Weiss and McIntyre, 1957.) 
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FIGURE 2. 


Some geometrical features of surfaces affected 
by mutually interfering fold systems are dis- 
cussed with reference to idealized examples in 
which the following conditions are assumed: 

(1) The first folding is evlindroidal and is the 
result of flexural slip (Fig. 4). The left-hand 
portion of Figure 4 is a map of a large fold with 
plunging axis (8) and a dipping axial plane. 
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—GEOMETRY OF A CYLINDROIDAL FoLp SysTEM def 
Figures 2a, 2b, 2-—projections of S; Figures 2d, 2e—projections of +S; Figures 2/, 2¢—projections 
of B. The lower projection in each column is the contoured equivalent of the one above. 
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In subfields I-IV (outlined on the map by 
broken lines), S is macroscopically planar, as | 
shown by the strong single maximum in each 


of the mS diagrams inset on the map. In each . 
of these fields weak folding about @ is shown by | ,,. 
elongation of the maximum into a partial girdle |, 
with axes respectively 8; to Bry. Where data | ,, 


from these four diagrams are compiled into a 
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gnoptic diagram (Fig. 4c), a pronounced girdle 
appears with axis 8. Likewise the four statisti- 
ally defined attitudes of S in fields I-IV (re- 
spectively S; to S;y) intersect in 8 as shown in 


a 


Figure 4g. Because deformation is by flexural 
slip, small-scale linear structures (small folds, 
mullions, B lineations, and so on) are present 
parallel to B = 8. In each subfield these lie in 
a maximum coincident with the corresponding 
3 (respectively 8;-8,;, in Figs. 4a-4d) and 
elongated in the plane of the. statistically 
defined S for that subfield (respectively S,—5; , 
in Figs. 4u-4d). Where plotted synoptically for 
the whole field, these B structures lie in a nearly 
equidimensional maximum (B in Fig. 4/) co- 
incident with @. 

(2) The geometry of superposition of a 
second generation of folds in S is investigated 
by geometrical operations in projection. The 
two cases considered are where the secon 
generation forms, first, by flexural slip of S, 
and second, by nonaffine slip on a set of sur- 
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faces (S’) oblique to 8. For both models of 
folding, axial planes of individual folds of the 
second generation are assumed to be subparal- 
lel. Such a relationship has been observed in 





b 


Figure 3.—RELATIONSHIP BETWEEN S’ AND 7S MAXIMA IN FOLDED Rocks 
Figure 3a shows symmetrical folds with closely appressed limbs; 3b shows ‘fopen”’ 
symmetrical folds. 


natural examples of superposed folds (sce 
Sander, 1948, p. 177-178; Clifford ef al., 1957, 
p. 17-19; Weiss and McIntyre, 1957, p. 596- 
600). 


Geometry of Superposition of Folds 


A second deformation of a body in which S 
is folded generally does not form a new genera- 
tion of folds with rectilinear subhorizontal fold 
axes. For each subfield in which S is planar 
alter the first folding, a rectilinear fold axis 
can form; but none of these axes is subhori- 
zontal because for the general case of a plunging 
first fold (Fig. 4), no horizontal attitudes of S 
remain anywhere in the field. The geometry of 
superposition of one generation of folds on 
another is shown in Figure 5; in Figure 5a, 
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EOMETRY OF A FIELD OF CyYLINDROIDAL FOLDING 
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upper, S is folded about a subhorizontal axis 
(8). The geometry of this folding is projected 
in Figure 5a, lower, in terms of several attitudes 
of S measured on opposite limbs and hinges of 
folds. Because folding is relatively open, two 
attitudes of S predominate, corresponding to 
the two kinds of fold limb. The axial planes of 
the folds (Sg) are shown vertical. Figure 50, 
upper, shows the geometry of the @ folds after 


superposition of a new generation of . folds 
(either flexural slip or slip folds) with sub- 


parallel axial planes (.S’). The axes of the second 
generation of folds (8’) lie in S’; but the orien- 
tation of B’ in a particular subfield is parallel 
to the intersection of S with S’ and cannot be 
rectilinear. Figure 5), lower, projects the geom- 
etry of this structure. No vestige remains of 
the 8 axis in Figure 5a, lower, because all atti- 
tudes of S in the field have changed, and pre- 
vious 6 intersections have been dispersed. The 
B’ axes and associated lineations and small 
folds are rectilinear in subfields where S was 
planar after the first folding, but for the field 
as a whole, there are two maxima of 9’ struc- 
tures, corresponding to the intersections of S’ 
with the two planar preferred orientations of S 
(Fig. 5a, lower), spreading into a girdle in the 
plane of S’, as shown by the diagrammatic 
contours in Figure 5b, lower (see Weiss and 
McIntyre, 1957, Fig. 12). These two maxima 
correspond to §’ axes with greatly different 
trend and appreciably different plunge, al- 
though they belong to the same generation of 
folds. 

The geometrical operations just described 
can be performed also on the plunging fold 
shown in Figure +. Superposition on this struc- 
ture of a second generation of folds with sub- 
parallel axial surfaces (S’) gives rise to the 
geometry shown in Figure 6. In subfields where 
S was planar (subfields I-IV), cylindroidal 
folds with axes respectively 8, to 8, y are formed 
(see diagrams inset on map). These axes are 
oriented differently in each subfield, reflecting 
the differently oriented intersections between 
S and S’. For the sake of simplicity, folds of the 
second generation are assumed to have closely 
appressed limbs, so that aS defines a single 
maximum within each girdle subnormal to .S’ 
(Fig. 3a). The pattern of preferred orientation 
of wS for the whole field does not define a 
rectilinear 6’. There must be a strong maxi- 
mum of wS at 7S’; but the presence of girdles 
in different planes for subfields I-IV produces 
an irregular pattern in which no (’ axis for the 
whole field can be determined, and in which the 
previously existent 8 maximum is distorted 
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or, where folding is advanced, completely di: 
persed. This condition is shown diagrammat. 
cally in Figure 67, in terms of the planes of th. 
four mS girdles intersecting in aS’. 
Small-scale structures parallel to 8’; to ¢', 
are B’ structures only where the second de 
formation is by flexural slip. Whatever the: 
kinematic significance, however, these stn. 
tures lie as shown in Figures 6a—6d re spectivel 
for the four subfields. The maximum in ea¢ 
subfield coincides with the corresponding ; 
axis. For the whole field, the small-scale stryc. 
tures lie in a girdle with maxima at the line 
of intersection of S’ with the planar S$ in sub. 
fields I to IV (Fig. 6)). 
A general rule that can be established from 
Figure 6 with respect to the scale of successive 
generations of folds is that the largest folds o' 
the generation, where cylindroidal, 
must be smaller than those of the first genera. 
tion, because they can be formed only on planar 
limbs of the first generation. In a field with 
two generations of folds, it is therefore to b 
expected that the largest or major structures 
are the first formed, the second-generation 
folds take the form of smaller-scale or minor 
structures superposed on the limbs of the large 
folds. Such structural relationships are wide- 
spread in the Scottish Highlands and occur in 
many other parts of the world. 
The geometry of S outlined holds equally 
well for folds formed by flexural slip and by 
slip, but, where the geometry of rotation of 
early formed B structures (in particular, linea- 
tions and small folds) is considered, significant 
theoretical differences can be seen for the two 
models of folding. 


second 


Geometry of Rotation of B Structures by Flexural 
Slip Folding 


The geometry of rotation of an older linea- 
tion by flexure of the surface in which it lies 
is well known as a result of the work of Sander 
(1948, p. 171-172). Let a B lineation lie ina 
planar S. Then, flexure of S about another axis 
8’ (Fig. 7a) rotates B as shown. On the attitudes 
of S resulting from the folding (5), Ss, and so 
on), B maintains a constant angle with #’ and 
so follows a small-circle' path on the projection 
with 8’ as its center (B,, Bs, and so on). In each 
of the subfields I to IV in Figure 6, the angle 
between B (= 8) and 8’ has a different value, 


' The term small circle is used in this paper to 
denote any path (other than a great circle) that is 
circular on the sphere of projection, regardless of 
the attitude of this circle. 
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so that the small circle of rotation of B struc- 
tures about §’ is different for each subfield. 
These circles are shown respectively by C, to 
C,y in Figures 6e-6h. In Figure 6k they are 
shown superposed. For the field as a whole 


Geometry of Rotation of B Structures 
by Slip Folding 


Figure 8a shows a body containing a simple 
open fold in S with B lineation parallel to its 





FIGURE 7.—GEOMETRY OF ROTATION OF A LINEAR StRucTURE (B) BY LATER FOoLpING 


Figure 7a shows flexural slip folding with axis 8 = B; 7b shows slip folding by nonaffine slip 
on S’ with axis Bp’ + B’. 


there is no simple path of movement of B struc- 
tures. Thus, the second deformation by flexural 
slip folding tends to disperse the initial maxi- 
mum of B structures (Fig. 4/) irregularly. In 
each subfield, B structures can be “unrolled”’ 
about the respective 6’ in the manner outlined 
by Sander (1948, p. 174-178); but no single 
axis of unrolling will serve for the field as a 
whole. Also, in order to reconstruct their initial 
preferred orientation, B axes should not be 
unrolled to a horizontal orientation as sug- 
gested by Sander, but only to the point of 
intersection of the small circles of rotation. 
This point coincides with the initial preferred 
orientation of B structures shown in Figure 4/. 
A simultaneous unrolling of the attitudes of S 
in which individual B structures lie should 
restore these to their preferred orientation 
before the second folding. 

It can be seen from Figure 6 that a second 
deformation by flexural slip must be, on the 
scale of the whole field, triclinic. In individual 
subfields, deformation can be monoclinic or 
orthorhombic, but the planes of symmetry for 
the subfields are mutually inclined, so that 
none hold for the field as a whole. 


axis (G). In Figure 8) the same body is shown 
after nonaffine slip on S’ (one of the vertical 
sides of the body) in direction a’, normal to 
B’ (one of the horizontal edges of the body). 
Different attitudes of S become slip-folded 
according to the geometry shown in Figure 3b, 
and cylindroidai folds form on the limbs of the 
8 folds with axes (8’) parallel to the intersec- 
tions of S’ with the effectively planar attitudes 
of S. The geometry of rotation of B in one such 
cylindroidal field is shown in Figure 7). As S 
is internally rotated to S, or Ss by opposite 
senses of slip, B moves respectively to B, and 
By. The geometry of this deformation has been 
discussed elsewhere (Weiss, 1955, Fig. 1). 
B moves along a great-circle path of the pro- 
jection passing through a’. This path remains 
constant for any field in which S’ is planar and 
B is rectilinear before deformation. 

The geometry of a field such as that in 
Figure 8 is projected in Figure 9. The folded S 
of Figure 8a is shown in Figure 9a. S; to S¢ 
correspond to six subfields in which S is effec- 
tively planar on the first fold with axis 6. The 
surface of nonaffine slip induced by the second 
deformation is S’. In each of the six subfields, 
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cvlindroidal slip folds are formed with axes 
(3': to 8’s respectively in Figs. 9b-9g) parallel 
to the intersections of the respective planar 
preferred orientations of S with S’. In this 


attitudes are plotted, and in 10c these inter- 
sections are contoured to give a $’ diagram. 
Plainly, the single maximum in this diagram is 
an oversimplification of geometry and conceals 





FiguRE 8.—GEOMETRY OF SUPERPOSITION OF SLIP FOLDS ON A CYLINDROIDAL FOLpD 


(Axis B = B) 


fashion the initially homogeneous field of 
Figure 9a becomes divided into six homogeneous 
subfields (F; to Fs) with geometry as shown 
respectively in Figures 9b-9g. In each subfield 
the poles of S (circled dots) come to lie in a 
girdle normal to the corresponding 6’. This 
geometry is independent of the direction of 
slip (a’) in S’. 

Internal rotation of S in each subfield is 
accompanied by internal rotation of minor 
linear structures (B) initially parallel to 8 
(Fig. 9a).-In each of the subfields F; to Fg, 
B rotates toward a’. The unique plane defined 
by B and a’ is constant for the whole field, so 
that the path of rotation of B is along the same 
great circle of the projection in all subfields 
and is independent of the orientation of 9’. 
The amount of rotation of B in any subfield 
depends upon the amount of internal rotation 
of the S in which it lies (crosses in Figs. 9b-9g). 

The data from Figure 9 are compiled into 
synoptic diagrams in Figure 10. Figure 10a 
shows all planar attitudes of S from Figures 9b 
to 9¢. In Figure 10) the intersections of these 


BY NONAFFINE SLIP ON S’ NorMAL to B’ 


the true relation of 6’ axes in the whole field. 
Only by preparing §’ diagrams for the six sub- 
fields in which folding is cylindroidal can the 
true geometry of S be established. Figure 10d 
is a plot of all small-scale B structures (dots) 
and all small-scale 8’ structures (circles) from 
Figures 9b to 9g. For the field as a whole, these 
small-scale structures lie in two distinct girdles 
intersecting in a’. Because the @’ structures lie 
in the plane of S’, the orientation of B’ (lying 
in S’ normal to a’) can be found. 

The results of this analysis can be applied to 
the field shown in Figure 6. The great circles of 
rotation of B must be the same for all subfields 
as shown by G, to G;y respectively in Figures 
6e-6h. These girdles intersect S’ at a’, thereby 
fixing the orientation of B’ in each subfield. 
The girdle of B structures for the field as a 
whole (Figure 6/) coincides with girdles for 
individual subfields. 

Comparison of Figures 6k and 6/ shows that 
geometry of rotation of older lineations by 
flexural slip folding is markedly different from 
geometry of rotation by slip folding. Unrolling 
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FIGURE 9.—DIAGRAMMATIC PROJECTION OF FIELD IN FIGURE 8 


of older lineations after slip folding is perhaps 
possible, because the path of movement is 
known; but, as can be seen from Figure 6/, no 
clue remains to the initial preferred orientation 


of B. 


Departures From Ideal Geometry Which May 
Be Expected In Nature 


The geometry discussed is ideally perfect. 
Obvious departures from this geometry can be 
predicted in natural fabrics, as follows: 

(1) It has been assumed in this study that 
the slip surface (S’) is planar. In structures of 
this kind at Loch Leven (Weiss and McIntyre, 


1957, p. 590) S’ surfaces vary slightly in orien- 
tation, so that poles of the axial planes lie ina 
girdle with one very strong maximum. This 
range in orientation can be explained in two 
ways: in a large field S’ is more correctly repre- 
sented by a set of intersecting surfaces with a 
limited but marked range in orientation. Where 
such a set of surfaces intersects pre-existent 
folds, as shown in Figures 11a—11d, the inter- 
sections of S’ with S lie in a sheaf-shaped girdle 
narrowest at B’ (Fig. 11d). A sheaf-shaped 
girdle of this kind observed at Loch Leven was 
interpreted in this fashion (Weiss and 
McIntyre, 1957, Fig. 13). According to the 
second explanation slip on S’ has developed 





“a7 


orien- 
> ina 
This 
1 two 
‘epre- 
‘ith a 
Vhere 
stent 
inter- 
rirdle 
aped 
) was 

and 
» the 
oped 





SUPERPOSED DEFORMATIONS 103 





Figure 10. 


into flexural slip of S’ (cf. Carey, 1954, p. 93- 
93). In this way an clement of external rota- 
tion about B’ is added to all structural elements. 
Such rotations follow small circles of the pro- 
jection and have the geometry shown in 
Figure 7a. 

(2) The field at the end of the first deforma- 
tion may be noncylindroidal, so that a second 
deformation by slip folding could not reproduce 


Synoptic DiAGRAMs FoR Dara From Ficure 9 

Figure 10a shows S planes (full lines) and S’ (broken line); 105 shows intersections of S planes; 10¢ 
shows intersections contoured, contours 1-2-3-4-5 per cent per 1 per cent area, 10d shows girdles of B 
dots) and 8’ (circles). 


the ideal geometry of Figure 7. Axial planes of 
B’ folds are still subparallel, but geometry of 
internal rotation of B structures is not constant 
throughout the field, because the initial pre- 
ferred orientation of B is not everywhere the 
same. This difficulty should be overcome by 
dividing such a field into smaller subfields. 

(3) Slip on S’ may not be constant in direc- 
tion throughout a field in which 5S’ is planar. 
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FicgurE 11.—GErOMETRY OF INTERSECTION OF A Group OF TaAUTOZONAL S’ SuRFACES WITH 
S on A CyLINDROIDAL 8 FOLp 
Figure 11a shows geometry of S in 8 fold (dots are 7S); 116 shows S’ surfaces with a limited range of 
orientation, intersecting in B’; 11¢ shows intersections of S and S’ (200); 11d shows ’ diagram, con- 


tours, 2!9-5-10-15-20 per cent per 1 per cent area. 


Geometry of internal rotation of B is different 
for different directions of slip (a’ axes). Axial 
pianes of B’ folds remain parallel to S’, and 
8’ axes remain parallel to the intersections of 
S and S’. Deformation is triclinic. 

(4) As De Sitter (1956, p. 97) has pointed 
out, a shortening normal to S’ is probably 
associated with slip on S’ in the natural equiva- 
lent of slip folding. For deformation to remain 
monoclinic, this shortening must be accom- 


panied by an elongation of the deformed body 
in the plane of S’, either in a’ or B’. Such sub- 
sidiary elongations do not significantly affect 
the geometry of internal rotation of S. They 
result in a more pronounced appression of fold 
limbs than could be achieved by the simplified 
model of slip folding. Similarly, the effects of 
these elongations on the geometry of rotation 
of B are minor so long as deformation is effec- 
tively homogeneous. Whether elongation is 
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along a’ or B’ is immaterial; the effects are the 
same—namely, that the plane of rotation of 
B structures approaches parallelism with S’ by 
rotation about a’. The geometry of rotation is 
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the following general conclusions can be 
drawn: 
1. Where one generation of folds is super- 


posed on an earlier generation, the geometry of 





Figure 12.—ROTATION OF PLANE OF ROTATION OF AN EARLIER LINEAR STRUCTURE (B) WHERE 
NONAFFINE SLIP ON S’ 1s ACCOMPANIED BY ELONGATION ALONG FE1THER B’ or a’ 


shown in projection in Figure 12. BP, is the 
plane of rotation of B by nonaffine slip on S’ 
without elongation along a’ or B’. Where such 
an elongation occurs, the plane of rotation of B 
rotates about a’, occupying successively such 
positions as BP» and BP3 as it approaches 
parallelism with S’. 

Any of these complicating factors may 
operate in natural fabrics, thereby contributing 
toward divergence from ideal geometry. 


CONCLUSIONS 


The ideal geometry of folding adopted in the 
analyses discussed must occur rarely in nature. 
Most folds in tectonites probably have an inter- 
mediate kinematic geometry combining ele- 
ments of the two models. From the analyses, 


S ceases to be expressible in terms of a recti- 
linear axial trend of folding, because the earlier 
folds are bent, and the later folds develop with 
variable trend and plunge. 

2. In general, cylindroidal folds of a later 
generation are on a smaller scale than those of 
an earlier generation, because they form only 
in fields where S remains planar after the earlier 
folding. 

3. The geometry of rotation of earlier linear 
structures by later folding depends upon the 
kinematics of the later folding. For most kinds 
of folding a lineation in projection follows a 
path intermediate between a small circle and a 
great circle. Kinematics of folding may be 
indicated by the resemblance of the natural 
path to either of the two ideal paths. 

4. The diagrams show that fields affected by 
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superposed folding show a high degree of geo- 
metrical regularity which is capable of being 
established and interpreted. 
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FAULTED RUNDLE (MISSISSIPPIAN) SECTION AT CROWSNEST PASS, 
ALBERTA, CANADA 


By SAMUEL J. NELSON 


The abnormally thick carbonate section of 
the Rundle group at Crowsnest Pass, Alberta, 
has interested geologists for many years. The 
Rundle here attains an apparent thickness of 
more than 5000 feet in contrast to near-by sec- 
tions most of which are less than 3000 feet. 
Although no obvious faults occur in this sec- 
tion, some geologists attribute the thickness to 
repetition, while others consider basinal deposi- 
tion as the answer. J. R. Cundill (1955) out- 
lined stratigraphic evidence illustrating fault- 
ing in the section, and the writer here presents 
paleontologic evidence in support of Cundill’s 
interpretation. 

The section discussed is exposed on the north 
side of Crowsnest Lake at Crowsnest Pass in 
the southern Canadian Rockies immediately 
east of the Alberta-British Columbia boundary. 
The Rundle here is a westward-dipping, super- 
ficially undisturbed sequence of carbonates 
underlain on the east by the Banff and overlain 
on the west by the Rocky Mountain formation. 

The Rundle group at Crowsnest Pass, as in 
most areas of the southern Canadian Rockies, 
can be divided into three formations which, 
in ascending order, are the Livingstone, Mount 
Head, and Etherington. In other areas these 
three formations generally aggregate between 
2000 and 3000 feet: the Livingstone (1000- 
1500 feet), of resistant light-colored crinoidal 
I'nestones and interbedded finer-grained car- 
bonates, of Osagean age; the Mount Head 
(600 feet), a recessive unit of Meramecian age 
consisting of dark limestone, dolomite, and 
shaly beds; and the Chesterian Etherington 
formation (700 feet or less) of resistant light- 
weathering dolomites and limestones. The 
reader is referred to Warren (1927), Beales 
(1950), Douglas (1953), Norris (1957), and 
Bostock, Mulligan, and Douglas (1957) for 
details of Rundle stratigraphy. 

Dawson (1886), who made the first reference 
to the Rundle of Crowsnest Pass, considered 
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that complex folding was the cause of the ex- 
traordinary thickness. Warren (1928; 1933) and 
McKay (1932) thought that their measured 
thicknesses of between 5000 and 5300 feet were 
caused by deposition and not structura! com- 
plications. Cundill (1955) concluded, from 
lithologic studies, that the succession was re- 
peated by three faults and possibly a fourth, 
and indicated that the true stratigraphic thick- 
ness is about 2700 feet.! The section was 
recently briefly referred to by Bostock, Mulli- 
gan, and Douglas (1957), who stated that the 
maximum thickness is 3900 feet. Faulting was 
not mentioned. 

In 1955 the writer with J. R. Cundill and 
two assistants spent 10 days making detailed 
fossil collections from the section. A study of 
faunas has shown that Cundill’s three main 
faults are probably correct, and that the fourth 
which he tentatively proposed is also present. 
Faunal evidence indicates that the stratigraphic 
thickness of the Rundle is between 2900 and 
3000 feet (Fig. 1). Thus Cundill’s excellent 
stratigraphic analysis is confirmed by paleon- 
tology. 

The writer wishes to express his gratitude to 
officials of The British American Oil Company, 
Ltd., Calgary, for financial backing of this 
project. Dr. P. S. Warren, of the University of 


! Cundill placed only the Livingstone and Mount 
Head formations in the Rundle group and stated 
that their combined stratigraphic thickness was 
about 2100 feet. If the Etherington (his Tunnel 
Mountain formation) is included in the group 
then the above figure of 2700 feet is obtained. 
Some confusion has been caused by the fact that 
his diagram shows a total thickness of 2750 feet 
for the Livingstone and Mount Head in contrast 
to the stated figure of 2100 feet. Cundill (Personal 
communication) mentioned that the latter figure 
was derived by assuming an arbitrary value for 
repetition on the fourth fault. Faunal evidence 
(Fig. 1) suggests that the figure is closer to 2300 
feet. 



















































































Daay Arjona “0g 





Se. ESBS SES SEES SSEBSBSEF LSS  BHEESBSEHSSSEsHESEEBESE SP SSSRESCELCSESHSE 
“_ &. >= aa WE a oF Cs NVOOU =s'a 6 = — a 6 — a oe a a — A= OVm SS Ve SS SSO . Fe — -.= = © > 
VLaaATY “SSVq LSANSMOAD AO (NVIddISSISSIJZ) dQOUD AIGNOY AUL AO NOLLVIANAOO—'] AAA 
oe 
Uo! yOuss0 y fe ee 
jjuog = 
UdUDM 1 Vv 
iy aida nina Bune) 0 SZ] 
(d3.N) S/swavoyo cnpnposd ¢ | ag, ee VRP 
(Mim) avgeyour7 4a’ os b q tea joussoy rae ; 
vos u vovouioyy7 ¢ } POaH 4unow fe-T I 
(Aowyrtsy) rrauys oyavoyousoyy/7 9 id | se: 
(944) unyouy ru wiryjhydoyopbhuy | | uoryOu04 x iy cs equmenn 
os! J = - 
WEE?) uoGursdyy yj IA II 7 ] 0S? - 
. \ x ; 3 
(1H) v0! jouw 404 mer LJ a 
: uinjunow 49204 Fes l ia 
i = 
Z Loc 
as sh | 
— oe = 
wd / T TI 0004 ooor 4 T 4 
si \y u 4 T TI 
7 x 00st = é ; 
\ \ [Thee 
se ’ ' a | 
} ae x . 
= me * 
7 CI I | 
= 19S'6¢'2 at PT ; oe € PS 4jNo4 => TSF eos ry 
TH ovo» i : 
) i Ur a i 
: a = oe Prt 005 
Z t'1'016'8.> # wi . jest 
{fF —- ~_ re 
fe) =r 2 eS ee ci ane oe a = 
é - <i _ — — vinee» oy 
| K/ ted 1295 419% eae | 
jo | i fess 
‘ | X == GA 
oi? / 61-O*= Tay 00%b oak dhiehG: Ancues GnedeD debits atietes  ccceee Sen - | 
, | ‘ I rm — Ss 2 | 
- b2 1 ones be ty 
« t | 
grz0 error al Tt -s a e ——_ —<—<— <<< ams om - - -~ L | 
| 4 ® nat bd 
be | | 
} yd | 
} } | 
| +4 | 
| 
| 
} L , | | 
| ITT} oss | | 
| a | 
| xe | 
| Sem | 
= | 
| ayy | 
| } 
| 5 PG, 4104 | 
| a epoaty ’ yuog 
| 
| | 


108 








NIVINQOW T3NNAL 





ALBERTA 


Group (MIssIssipP1iAN) OF CROWSNEST Pass, 


RUNDLE 


THE 


-—CorRRELATION OF 


1 


FIGURE 





SHORT NOTES 


\lberta, has guided the writer through most 
phases of this study. 
- Figure 1 shows the stratigraphic correlation 
of the fault slices diagnosed for the Rundle 
limestone of Crowsnest Lake with two well- 
known sections in Alberta—Mount Rae (Storm 
Creek), about 65 miles north of the Pass, and 
Tunnel Mountain, 45 miles farther northwest. 
Footages for the last section are those of Beales 
1950): those of Mount Rae are the writer’s. 
Portions of sections not pertinent to this dis- 
cussion are omitted; their deletion is indicated 
by jagged lines on columnar sections in the 
figure. 
Lithologic divisions and footages of the 
Crowsnest Pass Rundle (Fig. 1) are those of 
Cundill (1955). The base pf the Mount Head 
jormation in this paper has been raised 200 feet 
above that of Cundill so that the lithologies 
conform more closely with the formation at 
Tunnel Mountain and Mount Rae. Cundill’s 
‘Tunnel Mountain formation” is here called 
the Etherington formation (see Norris, 1957; 
Bostock, Mulligan, and Douglas, 1957) to 
conform with more recent nomenclature. 
Five diagnostic fossil assemblages are recog- 
nized in sections portrayed on Figure 1 and 
are designated as assemblages A to E, respec- 
tively. Spirifer banffensis Warren is the only 
jossil present in assemblage A. This easily 
recognizable species occurs 585 feet above the 
base of the Livingstone on Tunnel Mountain, 
suggesting a middle Livingstone correlation 
for assemblage A in fault slices 1 and 3. Species 
of assemblage B, geographically widespread in 
the Canadian Rockies, are characteristic of 
lower Mount Head strata. Cundill’s (1955) 
lithelogic evidence has been used in determining 
the fault position at the top of fault slice 1, for 
the stratigraphic interval between assemblage B 
in this slice and in superjacent slice 2 could 
represent the fossil assemblage range rather 
than repetition by faulting. Assemblage C is 
characteristic of the lower part of middle 
Mount Head although Echinoconchus biseriatus 
(Hall) and Dictyoclostus tenuicostus (Hall) are 
known to range some distance higher. Assem- 
blaze D consists of Lithostrotion arizelum 
(Crickmay). The holotype of this species is 
irom basal Mount Head strata at Lake Minne- 
wanka (Crickmay, 1955). Specimens almost 
identical with the type are abundant in upper- 
most Mount Head at Mount Rae. The writer 
considers that these stratigraphically higher 
specimens are either homeomorphs developing 
out of L. whitnevi Meek of middle Mount 
Head, or an expression of a reintroduced zonule 
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useful for local mapping. The stratigraphic 
position of the Crowsnest Pass L. arizelum with 
respect to assemblage C suggests that the 
containing strata probably correlate with 
uppermost rather than lowermost Mount Head. 
An interval of about 500 feet separates assem- 
blages D and E, while in other nearby sections 
this interval is 140-200 feet. Thus the enlarged 
interval suggests repetition by faulting. No 
diagnostic fossils occur above assemblage E, 
but the 400-foot section between E and _ the 
base of the Rocky Mountain formation is simi- 
lar to that of normal sections and suggests that 
faulting is absent. 
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AGE OF THE SAMARSKITE OF KISHENGARH, RAJASTHAN, INDIA 


By U. ASWATHANARAYANA 


There are three orogenic belts in Rajasthan 
also known as Rajputana), India—the 
Gneissic complex, the Aravalli orogenic belt, 
and the Delhi orogenic belt. “Each of the 
three orogenic cycles of Rajputana ends with 
adisplay of pegmatites, and, since the younger 
pegmatites also occur in the older rocks, it is 
not always easy to assign its proper age to any 
particular occurrence” (Holmes ef al., 1949, 
p.291). Apparently there are marked differences 
in the characteristics of pre-Delhi and post- 
Delhi pegmatites. The former are medium- 
grained, foliated, and disturbed by later move- 
ments; they contain few, if any, radioactive 
minerals. The post-Delhi pegmatites are coarse- 
grained, massive, and unfoliated and have 
yielded radioactive minerals. Holmes el ail. 
(1949, p. 293, 296) determined the ages of 
uraninite from Bisundni (735 + 5 m.y.) and 
monazite from Soniana (700-865 m.y.) in 
Rajasthan, both of which belong to the post- 
Delhi group. 

The present investigations extend the earlier 
work reported by the author (1956) and aim to 
date the samarskite and feldspar from the peg- 
matite of Bajrang mine located 13 miles east 
of Kishengarh town along the Kishengarh- 
Jaipur road (75° 2’; 26° 36’). This pegmatite 
is coarse-grained, massive, and contains radio- 
active minerals and is considered post-Delhi. 
It strikes N. 70° W. and dips 70° N. It consists 
of quartz, pink feldspar, usually microcline, 
muscovite, beryl, samarskite, and monazite. 
The feldspar does not occur in conspicuously 
large crystals. There is no evidence of segrega- 
tion of minerals in any particular zone. 

Spectrochemical and X-ray analyses have 
been made to supplement megascopic identifi- 
cation of the radioactive mineral (samarskite). 
The following are the spectrochemical data 
in terms of the elements and their concentra- 
tion in ppm.: 


Ba—145; Ce—4250; Co—240; Cr—29; Cu—40; 
La—400;  Mn—630; Nd—2000; Pb—2000; Sc— 
215; Th—18500; Zr—2650; Sr—320; V—121. 


¥Y and Yb are present as minor elements. The 
spectrochemical analysis was not conclusive. 
X-ray-diffraction analysis using goniometric 
set-up clearly indicated that the mineral is 
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noncrystalline and metamict. This observation 
coupled with the data on the d spacing of the 
X-ray powder pattern suggests that the mineral 
is samarskite 

Elaborate precautions were taken to avoid 
contamination in the process of powdering. In 
the case of feldspar, the diamond mortar and 
pestle, hammer, chisel, and stainless-steel sieve 
were first cleaned with triple-distilled water, 
immersed in nitric acid, rinsed with quadruple- 
distilled water, and tested for the absence of 
aciditv. The crushing is done in a dust-free 
room and on a table top vacuum-cleaned and 
covered with clean paper. 

The core of the mineral is separated from the 
rest, untouched by hand, crushed and sieved 
with a 60-mesh sieve. The minus 60 mesh 
powder is used in the investigation. 

Similar precautions were taken for powdering 
samarskite. The precautions described by 
Tilton et al. (1955) to avoid contamination 
from reagents, apparatus, and laboratory air 
have been scrupulously followed. Whenever 
possible, lead-free teflon beakers and hoods 
have been used in preference to pyrex glassware, 
and it is hoped that this step contributes to 
the accuracy of the data. The nature of the 
laboratory studies undertaken are given in the 
following schematic arrangement: 


Isotopic com- 


‘ position Two aliquots 
Unwashed— 
> from the same 
Dissolution . 
Total lead con- solution 
Feldspar centration 
Washed— Isotopic com- 
Dissolution position 
U concentration ‘ 
‘ ‘ Two aliquots 
Dissolution 1 Giga 
. ¢ ne same 
Samar- Pb concentra- : ; 
m . solution 
skite tion 
Dissolution | Pb isotopic 
composition 
The concentrations of uranium and lead 


in samarskite and the concentration of lead in 
the feldspar were determined by using isotope- 
dilution techniques (Tilton e al., 1955). Be- 
cause of the large quantities of lead in 
samarskite, about a milligram of the sample 
only was used for the analysis. 
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The mass-spectrometric procedures employed 


by Tilton ef al. (1955) for the isotopic analysis 
of uranium and lead have been utilized in the 
present investigation too. The isotopic com- 


position of lead was computed from about 
40 sets of peaks in each case. The concentra- 


TABLE 1.—SAMARSKITE—LEAD AND URANIUM 


Concentration of 


Isotopic composition 
total uranium 


Concentration of 
of toal lead Phrve 


Pb26/Pb24 = 5942 9.094 wg per 1.2 mg 


of sample 


119.2 wg per 1.2 
mg of sample 
16.28 


22.65 


Pb**6/Pb27 = 


Pb2¢/Ph2s = 


tions of lead and uranium were calculated from 
about 25 sets of peaks. The limits for the ac- 
ceptable data were the following: (1) 3 & 107"! 
amps. in ion beam for more than 2 hours, with 
less than about 0.5 per cent of intensity change 
per minute, and (2) the deflection on the data 
chart above 50 per cent of full scale (Chow 
and Patterson, unpublished). The probable 
absolute accuracy is within ~ + 0.5% for 
Pb**/Pb™ ratio, ~ +0.2% for Pb®*/Pb ratio, 
and ~ +0.3% for Pb®*/Pb™ ratio. 

The atomic percentages of Pb™, Pb’, 
Pb**, and Pb** were computed from the ratios 
in Table 1 and corrected for nonradiogenic 
Pb*® and Pb*’. The nonradiogenic lead ratios 
in the unwashed feldspar lead (Table 2) 
(Pb”*/Pb™ = 16.47; Pb*®?/Pb™! = 15.58) 
were used in the calculations. The isotopic 
composition of nonradiogenic lead is, however, 
not critical as the Pb®* atomic per cent is 
extremely low. The corrected atomic per cents 
were then converted into weights using the 
Pb**® content of the sample, as estimated by 
the use of isotope-dilution technique. 

Ages have been computed from the standard 
equations (Faul, 1954, p. 282, 294), and the 
following constants have been used in the 
calculations: 


A238 = 0.1537 K 10-9 Yr“; 
235 = 0.972 K 1079 Yr; 
K = 137.8 


The age spread is given below: 


Pb*®*/Pb6 age = 579 m.y. 
Pb”6/U2% age = 587 m.y. 
Pb*/U25 age = 578 m.y. 

The data in Table 2 can be used to calculate 
the hypothetical time at which the lead was 
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isolated from its uranium and thorium environ. 
ment and was introduced into the feldspar 
lattice. Assuming that the lead in the feldspar 
had a normal evolutionary history, one cap 
either calculate the time required to go from 
a normal lead today back into time into the 


TABLE 2.—TIsotoric Composition (ATOMIC Ratios) 
OF FELDSPAR LEAD 

Total le 

Sample Ph206/ P20 Ph26/Ph207, Ph206/ Pp20s pom ry 
ppm. 
Washed feld- 16.58 1.055 | 0.4267 | NLD. 

spar 

Unwashed 16.47 1.057 | 0.4285 68.3 


feldspar 


feldspar lead or the time required to go from 
the primordial earth lead to the feldspar lead. 
In the first calculation, two expressions from 
Collins ef al. (1953) were used to obtain two 
ages. The values for constants a (= 18.69) 
and 6 (= 38.82) in the above calculations 
were provided by T. J. Chow (unpublished), 
In the second calculation, the expression and 
constants of Patterson (1956, p. 231, 235) 
were used. 


The calculated ages are as follows: 


206 208 207 206 

— = 1200 m.y.;¢ = 200 my.;/ — ,— 

204 BOA =?" 904 ’ 204 
= 1400 my. 


Samarskite has given an unusually closely 
spaced age spread, and hence an age of 580 + 
20 m.y. can be given to it with confidence. It 
is believed that the errors in chemical deter- 
minations define the errors in the ages. Eckel- 
mann and Kulp (1957, p. 1121) listed lead 
isotopic data on five samarskites. The low 
abundance of Pb* in the present sample is o/ 
interest as it is among the lowest concentra- 
tions reported, if not the lowest. 

The gross discrepancy among the feldspar 
lead ages and the samarskite age can be best 
explained in the following way: The lead con- 
tained in the feldspar is abnormal in that its 
isotopic composition does not correspond to 
that of ore lead of the same age. The abnormal- 
ity consists in its having marked deficiency 
in Pb*®*® and marked excess of Pb®. 

The similarity in the isotopic composition 
of lead obtained from the acid-washed and 
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SHORT NOTES 


unwashed feldspars indicates that neither leach- 
ing subsequent to the formation of the feldspar 
nor the presence of trace quantities of galena 
could have modified the isotopic composition 
of feldspar lead. The unwashed feldspar is inac- 
tive when radiometrically assayed, but, consid- 
ering the statistical fluctuations in the back- 
sound, it is computed that the sample could 
(but need not necessarily) contain a maximum 
of 1 ppm of U or 3 ppm of Th. Thus, even if 
itis assumed that the sample is radioactive and 
the activity is due to thorium, the concentration 
of thorium is still about 15 times less than that 
required to account for the excess of Pb™, 

Data on the isotopic composition of lead 
isolated from three feldspars, drawn from 
three distinct environments, are available. 
Tilton ef al. (1955) reported that, although 
the feldspar component of the granite they 
analyzed has an abnormal excess of radiogenic 
Pb™* and Pb*, the feldspar in the pegmatite 
which intruded the granite has an isotopic com- 
position which corresponds to that of “normal” 
ore lead of the same age. Patterson (unpub- 
lished) found ‘‘normal’ lead in the Southern 
California batholith granite of Cretaceous age. 
The lead in the feldspar under study, with its 
deficiency of Pb and excess of Pb’, ap- 
parently constitutes a new and distinct type. 

The age data, given earlier, permit the 
following explanation: Both samarskite and 
ieldspar crystallized 580 m.y. ago. The feldspar, 
however, drew into its lattice lead evolved in a 
region which, for long periods in the past, con- 
tained abnormally high concentrations of 
thorium and low concentration of uranium rela- 
tive to lead. It remained “frozen” with this 
abnormal lead as there was practically no 
subsequent contribution of lead from the 
negligible amounts of U and/or Th contained 
in the feldspar. 

The age of the samarskite indicates that the 
Kishengarh pegmatite is very young Pre- 
cambrian or earliest Paleozoic; it is younger 
than the pegmatites of Bisundni and Ajmer- 
Merwara dated by Holmes ef al. (1949) and 
Holmes (1955). This finding is consistent with 
geological evidence as it is not necessary to 
suppose that pegmatites intruded at the same 
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time all along the immense synclinorium 


(Heron, 1953, p. 368). 
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